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ABSTRACT 

We review progress over the past decade in observations of large-scale star formation, with a focus 
on the interface between extragalactic and Galactic studies. Methods of measuring gas contents and 
star formation rates are discussed, and updated prescriptions for calculating star formation rates are 
provided. We review relations between star formation and gas on scales ranging from entire galaxies 
to individual molecular clouds. 

Subject headings: star formation, galaxies, Milky Way 



1. OVERVIEW 
1.1. Introduction 

Star formation encompasses the origins of stars and 
planetary systems, but it is also a principal agent of 
galaxy formation and evolution, and hence a subject at 
the roots of astrophysics on its largest scales. 

The past decade has witnessed an unprecedented 
stream of new observational information on star forma- 
tion on all scales, thanks in no small part to new facil- 
ities such as the Galaxy Evolution Explorer (GALEX), 
the Spitzer Space Telescope, the Herschel Space Obser- 
vatory, the introduction of powerful new instruments on 
the Hubble Space Telescope (HST), and a host of ground- 
based optical, infrared, submillimeter, and radio tele- 
scopes. These new observations are providing a detailed 
reconstruction of the key evolutionary phases and phys- 
ical processes that lead to the formation of individual 
stars in interstellar clouds, while at the same time ex- 
tending the reach of integrated measurements of star for- 
mation rates (SFRs) to the most distant galaxies known. 
The new data have also stimulated a parallel renaissance 
in theoretical investigation and numerical modelling of 
the star formation process, on scales ranging from in- 
dividual protostellar and protoplanetary systems to the 
scales of molecular clouds and star clusters, entire galax- 
ies and ensembles of galaxies, even to the first objects, 
which are thought to have reionized the Universe and 
seeded today's stellar populations and Hubble sequence 
of galaxies. 

This immense expansion of the subject, both in terms 
of the volume of results and the range of physical scales 
explored, may help to explain one of its idiosyncracies, 
namely the relative isolation between the community 
studying individual star-forming regions and stars in the 
Milky Way (usually abbreviated hereafter as MW, but 
sometimes referred to as "the Galaxy"), and the largely 
extragalactic community that attempts to characterise 
the star formation process on galactic and cosmologi- 
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cal scales. Some aspects of this separation have been 
understandable. The key physical processes that deter- 
mine how molecular clouds contract and fragment into 
clumps and cores and finally clusters and individual stars 
can be probed up close only in the Galaxy, and much of 
the progress in this subject has come from in-depth case 
studies of individual star-forming regions. Such detailed 
observations have been impossible to obtain for even rel- 
atively nearby galaxies. Instead the extragalactic branch 
of the subject has focused on the collective effects of star 
formation, integrated over entire star-forming regions, or 
often over entire galaxies. It is this collective conversion 
of baryons from interstellar gas to stars and the emer- 
gent radiation and mechanical energy from the stellar 
populations that is most relevant to the formation and 
evolution of galaxies. As a result, much of our empirical 
knowledge of star formation on these scales consists of 
scaling laws and other parametric descriptions, in place 
of a rigorous, physically-based characterization. Improv- 
ing our knowledge of large-scale star formation and its 
attendant feedback processes is essential to understand- 
ing the birth and evolution of galaxies. 

Over the past several years, it has become increasingly 
clear that many of the key processes influencing star for- 
mation on all of these scales lie at the interface between 
the scales within individual molecular clouds and those of 
galactic disks. It is now clear that the large scale SFR is 
determined by a hierarchy of physical processes spanning 
a vast range of physical scales: the accretion of gas onto 
disks from satellite objects and the intergalactic medium 
(Mpc) ; the cooling of this gas to form a cool neutral phase 
(kpc); the formation of molecular clouds 10 — lOOpc); 
the fragmentation and accretion of this molecular gas to 
form progressively denser structures such as clumps ('^1 
pc) and cores ('^ 0.1 pc); and the subsequent contraction 
of the cores to form stars (R©) and planets AU). The 
first and last of these processes operate on galactic (or 
extragalactic) and local cloud scales, respectively, but 
the others occur at the boundaries between these scales, 
and the coupling between processes is not yet well under- 
stood. Indeed it is possible that different physical pro- 
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cesses provide the "critical path" to star formation in dif- 
ferent interstellar and galactic environments. Whatever 
the answer to these challenging questions, however, Na- 
ture is strongly signalling that we need a unified approach 
to understanding star formation, one which incorporates 
observational and astrophysical constraints on star for- 
mation efficiencies, mass functions, etc. from small-scale 
studies along with a much deeper understanding of the 
processes that trigger and regulate the formation of star- 
forming clouds on galactic scales, which, in turn, set the 
boundary conditions for star formation within clouds. 

This review makes a modest attempt to present a con- 
solidated view of large-scale star formation, one which in- 
corporates our new understanding of the star formation 
within clouds and the ensemble properties of star forma- 
tion in our own MW into our more limited but broader 
understanding of star formation in external galaxies. As 
will become clear, the subject itself is growing and trans- 
forming rapidly, so our goal is to present a progress re- 
port in what remains an exciting but relatively immature 
field. Nevertheless a number of factors make this a timely 
occasion for such a review. The advent of powerful multi- 
wavelength observations has transformed this subject in 
fundamental ways since the last large observational re- 
view of the galaxy-scale aspects (Kennicutt 1998a), here- 
after denoted K98. Reviews have covered various obser- 
vational aspects of star formation within the MW (Evans 
1999; Bergin & Tafalla 2007; Zinnecker & Yorke 2007), 
but the most comprehensive reviews have been theoret- 
ically based (Shu et al. 1987; McKee & Ostriker 2007). 
We will take an observational perspective, with empha- 
sis on the interface between local and galactic scales. We 
focus on nearby galaxies and the MW, bringing in re- 
sults from more distant galaxies only as they bear on the 
issues under discussion. 

The remainder of this article is organized as follows. 
In the next subsection, we list some definitions and con- 
ventions that will be used throughout the paper. All 
observations in this subject rest on quantitative diagnos- 
tics of gas properties and star formation rates on various 
physical scales, and we review the current state of these 
diagnostics (§2, §3). In §4, we review those properties of 
Galactic star-forming regions which are most relevant for 
comparison to other galaxy-wide studies. In §5, we re- 
view the star-forming properties of galaxies on the large 
scale, including the MW; much of this section is effec- 
tively an update of the more extended review presented 
in K98. Section 6 updates our knowledge of relations be- 
tween star formation and gas (e.g., the Schmidt law) and 
local tests of these; relations. The review concludes in §7 
with our attempted synthesis of what has been learned 
from the confiuence of local and global studies and a look 
ahead to future prospects. 

We list in §1.3 some key questions in the field. Some 
have been or will be addressed by other reviews. Some 
are best answered by observations of other galaxies; some 
can only be addressed by observations of local star- 
forming regions in the MW. For each question we list 
sections in this review where we review the progress to 
date in answering them, but many remain largely unan- 
swered and are referred to the last section of this review 
on future prospects. 

Space limitations prevent us from citing even a fair 
fraction of the important papers in this subject, so we 



instead cite useful examples and refer the reader to the 
richer lists of papers that are cited there. 

1.2. Definitions and Conventions 

Here we define some terms and symbols that will be 
used throughout the review. 

The term "cloud" refers to a structure in the interstel- 
lar medium (ISM) separated from its surroundings by the 
rapid change of some property, such as pressure, surface 
density, or chemical state. Clouds have complex struc- 
ture, but theorists have identified two relevant structures: 
clumps are the birthplaces of clusters; cores are the birth- 
places of individual or binary stars (e.g., McKee & Os- 
triker 2007). The observational equivalents have been 
discussed (e.g., Table 1 in Bergin & Tafalla 2007), and 
cores reviewed (di Francesco et al. 2007; Ward-Thompson 
et al. 2007a), and we discuss them further in §2.2. 

We will generally explicitly use "surface density" or 
"volume density" , but symbols with S will refer to sur- 
face density, A'^ will refer to column density, and symbols 
with n ov p will refer to number or mass volume den- 
sity, respectively. The surface density of H I gas is rep- 
resented by Shi if He is not included and by S (atomic) 
if He is included. The surface density of molecular gas 
is symbolized by Ejj^ , which generally does not include 
He, or Emo/, which generally does include He. For MW 
observations, He is almost always included in mass and 
density estimates, and Timoi may be determined by ex- 
tinction or emission by dust (§2.3) or by observations 
of CO isotopologues (§2.4). When applied to an in- 
dividual cloud or clump, we refer to it as S(cloud) or 
E (clump). For extragalactic work, Ejj^ or S„io( are al- 
most always derived from CO observations, and inclu- 
sion of He is less universal. In the extragalactic context, 
^moi is best interpreted as related to a filling factor of 
molecular clouds until Y,moi ~ 100 M© pc~^, where the 
area filling factor of molecular gas may become unity 
(§2.4). Above this point, the meaning may change sub- 
stantially. The total gas surface density, with He, is 
Egos = Emoi + M^Hi, where jjL is the mean molecular 
weight per H atom; = 1.41 for MW abundances (Kauff- 
mann et al. 2008), though a value of 1.36 is commonly 
used. 

The conversion from CO intensity (usually in the 
J = 1 — > rotational transition) to column density of H2, 
not including He, is denoted A (CO), and the conversion 
from CO luminosity to mass, which usually includes the 
mass of He, is denoted aco • 

The term "dense gas" refers generally to gas above 
some threshold of surface density, as determined by ex- 
tinction or dust continuum emission, or of volume den- 
sity, as indicated by emission from molecular lines that 
trace densities higher than does CO. While not precisely 
defined, suggestions include criteria of a threshold surface 
density, such as T^moi > 125 Mq pc^^ (Goldsmith et al. 
2008; Lada et al. 2010; Heiderman et al. 2010), a volume 
density criterion (typically n > 10"* cm~^) (e.g., Lada 
1992), and detection of a line from certain molecules, 
such as HCN. The surface and volume density criteria 
roughly agree in nearby clouds (Lada et al. 2012) but 
may not in other environments. 

The star formation rate is symbolized by SFR or M*, 
often with units of Mq yr~ or M0 Myr~^, and its 
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meaning depends on how miich averaging over time and 
space is involved. The surface density of star forma- 
tion rate, E(SFR), has the same averaging issues. The 
efRciency of star formation is symbohzed by e when it 
means Mjf/(M* + Mdoud)-, as it usually does in Galac- 
tic studies. For extragalactic studies, one usually means 
Mt,/ Mgas, which we symbolize by e'. The depletion time, 
tdep = Another time often used is the "dynamical" 
time, tdyn, which can refer to the free-fall time (t//), the 
crossing time {tcross)^ or the galaxy orbital time {torb)- 
In recent years, it also has become common to compare 
galaxies in terms of their SFR per unit galaxy mass, or 
"specific star formation rate" (SSFR). The SSFR scales 
directly with the stellar birthrate parameter 6, the ra- 
tio of the SFR today to the average past SFR over the 
age of the galaxy, and thus provides a useful means for 
characterizing the star formation history of a galaxy. 

The far-infrared himinosity has a number of defini- 
tions, and consistency is important in converting them 
to M*. A commonly used definition integrates the dust 
emission over the wavelength range 3-1100 /xm (Dale & 
Helou 2002), and following those authors we refer to this 
as the total-infrared or TIR luminosity. Note however 
that other definitions based on a narrower wavelength 
band or even single-band infrared measurements are of- 
ten used in the literature. 

We refer to mass functions with the following conven- 
tions. Mass functions are often approximated by power 
laws in logM: 

dN/dlogM oc M-T, (1) 

In particular, the initial mass function (IMF) of stars is 
often expressed by equation 1 with a subscript to indicate 
stars (e.g., M*). 

Structures in gas are usually characterized by distribu- 
tions versus mass, rather than logarithmic mass: 

dN/dM oc M-", (2) 

for which a = 'y + 1. Terminology in this area is wildly 
inconsistent. To add to the confusion, some references 
use cumulative distributions, which decrease the index 
in a power law by 1. 

The luminosities and masses of galaxies are well repre- 
sented by an exponentially truncated power-law function 
(Schechter 1976): 

$(i) = ($7i*)(L/L*)"e(--^/-^*) (3) 

The parameter a denotes the slope of the power-law func- 
tion at low luminosities, and L* represents the luminos- 
ity above which the number of galaxies declines sharply. 
For the relatively shallow slopes (a) which are typical of 
present-day galaxies, L* also coincides roughly with the 
peak contribution to the total light of the galaxy popu- 
lation. As discussed in §5.2, the Schechter function also 
provides a good fit to the distribution function of total 
SFRs of galaxies, with the characteristic SFR in the ex- 
ponential designated as SFR*, in analogy to L* above. 

The term "starburst galaxy" was introduced by Weed- 
man et al. (1981), but nowadays the term is applied to a 
diverse array of galaxy populations. The common prop- 
erty of the present-day populations of starbursts is a SFR 



out of equilibrium, much higher than the long-term aver- 
age SFR of the system. No universally accepted quanti- 
tative definition exists, however. Some of the more com- 
monly applied criteria are SFRs that cannot be sustained 
for longer than a small fraction (e.g., < 10%) of the 
Hubble time, i.e., with gas consumption timescales of 
less than ^1 Gyr, or galaxies with disk-averaged SFR 
surface densities S(SFR) > 0.1 M© yr-ikpc"2(§5.2). 
Throughout this review we will use the term "quiescent 
star-forming galaxies" merely to characterize the non- 
starbursting galaxy population. Note that these local cri- 
teria for identifying starbursts are not particularly usefiil 
for high-redshift galaxies; a young galaxy forming stars 
at a constant SFR might more resemble a present-day 
starburst galaxy than a normal galaxy today. 

1.3. Questions 

1. How should we interpret observations of the main 
molecular diagnostic lines (e.g., CO, HCN) and 

millimeter-wave dust emission? How does this in- 
terpretation change as a function of metallicity, 
surface density, location within a galaxy, and star 
formation environment? (§2.4, §5.1, §7.3) 

2. How does the structure of the ISM, the structure 
of star-forming clouds, and the star formation itself 
change function of metallicity, surface density, 
location within a galaxy, and star formation envi- 
ronment? (§2.4, §5, §7.3). 

3. How do the mass spectra of molecular clouds and 

dense clumps in clouds vary between galaxies and 
within a galaxy? (§2.5, §5.1, §7.3) 

4. How constant is the IMF, and how are star for- 
mation rate measurements affected by possible 
changes in the IMF or by incomplete sampling of 
the IMF? (§2.5, §3.3, §4.2, §6.4) 

5. What are the limits of applicability of current star 
formation rate tracers? How are current mea- 
surements biased by dust attenuation or the ab- 
sence of dust, and how accurately can the effects of 

dust be removed? How do different tracers depend 
on metallicity, and what stellar mass ranges and 
timescales do they probe? (§3) 

6. How long do molecular clouds live and how can we 

best measure lifetimes? Do these lifetimes change 
systematically as functions of cloud mass, location 
in a galaxy, or some other parameter? (§4.3, §7.3) 

7. Do local observations provide any evidence for bi- 

modality in modes of star formation, for example 
distributed versus clustered, low-mass versus high- 
mass star formation? (§4, §7.3) 

8. On the scale of molecular clouds, what are the 
star formation efficiencies [e = M*/ (M* -|- Mdoud)], 
and star formation rates per unit mass (e' = 
M*/ Mdoud), and do these efficiencies vary system- 
atically as functions of cloud mass or other param- 
eters? (§4, §7.3) 
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9. How do spatial sampling and averaging affect the 
observed form of the star formation relations ex- 
pressed in terms of total, molecular, or dense gas 
surface densities? How well do M* diagnostics used 
in extragalactic studies work on finer scales within 
a galaxy? (§3.9, §6, §6.4, §7.3) 

10. Are there breakpoints or thresholds where either 
tracers or star formation change their character? 

(§2, §3, §7) 

2. ISM DIAGNOSTICS: GAS AND DUST 
2.1. ISM Basics 

The interstellar medium (ISM) is a complex envi- 
ronment encompassing a very wide range of properties 
(e.g., McKee & Ostriker 1977; Cox 2005). In the Milky 
Way, about half the volume is filled by a hot ionized 
medium (HIM) with n < 0.01 cm'^ and Tk > 10^ K, 
and half is filled by some combination of warm ionized 
medium (WIM) and warm neutral medium (WNM), with 
n ~ 0.1 — 1 cm~'^ and Tk of several thousand K (Cox 
2005). The neutral gas is subject to a thermal instabil- 
ity that predicts segregation into the warm neutral and 
cold neutral media (CNM), the latter with n > 10 cm~^ 
and Tk < 100 K (Field et al. 1969). The existence of 
the CNM is clearly revealed by emission in the hyper- 
fine transition of H I, with median temperature about 70 
K, determined by comparison of emission to absorption 
spectra (Heiles & Troland 2003). However, the strict 
segregation into cold and warm stable phases is not sup- 
ported by observations suggesting that about 48% of the 
WNM may lie in the thermally unstable phase (Heiles & 
Troland 2003). Finally, the molecular clouds represent 
the coldest {Tk ~ 10 K in the absence of star forma- 
tion), densest (n > 30 cm~^) part of the ISM. We fo- 
cus on the cold {Tk < 100 K) phases here, as they are 
the only plausible sites of star formation. The warmer 
phases provide the raw material for the cold phases, per- 
haps through "cold" flows from the intergalactic medium 
(Dekel & Birnboim 2006). In the language of this review, 
such flows would be "warm." 

In the next section, we introduce some terminology 
used for the molecular gas as a precursor to the following 
discussion of the methods for tracing it. We then discuss 
the use of dust as a tracer (§2.3) before discussing the 
gas tracers themselves (§2.4). 

2.2. Clouds and Their Structures 

The CNM defined in the preceding section presents a 
wide array of structures, with a rich nomenclature, but 
we focus here on the objects referred to as clouds, invok- 
ing a kind of condensation process. This name is perhaps 
even more appropriate to the molecular gas, which cor- 
responds to a change in chemical state. 

The molecular cloud boundary is usually defined by 
the detection, above some threshold, of emission from 
the lower rotational transitions of CO. Alternatively, a 
certain level of extinction of background stars is often 
used. In confused regions, velocity coherence can be 
used to separate clouds at different distances along the 
line of sight. However, molecular clouds are surrounded 
by atomic envelopes and a transition region in which 
the hydrogen is primarily molecular, but the carbon is 



mostly atomic (van Dishoeck & Black 1988). These 
regions are known as PDRs, which stands either for 
Photo-Dissociation Regions, or Photon-Dominated Re- 
gions (HoUenbach & Tielens 1997). More recently, they 
have been referred to as "dark gas." 

The structure of molecular clouds is complex, leading 
to considerable nomenclatural chaos. Cloud projected 
boundaries, defined by either dust or CO emission or 
extinction, can be characterized by a fractal dimension 
around 1.5 (Falgarone et al. 1991; Scalo 1990; Stutzki 
et al. 1998), suggesting an intrinsic 3-D dimension of 2.4. 
When mapped with sufficient sensitivity and dynamic 
range, clouds are highly structured, with filaments the 
dominant morphological theme; denser, rounder cores 
are often found within the filaments (e.g., Andre et al. 
2010; Men'shchikov et al. 2010; Molinari et al. 2010). 

The probability distribution function of surface den- 
sities, as mapped by extinction, can be fitted to a log- 
normal function for low extinctions. Ay < 2 — 5 mag 
(Lombardi et al. 2010), with a power law tail to higher 
extinctions, at least in some clouds. In fact, it is the 
clouds with the power-law tails that have active star for- 
mation (Kainulainen et al. 2009). Further studies in- 
dicate that the power-law tail begins at 'Smoi 40 — 80 
M© pc^^, and that the material in the power- law tail lies 
in imboimd clumps with mean density about 10^ cm~^ 
(Kainulainen et al. 2011). 

Within large clouds, and especially within the fila- 
ments, there may be small cores, which are much less 
massive but much denser. Theoretically, these are the 
sites of formation of individual stars, either single or 
binary (Williams et al. 2000). Their properties have 
been reviewed by di Francesco et al. (2007) and Ward- 
Thompson et al. (2007a). Those cores without evidence 
of ongoing star formation are called starless cores, and 
the subset of those that are gravitationally bound are 
called prestellar cores. Distinguishing these two requires 
a good mass estimate, independent of the virial theorem, 
and observations of lines that trace the kinematics. A 
sharp decrease in turbulence may provide an alternative 
way to distinguish a prestellar core from its surroundings 
(Pineda et al. 2010a). 

While clouds and cores are reasonably well defined, 
intermediate structures are more problematic. Theoreti- 
cally, one needs to name the entity that forms a cluster. 
Since large clouds can form multiple (or no) clusters, an- 
other name is needed. Williams et al. (2000) established 
the name "clump" for this structure. The premise was 
that this region should also be gravitationally bound, at 
least until the stars dissipate the remaining gas. Finding 
an observable correlative has been a problem (§2.5). 

Finally, one must note that this hierarchy of structure 
is based on relatively large clouds. There are also many 
small clouds (Clemens & Barvainis 1988), some of which 
are their own clumps and/or cores. The suggestion by 
Bok & Reilly (1947) that these could be the sites of star 
formation was seminal. 

2.3. Dust as an ISM Tracer 

In the Milky Way, roughly 1% of the ISM is found in 
solid form, primarily silicates and carbonaceous material 
(Draine 2003) (see also a very useful website^). Dust 

^ http://www.astro.princeton.edu/~draine/dust/dust.html 
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grain sizes range from 0.35 nm to around 1 /im, with ev- 
idence for grain growth in molecular clouds. The column 
density of dust can be measured by studying cither the 
reddening (using colors) or the extinction of starlight, 
using star counts referenced to an unobscured field. The 
relation between extinction (e.g., Ay) and reddening 
[e.g., E{B — V)] is generally parameterized by Ry = 
Ay I E{B ~V). In the diffuse medium, the reddening has 
been calibrated against atomic and molecular gas to ob- 
tain Ar(H I)+2Ar(H2) = 5.8 X 10^1 cm-2£;(B-F)(mag), 
not including helium (Bohlin et al. 1978). For Ry = 3.1, 
this relation translates to N{R I)+27V(H2) = 1.87 x 10^^ 
cm~^ Ay (msig). This relation depends on mctallicity, 
with coefficients larger by factors of about 4 for the LMC 
and 9 for the SMC (Weingartner & Draine 2001). 

The relation derived for diffuse gas is often applied in 
molecular clouds, but grain growth flattens the extinc- 
tion curve, leading to larger values of Ry, and a much 
flatter curve in the infrared (Flaherty et al. 2007; Chap- 
man et al. 2009). More appropriate conversions can be 
found in the curves for Ry = 5.5 in the website refer- 
enced above. There is evidence of continued grain growth 
in denser regions, which particularly affects the opacity 
at longer wavelengths, such as the mm/submm regions 
(Shirley et al. 2011). 

Reddening and extinction of starlight provided early 
evidence of the existence of interstellar clouds (Barnard 
1908). Lynds (1962) provided the mostly widely used 
catalog of "dark clouds" based on star counts, up to 
Ay ~ 6 mag, at which point the clouds were effectively 
opaque at visible wavelengths. The availability of in- 
frared photometry for large numbers of background stars 
from 2MASS, ISO, and Spitzer towsLvd nearby clouds has 
allowed extinction measurements up to Ay ^ 40 mag. 
These have proved very useful in characterizing cloud 
structure and in checking the gas tracers discussed in 
§2.4. Extinction mapping with sufficiently high spatial 
resolution can identify cores (e.g., Alves et al. 2007), but 
many axe not gravitationally bound (Lada et al. 2008), 
and they are likely to dissipate. 

At wavelengths with strong, diffuse Galactic back- 
ground emission, such as the mid-infrared, dust absorp- 
tion against the continuum can be used in a way paral- 
lel to extinction maps (e.g., Bacmann et al. 2000; Stutz 
et al. 2007). The main uncertainties in deriving the gas 
column density are the dust opacity in the mid-infrared 
and the fraction of foreground emission. 

The energy absorbed by dust at short wavelengths is 
emitted at longer wavelengths, from infrared to millime- 
ter, with a spectrum determined by the grain temper- 
ature and the opacity as a function of wavelength. At 
wavelengths around 1 mm, the emission is almost always 
optically thin and close to the Raylcigh- Jeans limit for 
modest temperatures. (For Td — 20K, masses, densities, 
etc. derived in the Rayleigh- Jeans limit at A = 1 mm 
must be multiplied by 1.5). Thus, the millimeter-wave 
continuum emission provides a good tracer of dust (and 
gas, with knowledge of the dust opacity) column density 
and mass. 

As a by-product of studies of the CMB, CORE, 
WMAP, and Planck produce large-scale maps of 
mm/submm emission from the Milky Way. With Planck, 
the angular resolution reaches 6 = 4.7' at 1.4 mm in 



an all-sky map of dust temperature and optical depth 
(Planck Collaboration et al. 2011b). Mean dust temper- 
ature increases from = 14 — 15 K in the outer Galaxy 
to 19 K aroimd the Galactic center. A catalog of cold 
clumps with = 7 — 17 K, mostly within 2 kpc of the 
Sun, many with surprisingly low column densities, also 
resulted from analysis of Planck data (Planck Collabo- 
ration et al. 2011c). Herschel surveys will be providing 
higher resolution [25" at 350 fim (Griffin et al. 2010) and 
6" at 70 /um (Poglitsch et al. 2010)] images of the nearby 
clouds (Andre et al. 2010), the Galactic Plane (Molinari 
et al. 2010), and many other galaxies (e.g., Skibba et al. 
2011; Wuyts et al. 2011). Still higher resolution is avail- 
able at A > 350 /im from ground-based telescopes at 
dry sites. Galactic plane surveys at wavelengths near 1 
mm have been carried out in both hemispheres (Aguirre 
et al. 2011; Schuller et al. 2009). Once these are put on 
a common footing, we should have maps of the Galac- 
tic plane with resolution of 20" to 30". Deeper sur- 
veys of nearby clouds (Ward- Thompson et al. 2007b), 
the Galactic Plane, and extragalactic fields are planned 
with SCUBA-2. 

While studies from ground-based telescopes offer 
higher resolution, the need to remove atmospheric fluc- 
tuations leads to loss of sensitivity to structures larger 
than typically about half the array footprint on the sky 
(Aguirre et al. 2011). While the large-scale emission is 
lost, the combination of sensitivity to a minimum column 
density and a maximum size makes these surveys effec- 
tively select structures of a certain mean volume density. 
In nearby clouds, these correspond to dense cores (Enoch 
et al. 2009). For Galactic plane surveys, sources mainly 
correspond to clumps but can be whole clouds at the 
largest distances (Dunham et al. 2011b). 

2.4. Gas Tracers 

One troublesome aspect of extragalactic studies has 
been the determination of the amount and properties of 
the gas. We discuss flrst the atomic phase, then the 
molecular phase, and finally tracers of dense gas. 

2.4.1. The Cold Atomic Phase 

The cold, neutral atomic phase is traced by the hyper- 
fine transition of hydrogen, occuring at 21 cm in the rest 
frame. This transition reaches optical depth of unity at 

7V(H I) = 4.57 X 1020cm-2(T,p,„/100K)(cr„/lkm s'^), 

(4) 

where Tgpin is the excitation temperature, usually equal 
to the kinetic temperature, Tk, and is the velocity 
dispersion (eq. 8.11 in Draine 2011). For dust in the 
diffuse ISM of the Milky Way, this column density cor- 
responds to Ay = 0.24 mag, and optical depth effects 
can be quite important. The warm, neutral phase is dif- 
ficult to study, as lines become broad and weak, but it 
has been detected (Kulkarni & Heiles 1988). A detailed 
analysis of both CNM and WNM can be found in Heiles 
& Troland (2003). 
With /(H I) / S^dv (Jy km s-^), 

M(H I)/M0 = 2.343 x 10^(1 + z)-^{DL/M-pc)^I{}i I), 

(5) 

not including helium (eq. 8.21 in Draine 2011, but note 
errata regarding {1 + z) factor, which was in the numer- 
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ator in the 2011 edition). A coefficient of 2.356 x 10'^ is 
in common use among H I observers (e.g., Walter et al. 
2008; Catinella et al. 2010; Giovanelli et al. 2005; Meyer 
et al. 2004), usually traced to Roberts (1975). The dif- 
ference is almost certainly caused by a newer value for 
the Einstein A value in Draine (2011). 

2.4.2. Molecular Gas 

The most abundant molecule is H2, but its spectrum 
is not a good tracer of the mass in molecular clouds. 
The primary reason is not, as often stated, the fact that 
it lacks a dipole moment, but instead the low mass of 
the molecule. The effect of low mass is that the ro- 
tational transitions require quite high temperatures to 
excite, making the bulk of the H2 in typical clouds invis- 
ible. Continuum optical depth at the mid-infrared wave- 
lengths of these lines is also an impediment. The ob- 
served rotational H2 lines originate in surfaces of clouds, 
probing 1% to 30% of the gas in a survey of galaxies (e.g., 
Roussel et al. 2007). 

Carbon monoxide (CO) is the most commonly used 
tracer of molecular gas because its lines are the strongest 
and therefore easiest to observe. These advantages are 
accompanied by various drawbacks, which MW observa- 
tions may help to illuminate. Essentially, CO emission 
traces column density of molecular gas only over a very 
limited dynamic range. At the low end, CO requires pro- 
tection from photodissociation and some minimum den- 
sity to excite it, so it does not trace all the molecular 
(in the sense that hydrogen is in H2) gas, especially if 
metallicity is low. At the high end, emission from CO 
saturates at modest column densities. A more complete 
review by A. Bolatto of the use of CO to trace molec- 
ular gas is scheduled for a future ARAA, so we briefly 
highlight the issues here. 

CO bccomc^s optically thick at very low total column 
densities. Using CO intensity to estimate the column 
density of a cloud is akin to using the presence of a brick 
wall to estimate the depth of the building behind it. The 
isotopologues of CO (^'^CO and C^'^O) provide lower op- 
tical depths, but at the cost of weaker lines. To use these, 
one needs to know the isotope ratios; these are fairly well 
known in the Milky Way (Wilson & Rood 1994), but they 
may differ substantially in other galaxies. In addition, 
the commonly observed lower transitions of CO are very 
easily thermalized {T^x = Tk) making them insensitive 
to the presence of dense gas. 

The common practice in stiidies of MW clouds is to 
use a combination of CO and ^"'CO to estimate optical 
depth and hence CO column density (A^(CO)). One can 
then correlate A^(CO) or its isotopes against extinction 
to determine a conversion factor (Dickman 1978; Frerk- 
ing et al. 1982). A recent, careful study of CO and ^^CO 
toward the Taurus molecular cloud shows that A^(CO) 
traces Ay up to 10 mag in some regions but only up to 4 
mag in other regions, and other issues cause problems be- 
low Av = 3 mag (Pineda et al. 2010b). Freeze-out of CO 
and conversion to CO2 ice causes further complications 
in cold regions of high column density (Lee et al. 2003; 
Pontoppidan et al. 2008; Whittet et al. 2007), but other 
issues likely dominate the interpretational uncertainties 
for most extragalactic CO observations. 

Since studies of other galaxies generally have only CO 
observations available, the common practice has been to 



relate H2 column density (A^(H2)) to the integrated in- 
tensity of CO [/(CO) = '/ Tdv (K km s^^)] via the infa- 
mous "X-factor": N(ll2)= X{CO)I{CO). For example, 
Kennicutt (1998b) used X(CO) = 2.8 x 10^° cm'^ (K 

km s-i)-i. 

Pineda et al. (2010b) find X{CO) = 2.0 x lO^" cm'^ 
(K km s~^)~^ in the parts of the cloud where both CO 
and -'^^CO are detected along individual sight lines. CO 
does emit beyond the area where ^''CO is detectable and 
even, at a low level, beyond the usual detection thresh- 
old for CO. Stacking analysis of positions beyond the 
usual boundary of the Taurus cloud suggest a factor 
of two additional mass in this transition region (Gold- 
smith et al. 2008). Such low intensity, but large area, 
emission could contribute substantially to observations 
of distant regions in our galaxy and in other galaxies. 
For the outer parts of clouds, where only CO is de- 
tected, and only by stacking analysis, X{GO) is six times 
higher (Pineda et al. 2010b). Even though about a quar- 
ter of the total mass is in this outer region, they argue 
that the best single value to use is X{CO) = 2.3 x 10^° 
cm~^ (K km s~^)~^, similar to values commonly adopted 
in extragalactic studies. When possible, we will use 
X{CO) = 2.3 x 10^0 cm-2 (K km s-^)-^ for this paper. 
In contrast, Liszt et al. (2010) argue that A'(CO) stays 
the same even in diffuse gas. If they are right, much of 
the integrated CO emission from some parts of galaxies 
could arise in diffuse, unbound gas. 

Comparing to extinction measures in two nearby 
clouds extending up to Ay = 40 mag, Heidcrman et al. 
(2010) found considerable variations from region to re- 
gion, with /(CO) totally insensitive to increased extinc- 
tion in many regions. Nonetheless, when averaged over 
the whole cloud, the standard X{CO) would cause errors 
in mass estimates of ±50%, systematically overestimat- 
ing the mass for I^moi < 150 Mq pc~^ and underestimat- 
ing for E„joi > 200 M0 pc~^. Simulations of turbulent 
gas with molecule formation also suggest that /(CO) can 
be a very poor tracer of Ay (Glover et al. 2010). Overall, 
the picture at this point is that measuring the extent of 
the brick wall will not tell you anything about particu- 
larly deep or dense extensions of the building, but if the 
buildings are similar, the extent of the wall is a rough 
guide to the total mass of the building (cf. Shetty et al. 
2011b). 

Following that thought, the luminosity of CO 
J = 1 — ^ is often used as a measure of cloud mass, 
with the linewidth reflecting the virial theorem in a 
crude sense. In most extragalactic studies, individual 
clouds arc not resolved, and the luminosity of CO is 
essentially a cloud counting technique (Dickman et al. 
1986). Again, this idea works only if the objects be- 
ing counted have rather uniform properties. In this ap- 
proach, M = acoL{CO). For the standard value for 
X{CO) discussed above (X(CO) = 2.3 x 10^° cm'^ (K 
km s~^)~-^) aco = 3.6 Mq(K km s~^pc^)~-^ without he- 
lium, and aco = 4.6 M0(K km s~^pc^)~^ with correc- 
tion for helium. 

The effects of cloud temperature, density, metallicity, 
etc. on mass estimation from CO were discussed by Mal- 
oney & Black (1988), who concluded that large variations 
in the conversion factor are likely. Ignoring metallicity ef- 
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fects, they would predict X{CO) oc ttP-^Tk ^ , where n is 
the mean density in the cloud. Shetty et al. (2011b) find 
a weaker dependence on T^, X{GO) oc Tk~^'^ ■ These 
scalings rests on some arguments that may not apply in 
other galaxies. Further modeling with different metal- 
licities (Shetty et al. 2011a; Feldmann et al. 2012) has 
begun to provide some perspective on the range of be- 
haviors likely in other galaxies. 

There is observational evidence for changes in aco in 
centers of galaxies. In the MW central region, Oka et al. 
(1998) suggest a value for X{CO) a factor of 10 lower 
than in the MW disk. The value of aco appears to 
decrease by factors of 2-5 in nearby galaxy nuclei (e.g., 
Meier & Turner 2004; Meier et al. 2008), and by a factor 
of about 5-6 in local ULIRGs (Downes & Solomon 1998; 
Downes et al. 1993), and probably by a similar factor in 
high-redshift, molecule-rich galaxies (Solomon & Vanden 
Bout 2005). 

A compilation of measurements indicates that aco de- 
creases from the usually assumed value with increasing 
surface density of gas once Y^^oi > 100 pc~^ (Fig. 1, 
Tacconi et al. 2008). This Y,moi corresponds to Ay ~ 10 
mag, about the point where CO fails to trace column 
density in Galactic clouds, and the point where the beam 
on another galaxy might be filled with CO-emitting gas. 
Variation in aco has a direct bearing on interpretation 
of starbursts (§6), and it is highly unlikely that aco is 
simply bimodal over the full range of star forming en- 
vironments. Indeed, Narayanan et al. (2012) use a grid 
of model galaxies to infer a smooth function of /(CO): 
aco = min[6.3,10.7 x 7(CO)"°-^^]Z'-0-65, where Z' is 
the metallicity divided by the solar value. 

Other methods have been used to trace gas indi- 
rectly, including gamma-ray emission and dust emis- 
sion. Gamma-rays from decay of pions produced from 
high-energy cosmic rays interacting with baryons in the 
ISM directly trace all the matter if one knows the cos- 
mic ray flux (e.g., Bloemen 1989). Recent analysis of 
gamma-ray data from Fermi have inferred a low value 
for X{GO) of (0.87 ± 0.05) x 10^" cm^^ (K km s-^)-^ 
in the solar neighborhood (Abdo et al. 2010). This value 
is a factor of three lower than other estimates and early 
results from dust emission, using Planck, which found 
X{CO) = (2.54±0.13) X 1020 cm'^ (Kkms"!)"! for the 
solar neighborhood. (Planck Collaboration et al. 2011b). 

Abdo et al. (2010) also inferred the existence of "dark 
gas", which is not traced by either HI or CO, sur- 
rounding the CO-emitting region in the nearby Gould 
Belt clouds, with mass of about half of that traced by 
CO. Planck Collaboration et al. (2011b) have also indi- 
cated the existence of "dark gas" with mass 28% of the 
atomic gas and 118% of the CO-emitting gas in the solar 
neighborhood. The dark gas appears around Ay = 0.4 
mag, and is presumably the CO-poor, but H2-containing, 
outer parts of clouds (e.g., Wolfire et al. 2010), which can 
be seen in [CII] emission at 158 /xm and in fluorescent H2 
emission at 2 /xm (e.g., Pak et al. 1998). It is dark only 
if one restricts attention to H I and CO. At low metallic- 
ities, the layer of C+ grows, and a given CO luminosity 
implies a larger overall cloud (see Fig. 14 in Pak et al. 
1998). However, Planck Collaboration et al. (2011b) also 
suggest that the dark gas may include some gas that is 
primarily atomic, but not traced by H I emission, owing 



to optical depth effects. 

Dust emission is now being used as a tracer of gas for 
other galaxies. Maps of H I gas are used to break the 
degeneracy between the gas to dust ratio 5gdr and a^g 
in the following equation: 

Sgdr^dust = aco/(CO) + (6) 

Leroy et al. (2011) minimize the variation in 6gdr to find 
the best fit to a^Q in several nearby galaxies. They 
determine S<i„st from Spitzer images of 160 /um emis- 
sion, with the dust temperature (or equivalently the ex- 
ternal radiation field) constrained by shorter wavelength 
far-infrared emission. This method should measure ev- 
erything with substantial dust that is not H I, so in- 
cludes the so-called "dark molecular gas." They find that 
a^o = 3 — 9 Mq pc~2 (K km s~^)~^, similar to val- 
ues based on MW studies, for galaxies with metallici- 
ties down to 12-|-[0/H]~ 8.4, below which a'^Q increases 
sharply to values like 70. Their method does not address 
the column density within clouds, but only the average 
surface density on scales larger than clouds. For galaxies 
not too different from the MW, the standard a^g '^iil 
roughly give molecular gas masses, but in galaxies with 
lower metallicity, the effects can be large. An extreme 
is the SMC, where a'co ~ 220 M© pc-^ (K km s'^)-^ 
(Bolatto et al. 2011). Maps of dust emission at longer 
wavelengths should decrease the sensitivity to dust tem- 
perature. 

2.4.3. Dense Gas Tracers 

Lines other than CO J = 1 — > generally trace warmer 
(e.g., higher J CO lines) and/or denser (e.g., HCN, 
CS, ...) gas. Wu et al. (2010a) showed that the line 
luminosities of commonly used tracers (CS J = 2 ^ 1, 
J = 5 4, and J = 7 6; HCN J=l-)-0, J = 3-)-2) 
correlate well with the virial mass of dense gas; indeed 
most follow relations that are close to linear. Since these 
lines are optically thick, the linear correlation is some- 
what surprising, but it presumably has an explanation 
similar to the cloud-counting argument offered above for 
CO. Collectively, we refer to these lines as "dense gas 
tracers", but the effective densities increase with J and 
vary among species (see Evans 1999; Reiter et al. 2011b 
for definitions and tables of effective and critical densi- 
ties). As an example for HCN J = 1 — >■ 0, 

log(iHCNl-o) = 1-04 X log(Mvir(-RHCNl-o)) - 1-35 (7) 

from a robust estimation fit to data on dense clumps in 
the MW (Wu et al. 2010a). The line luminosities of other 
dense gas tracers also show strong correlations with both 
virial mass (Wu et al. 2010a) and the mass derived from 
dust continuum observations (Reiter et al. 2011b). 

These tracers of dense gas are of course subject to 
the same caveats about sensitivity to physical conditions, 
such as metallicity, volume density, turbulence, etc., as 
discussed above for CO. In fact, the other molecules 
are more sensitive to photodissociation than is CO, so 
they will be even more dependent on metallicity (e.g., 
Rosolowsky et al. 2011). Dense gas tracers should be- 
come much more widely used in the ALMA era, but cau- 
tion is needed to avoid the kind of over-simplification 
that has tarnished the reputation of CO. Observations 
of multiple transitions and rare isotopes, along with re- 
alistic models, will help (e.g., Garcia-Burillo et al. 2012). 
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2.4.4. Summary 

What do extragalactic observations of molecular trac- 
ers actually measure? None are actually tracing the 
surface density of a smooth medium, with the possi- 
ble exception of the most extreme starbursts. The lu- 
minosity of CO is a measure of the number of emit- 
ting structures times the mean luminosity per structure. 
Dense gas tracers work in a similar way, but select gas 
above a higher density threshold than does CO. For CO 
J = 1 — > 0, the structure is the CO emitting part of a 
cloud, which shrinks as metallicity decreases. For dense 
gas tracers, the structure is the part of a cloud dense 
enough to produce significant emission, something like a 
clump if conditions are like the MW, but even more sen- 
sitive to metallicity. Using a conversion factor, one gets 
a "mass" in those structures. 

Estimates of the conversion factors for CO vary by fac- 
tors of 3 for the solar neighborhood and at least an or- 
der of magnitude at low metallicity and by at least half 
an order of magnitude for our own Galactic Center and 
for extreme systems like ULIRGs. Dividing by the area 
of the galaxy or the beam, one gets "Emoj", which re- 
ally should be thought of as an area filling factor of the 
structures being probed times some crude estimate of 
the mass per structure. That estimate depends on con- 
ditions in the structures, such as density, temperature, 
and abundance and on the external radiation field and 
the metallicity. Once the L{GO) translates to surface 
densities above that of individual clouds {Timoi >~ 100 
M0 pc~^), the interpretation may change as the area fill- 
ing factor can now exceed unity. Large ranges of velocity 
in other galaxies (if not exactly face-on) allow /(CO) to 
still count clouds above this threshold, but the clouds 
may no longer be identical. The full range of inferred 
molecular surface densities in other galaxies inferred from 
CO, assuming a constant conversion factor, is a factor of 
10^ (§6); given that CO traces Y,moi only over a factor 
of 10 at best in well-studied clouds in the MW (Pineda 
et al. 2010b; Heiderman et al. 2010), this seems a dan- 
gerous extrapolation. 

2.5. Mass Functions of Stars, Clusters, and Gas 
Structures 

The mass functions of the structures in molecular 
clouds (§2.2) are of considerable interest in relation to 
the origin of the mass functions of clusters (clumps) and 
stars (cores). Salpeter (1955) fit the initial mass func- 
tion (IMF) of stars to a power law in logM* with ex- 
ponent 7* = 1.35, or a* — 2.35 (see §1.2 for defini- 
tions) . More recent determinations over a wider range of 
masses, including brown dwarfs, indicate a clear flatten- 
ing at lower masses, either as broken power laws (Kroupa 
et al. 1993) or a log-normal distribution (Miller & Scalo 
1979; Chabrier 2003). A detailed study of the IMF in the 
nearest massive yoking cluster, Orion, with stars from 0.1 
to 50 Mq with a mean age of 2 Myr, shows a peak be- 
tween 0.1 and 0.3 M©, depending on evolutionary mod- 
els, and a deflcit of brown dwarfs relative to the fleld IMF 
(Da Rio et al. 2012). Steeper IMFs in massive elliptical 
galaxies have been suggested by van Dokkum & Con- 
roy (2011). Zinnecker & Yorke (2007) summarize the 
evidence for a real (not statistical) cut-off in the IMF 
around 150 M© for star formation in the MW and Mag- 



ellanic Clouds. Variations in the IMF with environment 
are plausible, but convincing evidence for variation re- 
mains elusive (Bastian et al. 2010). 

The masses of embedded young clusters (Lada & 
Lada 2003), OB associations (McKee & WilUams 1997), 
and massive open clusters (Elmegreen & Efremov 1997; 
Zhang & Fall 1999), follow a similar relation with 
Jciuster ~ 1- Studies of clusters in nearby galaxies have 
also found Jcluster = 1 ± 0.1, with a possible upper mass 
cutoff or turn-over around 1 — 2 x 10^ Mq (Gieles et al. 
2006), perhaps dependent on the galaxy (but see Chan- 
dar et al. 2011). In contrast, the most massive known 
open clusters in the MW appear to have a mass of 
(6 - 8) X 10'' Mq (Portegies Zwart et al. 2010; Davies 
et al. 2011; Homeier & Alves 2005). 

The mass functions of clusters and stars are presum- 
ably related to the mass functions of their progenitors, 
clouds or clumps, and cores. For comparison to struc- 
tures in clouds, we will use the distributions versus mass, 
rather than logarithmic mass, so our points of compari- 
son will be a* = 2.3 and aduster = 2. From surveys of 
emission from CO J = 1 — >^ in the outer galaxy, where 
confusion is less problematic, Heyer et al. (2001) found 
a size distribution of clouds, dN/dr oc 7.-3-2±o.i^ with no 
sign of a break from the power-law form from 3 to 60 pc. 
The mass function, using the definition in equation 2, was 
fitted by adoud = 1.8 ±0.03 over the range 700 to 1 x 10^ 
Mq . Complementary surveys of the inner Galaxy found 
a-cioud = 1.5 up to M^ax = 10^-^ M© (Rosolowsky 2005). 
The upper cut-off around 3-6 x 10® M© appears to be 
real, despite issues of blending of clouds (Williams & 
McKee 1997; Rosolowsky 2005). Using ^^GO J =l^Q, 
which selects somewhat denser parts of clouds, Roman- 
Duval et al. (2010) found aia^o = ^''^^ ^ ^'^^ 
Mcioud > 10^ Mq. Clumps within clouds, identified 
by mapping ^^CO or C^^O, have a similar value for 
Oiciump = 1.4—1.8 (Kramer ct al. 1998). All results agree 
that most clouds are small, but, unlike stars or clusters, 
most mass is in the largest clouds (adoud < 2). 

Studies of molecular cloud properties in other galaxies 
have been recently reviewed (Fukui & Kawamura 2010; 
Blitz et al. 2007). Mass fimctions of clouds in nearby 
galaxies appear to show some differences, though sys- 
tematic effects introduce substantial uncertainty (Wong 
ct al. 2011; Sheth et al. 2008). Evidence for a higher 
value of adoud has been found for the LMC (Fukui ct al. 
2001; Wong et al. 2011) and M33 (Engargiola et al. 2003; 
Rosolowsky 2005). Most intriguingly, the mass function 
appears to be steeper in the interarm regions of M51 
than in the spiral arms (E. Schinnerer and A. Hughes, 
personal communication), possibly caused by the aggre- 
gation of clouds into larger structures within arms and 
disaggregation as they leave the arms (Koda et al. 2009). 

Much work has been done recently on the mass func- 
tion of cores, and Herschel surveys of nearby clouds will 
illuminate this topic (e.g., Andre ct al. 2010). At this 
point, it seems that the mass function of cores is clearly 
steeper than that of clouds, much closer to that of stars 
(Motte et al. 1998; Enoch et al. 2008; Sadavoy et al. 
2010). A turnover in the mass function appears at a 
mass about 3-4 times that seen in the stellar IMF (Alves 
et al. 2007; Enoch et al. 2008). The similarity of acore 
to a* supports the idea that the cores are the immediate 
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precursors of stars and that the mass function of stars 
is set by the mass function of cores. In this picture, the 
offset of the turnover implies that about 0.2 to 0.3 of 
the core mass winds up in the star, consistent with sim- 
ulations that include envelope clearing by outflows (e.g., 
Dunham et al. 2010). Various objections and caveats to 
this appealing picture have been registered (e.g.. Swift 
& Wilhams 2008; Reid et al. 2010; Clark ct al. 2007; 
Hatcheh & Fuller 2008). 

Note that no such tempting similarity exists between 
the mass function of clusters and the mass function of 
clumps defined by '^'^CO maps, questioning whether that 
observational definition properly fits the theoretical con- 
cept of a chimp as a cluster-forming entity. Structures 
marked by stronger emission in ^^CO or C^^O (Rath- 
borne et al. 2009) are not always clearly bound gravita- 
tionally. Ground and space-based imaging of submillime- 
ter emission by dust (e.g., Men'shchikov et al. 2010) offer 
promise in this area, but information on velocity struc- 
ture will also be needed. The data so far suggest that 
the structure is primarily filamentary, more like that of 
clouds than that of cores. CS, HCN and other tracers of 
much denser gas (§2.4), along with dust continuum emis- 
sion (§2.3) have identified what might be called "dense 
clumps" , which do appear to have a mass function simi- 
lar to that of clusters (Shirley et al. 2003; Beltran et al. 
2006; Reid & Wilson 2005), 

The comparison of mass functions supports the idea 
(McKee & Ostriker 2007) that cores are the progenitors 
of stars and dense clumps are the progenitors of clusters, 
with clouds less directly related. Upper limits to the 
mass function of clumps would then predict upper limits 
to the mass function of clusters, unless nearby clumps 
result in merging of clusters. Because most clouds and 
many clumps are more massive than the most massive 
stars, it is less obvious that an upper limit to stellar 
masses results from a limit on cloud or clump masses. If 
clump masses are limited and if the mass of the most mas- 
sive star to form depends causally (not just statistically) 
on the mass of the clump or cloud, the integrated galac- 
tic IMF (IGIMF) can be steepened (Kroupa & Weidner 
2003; Weidner & Kroupa 2006). We discuss the second 
requirement in §4.2 and consequences of this controver- 
sial idea in §3.3 and §6.4. 

3. SFR DIAGNOSTICS: THE IMPACT OF 
MULTI- WAVELENGTH OBSERVATIONS 

The influx of new observations over the past decade 
has led to major improvements in the calibration and 
validation of diagnostic methods for measuring SFRs in 
galaxies. Whereas measuring uncertainties of factors of 
two or larger in SFRs were commonplace ten years ago, 
new diagnostics based on multi-wavclcngth data arc re- 
ducing these internal uncertainties by up to an order of 
magnitude in many instances. These techniques have 
also reduced the impact of many systematic errors, in 
particular uncertainties due to dust attenuation, though 
others such as the IMF (§2.5) remain important limiting 
factors. 

A detailed discussion of SFR diagnostics and their cali- 
brations was given in K98, and here we highlight progress 
made since that review was published. In §3.8, we com- 
pile updated calibrations for the most commonly used 
indices, based on updated evolutionary synthesis models 



and IMF compared to K98. The new challenges which 
come with spatially-resolved measurements of galaxies 
are discussed in §3.9. 

3.1. Star Counting and CMD Analysis 

The most direct way to measure star formation rates is 
to count the number of identifiable stars of a certain age. 
Ideally, one would have reliable masses and ages for each 
star and the mean star formation rate would be given by 

(M,)= ^ N{M^,U)M,/U (8) 

M.=Mi 

where N{Mi,,tt,) is the number of stars in a mass bin 
characterized by mass and a lifetime bin (the time 

since formation) characterized by , and the star forma- 
tion rate would be averaged over the longest value of in 
the sum. In practice, complete information is not avail- 
able, and one needs to limit the allowed lifetimes to find 
the star formation rate over a certain period. In particu- 
lar, for nearby clouds in the MW, nearly complete lists of 
young stellar objects (YSOs) with infrared excesses are 
available. If one assumes that the YSOs sample the IMF 
in a typical way, one can derive the mean mass of stars 
((M*)), and the equation becomes 

(M*) = N{YSOs){M,)/te^,^ss- (9) 

With currently favored IMFs, (M,) = 0.5 Mq. The main 
source of uncertainty is in texcessi the duration of an 
infrared excess (§4.1). 

For young clusters one can determine mean agcis from 
fitting isochrones to a color-magnitude diagram, and 
measure M* assuming coeval formation as long as some 
stars have not yet reached the main sequence. Some clus- 
ters have measurements of stars down to very low masses 
(e.g., the Orion Nebula Cluster, Hillenbrand 1997), but 
most need corrections for unseen low- mass stars. For 
older clusters one can derive the total mass, but not (di- 
rectly) the duration of star formation, which has to be 
assumed. 

In principle similar techniques can be applied on a 
galaxy-wide basis. Beyond the Magellanic Clouds, cur- 
rent instruments cannot resolve individual YSOs in most 
regions, and as a result most studies of SFRs in galax- 
ies are based either on measurements of massive 0-type 
and/or Wolf-Rayet stars (see K98) or on measurements 
of the entire color- magnitude diagram (CMD). Consider- 
able progress has been made recently in using spatially- 
resolved mapping of CMDs to reconstruct spatially and 
age-resolved maps of star formation in nearby galaxies. 
Many of these use the distribution of blue helium-burning 
stars in the CMD (e.g., Dohm-Palmer et al. 2002), but 
more recent analyses fit the entire upper CMDs to syn- 
thetic stellar populations to derive estimates of spatially 
resolved stellar age distributions with formal uncertain- 
ties (e.g., Dolphin 2002; Weisz et al. 2008). This tech- 
nique does not have sufficient age resolution to determine 
accurate SFRs for ages less than ^10 Myr, but when ap- 
plied to high quality datasets such as those which can be 
obtained with HST, they can provide sufficient age reso- 
lution to constrain the temporal behaviors and changes in 
the spatial distributions of formation over the past 100 
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Myr or longer (e.g., Weisz at al. 2008; Williams et al. 
2011). 

3.2. Ultraviolet Continuum Measurements: The 
Impact of GALEX 

The near-ultraviolet emission of galaxies longward of 

the Lyman-continuum break directly traces the photo- 
spheric emission of young stars and hence is one of the 
most direct tracers of the recent SFR. For a conventional 
IMF, the peak contribution to the integrated UV lu- 
minosity of a young star cluster arises from stars with 
masses of several solar masses. Consequently this emis- 
sion traces stars formed over the past 10-200 Myr, with 
shorter timescales at the shortest wavelengths (see Table 
!)• 

For extragalactic studies, this subject has been revolu- 
tionised by the launch of the Galaxy Evolution Explorer 
(GALEX) mission (Martin et al. 2005a). It imaged ap- 
proximately two thirds of the sky in far-ultraviolet (FUV; 
155 nm) and near-ultraviolet (NUV; 230 nm) channels to 
limiting source fluxes tuab ~ 20.5, and obtained deeper 
full-orbit or multi-orbit imaging {niAB > 23) for selected 
galaxies and fields such as those in the Sloan Digital 
Sky Survey (SDSS). The spatial resolution (4.5" to 6" 
FWHM) of the imager makes it an especially powerful 
instrument for integrated measurements of distant galax- 
ies and resolved mapping of the nearest external galaxies. 
Although most scientific applications of GALEX data to 
date have arisen from its imaging surveys, a series of 
spectroscopic surveys and pointed observations of vary- 
ing depths were carried out as well (Martin et al. 2005a). 

The main impacts of GALEX for this subject are sum- 
marized in §5. Broadly speaking its largest impacts were 
in providing integrated UV fiuxes (and hence SFR es- 
timates) for hundreds of thousands of galaxies, and in 
exploiting the dark sky from space at these wavelengths 
to reveal star formation at low surface brightnesses and 
intensities across a wide range of galactic environments. 

Other spaceborne instruments have also provided im- 
portant datasets in this wavelength region, including the 
XMM Optical Monitor (Mason et al. 2001) and the Swift 
UV/Optical Telescope (Roming et al. 2005). Ahhough 
these instruments were primarily designed for follow-\ip 
of X-ray and gamma-ray observations, they also have 
been used to image nearby and distant galaxies, with 
the advantage of higher spatial resolution (~1" FWHM). 
Several important studies have also been published over 
the past decade from observations made with the Ultravi- 
olet Imaging Telescope on the ASTRO missions (Stecher 
et al. 1997; Marcum et al. 2001). The Hubble Space 
Telescope continues to be a steady source of observa- 
tions (mainly in the NUV) for targeted regions in nearby 
galaxies. 

The primary disadvantage of the ultraviolet is its se- 
vere sensitivity to interstellar dust attenuation. The 
availability of new data from GALEX and other facili- 
ties has stimulated a fresh look at this problem. If the 
intrinsic color of the emitting stellar population is known 
a priori^ the FUV-NUV color or the UV spectral slope 
(usually denoted /3) can be used to estimate the dust 
attenuation, and numerous calibrations have been pub- 
hshed (e.g., Calzetti et al. 1994; Kong et al. 2004; Seibert 
et al. 2005; Johnson et al. 2007; Salim et al. 2007; Treyer 
et al. 2007; Hao et al. 2011). The accuracy of these pre- 



scriptions rests heavily on the presumed (but uncertain) 
intrinsic colors, the shape of the dust extinction curve, 
and the complicated effects of geometry and scattering 
when averaging over a large physical region (e.g., Gordon 
et al. 2001). 

The abundance of high quality far-infrared observa- 
tions of nearby galaxies has made it possible to test how 
well the UV colors correlate with independent estimates 
of the dust attenuation from the IR/UV flux ratios. Ear- 
lier observations of bright starburst galaxies suggest a 
tight relation between the logarithmic IR/UV ratio (of- 
ten termed the "IRX") and UV color (e.g., Meurer et al. 
1999), and indeed this "IRX-/3 relation" provided the 
primary means for calibrating the color versus attenua- 
tion relation. Subsequent observations of a wider range 
of galaxies, however, has revealed a much larger scatter 
in the relation (e.g., Boquien et al. 2012). When galaxies 
of all types are considered, the scatter in actual FUV at- 
tenuations for a fixed FUV NUV color can easily be two 
orders of magnitude. Even when the comparison is re- 
stricted to galaxies with intrinsically high SFRs (as might 
be observed at high redshift), the uncertainties can easily 
reach an order of magnitude. As a result, more weight 
is being applied to alternative estimates of attenuation 
corrections and SFRs based on combinations of UV with 
IR measurements (§3.7). 

3.3. Emission-Line Tracers 

The remaining SFR indicators to be discussed here rely 
on measuring starlight that has been reprocessed by in- 
terstellar gas or dust, or on tracers related to the death 
of massive stars. The most widely applied of these are 
the optical and near-infrared emission lines from ionized 
gas surrounding massive young stars. For a conventional 
IMF, these lines trace stars with masses greater than 
^15 Mq, with the peak contribution from stars in the 
range 30 40 Mq. As such, the lines (and likewise the 
free-free radio continuum) represent a nearly instanta- 
neous measure of the SFR, tracing stars with lifetimes of 
~3-10 Myr (Table 1). 

The application of Ha and other emission line SFR 
tracers has expanded dramatically in the last decade, 
through very large spectroscopic surveys of samples of lo- 
cal and distant galaxies, narrow-band emission-line imag- 
ing surveys, and large imaging surveys of nearby galax- 
ies designed to study spatial variations in the SFR. Ad- 
vances in near-infrared instrumentation have also led to 
the first systematic surveys in the Paschen and Brackett 
series lines, as well as for Ha observed at high redshifts. 
Results from several of these surveys are discussed and 
referenced in §5. 

The Ha emission line remains the indicator of choice 
for observations of local and distant galaxies alike, h\it 
for moderate redshifts, the bluer visible lines have been 
applied, in particular the [O II] forbidden line doublet 
at 372.7 nm. This feature is subject to severe systematic 
uncertainties from dust attenuation and excitation varia- 
tions in galaxies. K98 published a single SFR calibration 
for the line, but subsequent analyses (e.g., Jansen et al. 
2001; Kewley et al. 2004; Moustakas et al. 2006) have 
shown that the systematic effects must be removed or at 
least calibrated to first order for reliable measurements. 
Even then, the uncertainties in the [O II]-based SFRs are 
much larger than for Ha. 
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Over the past decade, increasing attention has been 
given to measurements of the redshifted Lya hne {\rest 
= 121.6 nm) as a tracer of star-forming galaxies, espe- 
cially at high redshifts where it provides unique sensi- 
tivity to both low-mass star-forming galaxies and inter- 
galactic gas clouds or "blobs" (e.g., Ouchi et al. 2009, 

2010) . The strength of the line (8.7 times stronger than 
Ha for Case B recombination) makes it an attractive 
tracer in principle, but in realistic ISM environments the 
line is subject to strong quenching from the combination 
of resonant trapping and eventual absorption by dust, 
usually quantified in terms of a Lya escape fraction. As 
a result Lya surveys to date have been mainly used for 
identifying large samples of distant star-forming galaxies. 
Applying the line as a quantitative SFR tracer requires 
an accurate measurement of the escape fraction. Sev- 
eral ongoing studies are quantifying this parameter by 
comparing Lya fluxes of galaxies with independent SFR 
tracers such as Ha or the UV continuum (e.g., Atck ct al. 
2009; Scarlata et al. 2009; Blanc et al. 2011; Hayes et al. 

2011) . These show that the escape fraction varies wildly 
between galaxies with a range of more than two orders of 
magnitude (order 0.01 to 1), but also increases systemat- 
ically with rcdshift. It is possible that Lya will prove to 
be a powerful SFR tracer for the highest-rcdshift objects, 
but in view of the large scatter and systematic imcertain- 
ties associated with its use we have chosen not to include 
a SFR calibration of the line in this review. 

Other workers have investigated the efficacy of the in- 
frared fine-structure cooling lines, which arise in H II 
regions or PDRs, as quantitative SFR tracers. Ho & 
Keto (2007) compiled data from ISO and Spitzer on the 
[Ne Il]12.8/z.m and [Ne III] 15.6 ^m lines, and found that 
the sum of the line fluxes correlates well with hydrogen 
recombination line fluxes, with a scatter of ~0.3 dex. 
In a similar vein, Boselli ct al. (2002) and Rodriguez- 
Fernandez et al. (2006) have investigated the applica- 
bility of the [C II]158/im line as a SFR measure. They 
found good general correspondance between the [C II] lu- 
minosity and other measures of the SFR such as ionized 
gas and dust emission, but the scatter in the relation- 
ships is at least a factor of ten. Even larger variations in 
L[Cii]/-^iR were found in a more diverse sample of star- 
forming galaxies by Malhotra et al. (1997). The availabil- 
ity of a rich new set of [C II] observations from Herschel, 
combined with the detection of redshifted [C II] emission 
in submillimeter galaxies (SMGs) from ground-based in- 
struments, has sparked a resurgence of interest in this 
application. 

The largest systematic errors affecting Ha-based SFRs 
are dust attenuation and sensitivity to the population of 
the upper IMF in regions with low absolute SFRs. For 
regions with modest attenuations the ratios of Balmer 
recombination lines (Balmer decrement) can be used to 
correct for dust, and this method has been applied in 
a number of large spectrophotometric surveys of nearby 
galaxies (e.g., Kewley et al. 2002; Brinchmann et al. 2004; 
Moustakas et al. 2006). The Balmer decrements offer 
only approximate corrections for attenuation because of 
variations on scales smaller than the resolution. These 
variations may lead to an underestimate of the extinc- 
tion because lines of sight with low extinction are more 
heavily weighted within the beam. This problem can 



be addressed partly by adopting a reddening law which 
compensates in part for these departures from a pure 
foreground scheme geometry (e.g.. Chariot & Fall 2000). 
The attenuation of the emission lines is found to be sys- 
tematically higher than that of the continuum starlight 
at the same wavelengths (e.g., Calzetti et al. 1994), which 
presumably reflects the higher concentrations of dust in 
the young star-forming regions. 

As is the case with UV-based SFRs, the availability 
of far-infrared maps and luminosities for nearby galaxies 
has also made it possible to calibrate multi-wavelength 
methods for applying dust attenuation corrections to 
these measurements (§3.7). They reveal that the Balmer 
decrement provides reasonably accurate attenuation cor- 
rections in normal galaxies, where attenuations are mod- 
est (typically 0-1 mag at Ha), and care is taken to cor- 
rect the emission-line fluxes for underlying stellar absorp- 
tion. The Balmer decrement method for estimating dust 
attenuation breaks down badly, howc;ver in circumnu- 
clear starbursts or other dusty galaxies (e.g., Moustakas 
et al. 2006). With the advent of large- format integral- 
field spectrographs there is promise of applying Balmer 
decrement measurements on a spatially-resolved basis in 
galaxies (e.g., Blanc 2010; Sanchez ct al. 2012). 

The accuracy of SFRs derived from emission lines will 
also degrade in regions where the SFR is so low that 
one enters the regime of small number statistics in the 
population of massive ionizing stars. If the IMF itself 
were completely blind to the SFR, we would expect such 
effects to become apparent below Ha luminosities of or- 
der 10^*^ergss"\ or SFRs of order 0.001 Mgyr^^ (e.g., 
Cervino et al. 2003). At the very least, this effect will 
produce a much larger scatter in ionizing flux per unit 
SFR in this regime. In low-SFR regions, this sampling 
noise can be exacerbated by the short lifetimes of the 
ionizing stars, producing large temporal fluctuations in 
Ha emission even for a fixed longer-term SFR. As a re- 
sult, other tracers (e.g., FUV emission) tend to provide 
more accurate and sensitive measurements of the SFR at 
low star formation levels. 

Can these effects cause systematic errors in the SFRs? 
Pflamm-Altenburg et al. (2007, 2009) have investigated 
the effect of an IGIMF (§2.5) on SFR tracers, and shown 
that the systematic depletion of massive stars in low-SFR 
environments could cause Ha to substantially underesti- 
mate the actual SFR. Interestingly a systematic deficit 
of Ha emission in dwarf galaxies with low SFRs and in 
low SFR density regions is observed (Sullivan et al. 2000; 
BeU & Kennicutt 2001; Lee et al. 2009a; Meurer et al. 
2009). Recent work suggests however that the system- 
atic dependence of the Ha/UV ratio may be produced 
instead by temporal variations in SFRs, without having 
to resort to modifying the IMF itself (Fumagalli et al. 
2011; Weisz et al. 2012). 

3.4. Infrared Emission: The Impact of Spitzer and 
Herschel 

Interstellar dust absorbs approximately half the 
starlight in the universe and re-emits it in the infrared, so 
measurements in the IR arc essential for deriving a com- 
plete inventory of star formation. This section focusses 
on the transformational results which have come from 
the Spitzer Space Telescope (Werner et al. 2004) and the 
Herschel Space Observatory (Pilbratt et al. 2010). A re- 
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view in this journal of cxtragalactic science from Spitzcr 
can be found in Soifer et al. (2008). 

Three other space missions are beginning to influence 
this subject: the ^iiTyli?/ mission (Murakami et al. 2007), 
the Wide-field Infrared Survey Explorer (WISE) mission 
(Wright et al. 2010), and the P/ancfc mission (Planck Col- 
laboration et al. 2011a). These observatories conducted 
all-sky surveys, with AKARI imiiging in the 2.4—160 /xm 
region, WISE in the range of 3.4 — 22 /im, and Planck 
at 350 — 850 /xm (with several bands extending to longer 
wavelengths). Much of the science from these missions 
is just beginning to emerge, but they will provide very 
important results in this subject in the coming decade. 

While; early applications of dust-based SFR mea- 
surements effectively (and necessarily) assumed a one- 
component dust model, subsequent observations show 
that the dust emission is comprised of distinct compo- 
nents, each with different spatial distributions and cou- 
plings to the young stars. At wavelengths of ~5-20/im, 
the emission is dominated by molecular bands arising 
from polycyclic aromatic hydrocarbons (PAHs). Long- 
ward of A ~ 20 /xm, the emission is dominated by thermal 
continuum emission from the main dust grain population. 
Emission from small grains transiently heated by intense 
radiation fields (usually in or near star-forming regions) 
is important out to about 60 /im, whereas at longer wave- 
lengths, emission from larger grains with steady state 
temperatures dominates (Draine 2003). 

The distribution of these different emission compo- 
nents is illustrated in Figure 2, which shows Spitzer 
and Hcrschcl images of the nearby star-forming galaxy, 
NGC6946, at wavelengths ranging from 3.6-500 /xm. At 
24 /xm, the emission peaks around the youngest star- 
forming regions and H II regions, with a more diffuse 
component extending between these regions. As one pro- 
gresses to longer wavelengths, the prominence of the dif- 
fuse component increases. Recent measurements with 
Herschel show that this is mainly a physical change, and 
not an artifact of lower spatial resolution at longer wave- 
lengths (Boquien et al. 2011). This diffuse emission is 
analogous to the "infrared cirrus" emission observed in 
our own Galaxy. Interestingly PAH emission appears to 
correlate the most strongly with the longer-wavelength 
component of the thermal dust emission (e.g., Bcndo 
ct al. 2008), though it often also appears as resolved 
shells around the young star-forming regions (Helou et al. 
2004). 

These variations in the morphologies in the different 
dust emission components translate into considerable 
variations in the dust SEDs within and between galax- 
ies (e.g., Dale k Helou 2002; Dale et al. 2005; Smith 
et al. 2007), and as a consequence the conversion of in- 
frared luminosities into SFRs must change for different 
IR wavelengths. Most early applications of the dust emis- 
sion as a SFR tracer were based on the integrated total- 
infrared (TIR) emission. This parameter has the phys- 
ical advantage of effectively representing the bolometric 
luminosity of a completely dust-enshrouded stellar pop- 
ulation. The TIR-based SFR calibration derived in K98, 
applicable in the limits of complete dust obscuration and 
dust heating fully dominated by young stars, is still in 
widespread use today. For most galaxies however this 
complete wavelength coverage will not be available, so 
many workers have calibrated monochromatic SFR in- 



dices, usually tuned to one of the Spitzer or Herschel 
bands, including 24 /im (e.g., Wu ct al. 2005a; Alonso- 
Herrero et al. 2006; Calzetti et al. 2007; Relano et al. 
2007; Rieke ct al. 2009, and references therein), and 70 
and 160 jum (Calzetti et al. 2010a). The latter paper 
contains an excellent discussion comparing the various 
calibrations in the literature at the time. 

As with all of the SFR indicators, the dust emission is 
subject to important systematic effects. Just as the UV 
and visible tracers miss radiation that has been attenu- 
ated by dust, the infrared emission misses the starlight 
that is not absorbed by dust (e.g., Hirashita ct al. 2001). 
As discussed earlier, dust attenuates only about half of 
the integrated starlight of galaxies on average, so the 
infrared emission will systematically underestimate the 
SFR if the missing fraction of star formation is not in- 
corporated into the calibrations. This "missing" unat- 
tenuated component varies from essentially zero in dusty 
starburst galaxies to nearly 100% in dust-poor dwarf 
galaxies and metal-poor regions of more massive galax- 
ies. Another major systematic error works in the op- 
posite direction; in most galaxies, evolved stars (e.g., 
ages above ~100-200 Myr) contribute significantly to the 
dust heating, which tends to cause the IR luminosity to 
overestimate the SFR. The fraction of dust heating from 
young stars varies by a large factor among galaxies; in 
extreme circumnuclear starburst galaxies or individual 
star-forming regions, nearly all of the dust heating arises 
from young stars, whereas in evolved galaxies with low 
specific SFRs, the fraction can be as low as ~10% (e.g., 
Sauvage & Thuan 1992; Waltcrbos & Grcenawalt 1996; 
Cortese et al. 2008). In practical terms this means that 
the conversion factor from dust luminosity to SFR even 
in the limit of complete dust obscuration — is not fixed, 
but rather changes as a function of the stellar popula- 
tion mix in galaxies. The difference in conversion fac- 
tor between starbursts and quiescient galaxies with con- 
stant SFR, for example, is about a factor of 1.3-2 (e.g., 
Sauvage & Thuan 1992; Kennicutt et al. 2009; Hao et al. 
2011). 

The calibration of the mid-IR PAH emission as a quan- 
titative SFR tracer deserves special mention. This index 
is of particular interest for studies of galaxies at high red- 
shift, because the observed-frame 24 /im fluxes of galaxies 
at z = 1 — 3 tend to be dominated by redsliifted PAH 
emission. A number of studies have shown that the PAH 
luminosity scales relatively well with the SFR in metal- 
rich luminous star-forming galaxies (e.g., Rousscl ct al. 
2001; Pcetcrs et al. 2004; Forster Schreiber et al. 2004; 
Wu et al. 2005a; Farrah et al. 2007; Calzetti et al. 2007), 
but the PAH bands weaken dramatically below metal 
abundances of approximately 1/4 to 1/3 (e.g.. Mad- 
den 2000; Engelbracht et al. 2005; Calzetti et al. 2007; 
Smith ct al. 2007) , rendering them problematic as quan- 
titative SFR tracers in this regime. 

The best way to overcome these systematic biases is 
to combine the IR measurements with UV or visible- 
wavelength SFR tracers, to measure the unattcnuated 
starlight directly and to constrain the dust-heating stel- 
lar population (§3.7). However, in cases where only IR 
observations are available one can attempt to incorpo- 
rate corrections for these effects into the SFR calibra- 
tions themselves, and this approach has been taken by 
most authors (e.g., Calzetti et al. 2010b and references 
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therein). Fortunately, for most galaxies with moderate 
to high specific SFRs, the effects of partial dust attenu- 
ation and cirrus dust heating by evolved stars appear to 
roughly compensate for each other. For example, Kewley 
et al. (2002) showed that for a sample of spiral galaxies 
with integrated emission-line spectra, the TIR-based cal- 
ibration of K98 for dusty starburst galaxies was in good 
agreement with SFRs based on attenuation-corrected Ha 
line fluxes. However, one must avoid applying these IR- 
based recipes in environments where they are bound to 
fail, for example, in low-metallicity and other largely 
dust-free galaxies or in galaxies with low specific SFRs 
and a strong radiation field from more evolved stars. 

3.5. Radio Continuum Emission 

The centimeter-wavelength radio continuum emission 
of galaxies consists of a relatively flat-spectrum, free- 
free component, which scales with the ionizing lumi- 
nosity (subject to a weak electron temperature depen- 
dence) and a steeper spectrum synchrotron component, 
which overwhelmingly dominates the integrated radio 
emission at < 5 GHz. The free- free component can be 
separated with multi-frequency radio measurements or 
high-frequency data (e.g., Israel & van der Hulst 1983; 
Niklas et al. 1997; Murphy et al. 2011), to provide a 
photoionization-based measure of the SFR, without the 
complications of dust attenuation which are encountered 
with the Balmer lines. 

At lower frequencies the integrated emission is dom- 
inated by the synchrotron emission from charged par- 
ticles produced by supernovae. A SFR calibration has 
not been derived from first principles, but observations 
have repeatedly confirmed a tight correlation between 
this non-thermal emission and the far-infrared emission 
of galaxies, which favors its application as a SFR tracer 
(e.g., Helou et al. 1985; Condon 1992). Moreover, im- 
provements to receiver technology with the Expanded 
Very Large Array (EVLA) and other instruments have 
made the radio continuum a primary means of identify- 
ing star-forming galaxies at high redshift and estimating 
their SFRs. As a result it is appropriate to include it in 
this discussion of SFR tracers. 

Current calibrations of the radio continuum versus 
SFR relation are bootstrapped from the far-infrared cali- 
brations, using the tight radio-IR correlation. The steep 
synchrotron spectrum makes this calibration strongly 
wavelength dependent, and most are referenced to 
1.4 GHz (e.g., Yun et al. 2001; Condon 1992; Bell 2003). 
The calibration adopted here in §3.8 is derived in similar 
fashion, but adapted to the Kroupa IMFs. 

As described above, the IR-based SFR calibrations 
break down severely in faint galaxies with low dust con- 
tents, yet the radio- IR correlation remains tight and 
nearly linear over the entire luminosity range. How can 
this arise? The likely explanation can be found in Bell 
(2003) and references therein, where it is shown that the 
decrease in dust opacity in low-mass galaxies is accom- 
panied by a decline in synchrotron emission relative to 
other tracers of the SFR. This is seen most directly as a 
decline in the ratio of non-thermal radio emission (still 
dominant in the 15 GHz region) to the free- free ther- 
mal radio emission. Since the thermal radio emission is 
directly coupled to the stellar ionization rate and SFR, 
the relative decline in synchrotron luminosity must re- 



flect a physical decline per imit SFR. If correct then the 
continuity of the radio IR relation to low luminosities is 
the result of a "cosmic conspiracy" (Bell 2003) , and one 
should beware of applying the method in galaxies fainter 
than ~0.1L*. 

3.6. X-ray Emission 

Over the past decade, the integrated hard X-ray emis- 
sion of galaxies has been increasingly applied as a SFR 
tracer. The component of X-ray emission that does not 
arise from AGN accretion disks is dominated by massive 
X-ray binaries, supernovae and supernova remnants, and 
massive stars, all associated with young stellar popula- 
tions and recent star formation. Furthermore the ob- 
served 2-10 keV fluxes of galaxies are observed to be 
strongly correlated with their infrared and non-thermal 
radio continuum fluxes (e.g., Bauer et al. 2002; Ranalli 
et al. 2003; Symeonidis et al. 2011), strengthening the 
link to the SFR. 

Since there is no way to calibrate the relation between 
X-ray luminosity and SFR from first principles, the cal- 
ibration is usually bootstrapped from the infrared or ra- 
dio. Ranalli et al. (2003) derived such a calibration for 
integrated 2-10 keV X-ray luminosities, referenced to the 
K98 calibrations and IMF, and this relation is still widely 
applied today. Persic et al. (2004) derived an alternate 
calibration in terms of the hard X-ray binary luminosity 
alone, which is useful for nearby resolved galaxies. Cal- 
ibrations of the SFR in terms of X-ray luminosity and 
stellar mass have been published by Colbert et al. (2004) 
and Lehmer et al. (2010). For simplicity, we have listed 
in Table 1 the widely-applied relation of Ranalli et al. 
(2003), but adjusted to the Kroupa IMF used for the 
other calibrations. 

3.7. Composite Multi-Wavelength Tracers 

Large multi- wavelength surveys of galaxies allow tests 
and calibrations of SFR indices that combine informa- 
tion from more than one tracer and exploit the comple- 
mentary strengths of diff'erent wavelengths. Currently 
this capability is mainly limited to nearby galaxies, but 
with the expansion of far-infrared to millimeter surveys 
of high-redshift galaxies, opportunities to apply multi- 
wavelength diagnostics to distant galaxies should expand 
in the coming decade (e.g., Overzier et al. 2011; Reddy 
et al. 2012). 

The most widely explored of these methods have com- 
bined UV (usually FUV) observations with infrared mea- 
surements to construct dust-corrected SFRs, using an 
approximate energy-balancing approach. In its simplest 
form one can use a linear combination of UV and IR lu- 
minosities to correct the UV fluxes for dust attenuation: 

iuv(corr) = iuv (observed) -\-r]Lin (10) 

where the luminosities are usually calculated from flux 
densities using the deflnition L = v L^, and the coeffi- 
cient T] is dependent on the bandpasses chosen for the 
UV and IR measurements. The most common form of 
this correction uses GALEX FUV (155 nm) and total- 
infrared luminosities, hence: 

-Lfuv (corr) = Lfuv (observed) -|- r/ Xtir ( 1 1 ) 
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Other prescriptions sometimes adopt a higher order poly- 
nomial dependence on Lm (e.g., Buat et al. 2005), but 
for brevity we only discuss the linear combinations here. 
In most cases, rj < 1, because only part of the dust- 
heating radiation is contained in the FUV band, and in 
many galaxies there is significant dust heating and TIR 
emission arising from stars other than the UV-emitting 
population (IR cirrus). The coefficient rj can be cali- 
brated theoretically using evolutionary synthesis models, 
or empirically, using independent measurements of dust- 
corrected SFRs (e.g., Treyer et al. 2010; Hao et al. 2011). 

Early applications of this method were largely re- 
stricted to luminous starburst galaxies and star-forming 
regions, for which UV and IR data could be obtained 
prior to the advent of GALEX, Spitzer, and Herschel 
(e.g., Buat et al. 1999; Meurer et al. 1999; Gordon et al. 
2000), yielding values of r/ ~ 0.6. Subsequent analy- 
ses extending to normal star- forming galaxies (e.g., Hi- 
rashita et al. 2003; Kong et al. 2004; Burgarella et al. 
2005; Buat et al. 2005; Treyer et al. 2010; Hao et al. 
2011) typically produce values of ?] which are some- 
what lower (e.g., 0.46 for Hao et al. 2011), which al- 
most certainly reflects the larger contribution to dust 
heating from starlight longward of the FUV. Considering 
the wide differences in dust heating populations between 
starburst and normal galaxies, however, this difference 
of ~30% is hardly crippling, especially when it reduces 
dramatically the much larger systematic errors from UV 
or IR-based SFRs alone. 

A similar energy-balancing approach has been ap- 
plied to derive dust-corrected emission-line luminosities 
of galaxies (Calzetti et al. 2007; Zhu et al. 2008; Kenni- 
cutt et al. 2009; Treyer et al. 2010). Here the equivalent 
parameter to r] above is calibrated using independent es- 
timates of the Ha attenuation, usually the recombination 
line decrement as measured from the Paa/Ha or Ha/H/3 
ratio. Kennicutt et al. (2009) have shown that this ap- 
proach can be more broadly applied to estimate atten- 
uation corrections for other emission lines (e.g., [O II]), 
and it can be calibrated using TIR fluxes, single-band 
IR fluxes, or (;von 1.4 GHz radio continuum fluxes. An 
illustration of the effectiveness of the method is shown 
in Figure 3, taken from Kennicutt et al. (2009). The top 
panel shows the consistency between an IR-based tracer 
(24 iim. luminosity) used alone and the dust-corrected Ha 
luminosity (via the integrated Balmer decrement). There 
is a general correlation but considerable scatter and a 
pronounced non-linearity, which reflects a general corre- 
lation of average dust attenuation with the SFR itself. 
The lower panel plots the best fitting linear combination 
of Ha and IR (in this case 24 /xm) luminosities against the 
Balmer-decrement-corrected Ha luminosities. The scat- 
ter between the tracers now is reduced sevcralfold, and 
the non-linearity in the comparison with IR luminosities 
alone is essentially removed. 

These composite SFR indicators are not without sys- 
tematic uncertainties of their own. Questions remain 
about the systematic reliability of the independent at- 
tenuation calibrations (in particular with the effects of 
dust geometry on the Balmer decrements), and most of 
the recipes still retain a dependence on stellar population 
age, via the infrared term; for example the best fitting 
values of rj for H II regions and young starbursts differ 
from those of normal star forming galaxies by a factor 



of ^1.5, consistent with expectations from evolutionary 
synthesis models, as discussed in §3.8 (Kennicutt et al. 
2009). Nevertheless they represent a major improvement 
over single-wavelength tracers, and provide a valuable 
testing ground for calibrating and exploring the uncer- 
tainties in the monochromatic indicators (e.g., Calzetti 
et al. 2010b). 

Table 2 lists examples of dust attenuation corrections 
using combinations of FUV and Hq: fluxes with various 
infrared and radio tracers, taken from Hao et al. (2011) 
and Kennicutt et al. (2009), respectively. These can 
be applied in combination with the monochromatic SFR 
zeropoint calibrations listed in Table 1 (§3.8) to derive 
dust-corrected SFR measurements. 

3.8. An Updated Compendium of Integrated SFR 

Calibrations 

K98 presented calibrations for SFRs derived from UV 
continuum, TIR, Ha emission-line, and [O II] emission- 
line luminosities, which have come into common usage in 
the field. The considerable expansion of the subject to 
other wavelengths and SFR diagnostics since that time 
motivates a revisit of these calibrations. 

Nearly all of these calibrations are based on evolu- 
tionary synthesis models, in which the emergent SEDs 
are derived for synthetic stellar populations with a pre- 
scribed age mix, chemical composition, and IMF. The 
K98 calibrations employed a mix of models from the lit- 
erature, and assumed a single power-law IMF (Salpeter 
1955), with mass limits of 0.1 and 100 M©. This IMF 
gave satisfactory SFR calibrations relative to that of 
more realistic IMFs for Ha, but for other wavelengths, 
the relative calibrations using different tracers are sensi- 
tive to the precise form of the IMF. Today most work- 
ers calibrate SFR tracers using a modern IMF with a 
turnover below ^1 Mq, for example the IMF of Kroupa 
& Weidner (2003), with a Salpeter slope (a* = —2.35) 
from 1-100 M© and a* = -1.3 for 0.1-1 Mq. The cal- 
ibrations presented here use this IMF, but the IMF fit 
from Chabricr (2003) yields nearly identical results (e.g., 
Chomiuk & Povich 2011). The past decade has also seen 
major improvements in the stellar evolution and atmo- 
spheric models which are used to generate the synthetic 
SEDs. The results cited here use the Starburst99 models 
of Leitherer et al. (1999), which are regularly updated in 
the on-line version of the package. 

Table 1 presents in compact form calibrations for a 
suite of SFR tracers in the form: 

log (Mq yr-1) = log - log C, (12) 

The table lists for each tracer the units of luminosity 

(i,j.), the logarithmic SFR calibration constant Cx, and 
the primary rcfcrcncc(s) for the calibration. For meth- 
ods presented in K98, we also list the scaling constant 
between the new (Kroupa IMF, SB99 model) SFRs and 
those from K98. As in K98, the recipes for infrared dust 
emission are in the limit of complete dust attenuation 
and continuous star formation over a period of 0-100 
Myr, which is appropriate for a typical dust-obscured 
starburst galaxy. Virtually all of the SFR calibrations 
presented in Table 1 are taken from the literature, and 
we strongly encourage users of these calibrations to cite 
the primary sources, whether or not they choose to cite 
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this paper as well. 

The second column of Table 1 lists the approximate 
age sensitivity of the different star formation tracers. 
These were estimated using the Starburst99 models in 
the approximation of constant star formation. The sec- 
ond number lists the mean stellar age producing the rele- 
vant emission, while the third column lists the age below 
which 90% of the relevant emission is produced. For 
the dust emission in the infrared these ages can only be 
estimated, because they depend on the detailed star for- 
mation history over periods of up to 100 Myr and longer, 
and they are also convolved with the level of dust attenu- 
ation as a hmction of stellar age. The numbers are given 
for the assumed starburst timescales above; for normal 
galaxies the 90th percentile age can be 500 Myr or longer. 

All of the results were calculated for solar metal abun- 
dances, and readers should beware that all of the cal- 
ibrations are sensitive to metallicity. These have been 
estimated by several authors using evohitionary synthe- 
sis models (e.g.. Smith et al. 2002; Raiter et al. 2010), but 
the precise dependences are sensitive to the details of the 
input stellar models, in particular how effects of stellar 
rotation are modelled. The models cited above (which do 
not adopt a dependence of rotation on metallicity) show 
that a decrease in metal abundance by a factor of ten in- 
creases the FUV luminosity of a fixed mass and IMF pop- 
ulation by ~ 0.07±0.03 dex (for the IMFs assumed here), 
while the ionizing luminosity is more sensitive, increasing 
by ~ 0.4 ±0.1 dex for a tenfold decrease in Z/Zq. The 
change in infrared luminosities for completely obscured 
regions should roughly track that for the FUV luminos- 
ity, but in most low-metallicity environments, the dust 
opacity will be severely reduced, producing a sharp fall 
in infrared emission for a given SFR. 

3.9. The Challenge of Spatially-Resolved Star 
Formation Rates in Galaxies 

Nearly all of the diagnostic methods described up to 
now have been designed for measuring integrated SFRs 
of galaxies or for regions such as starbursts, containing 

thousands (or more) of massive stars. The resulting SFR 
prescriptions implicitly assume that local variations in 
stellar age mix, IMF population, and gas/dust geometry 
largely average out when the integrated emission of a 
galaxy is measured. 

With the advent of high-resolution maps of galaxies in 
the UV, IR, emission lines, and radio (and integral-field 
spectroscopic maps in the visible and near- infrared), one 
would like of course to extend these diagnostic methods 
to create fully-sampled spatially-resolved "SFR maps" 
of galaxies. This extrapolation to much smaller regions 
within galaxies (or to galaxies with extremely low SFRs), 
however, is not at all straightforward. At smaller linear 
scales, nearly all of the statistical approximations cited 
above begin to break down. 

First, when the SFR in the region studied drops be- 
low -0.001-0.01 Mq yr-i (depending on the SFR tracer 
used), incomplete sampling of the stellar IMF will lead 
to large fluctuations in the tracer luminosity for a fixed 
SFR. These begin to become problematic for luminosi- 
ties of order 10"^* — 10^^ ergss^^ for emission-line, UV, or 
IR tracers, and they are especially severe for the ionized 
gas tracers, which are most sensitive to the uppermost 
parts of the stellar IMF. For actively star forming normal 



disk galaxies, this onset of stochasticity typically occurs 
on spatial scales of order 0.1 1 kpc, but the region can 
be considerably larger in galaxies, or parts of galaxies, 
with lower SFRs. Note that this breakdown on local 
scales occurs regardless of whether the IMF itself varies 
systematically. 

As a prime example of this stochasticity, consider what 
astronomers in M51 would observe if they examined the 
solar neighborhood in the Milky Way from a distance 
of 10 Mpc. In a pixel of 100 pc radius centered on the 
Sun they would observe no molecular gas and no star 
formation. If they degraded their resolution to a radius 
of 300 pc, they would pick up all the Gould Belt clouds, 
but no localized Ha emission. Within a slightly larger 
radius of 500 pc, their beam would include Orion, with 
its O stars and H II region. (For this example, we con- 
sider only the stars and emission from the Orion Nebula 
Cluster (ONC) stars and emission, for which we have 
good numbers, not the full Orion clouds and OB associ- 
ations.) This roughly doubles the number of YSOs. Our 
M51 observers would measure SmoZ — 0.11 Mq pc^^, 
but if they applied the relation between SFR and Hq lu- 
minosity from K98, they would derive S(SFR) a factor 
of 10 lower than the actual SFR in Orion, based on the 
actual stellar content of the ONC and a relatively long 
timescalc of 3 Myr (Chomiuk & Povich 2011), and they 
would miss the total SFR within 500 pc by a factor of 20 
because the other clouds produce no O stars. This severe 
underestimate for Orion itself from Ha emission results 
because the Orion chistcr is too small to fully populate 
the IMF; its earliest spectral type is 07 V, and so the star 
cluster produces relatively little ionizing luminosity rel- 
ative to its total mass and SFR. If the distant observers 
used the total infrared emission of Orion instead, they 
would still underestimate its SFR by a factor of 8 (Lada 
et al. 2012) (a factor of 10 with the newer conversions in 
Table 1). 

Secondly, when the spatial resolution of the SFR 
measurements encompass single young clusters, such 
as Orion, the assumption of continuous star formation 
which is embedded in the global SFR recipes breaks down 
severely. In such regions that emission at all wavelengths 
will be dominated by a very young population with ages 
of typically a few Myr, which will be shorter than the av- 
eraging times assumed for all but the emission-line trac- 
ers. These changes affect both the relative limiinosities of 
star formation tracers in different bands and of course the 
interpretation of the "star formation rate" itself. Tech- 
nically speaking the luminosities of young star-forming 
clusters only provide information on the masses of the 
regions studied, and converting these to SFRs requires 
independent information on the ages and/or age spreads 
of the stars in the region. 

A third measuring bias will set in if the resolution of 
the measurements becomes smaller than the Stromgren 
diameters of H II regions or the corresponding dust emis- 
sion nebulae, which tend to match or exceed those of 
the HII regions (e.g., Watson et al. 2008; Prescott et al. 
2007). This scale varies from <100pc in the Milky Way 
to 200 500 pc in actively star- forming galaxies such as 
NGC 6946 (Figure 2). On these small scales indirect trac- 
ers of the SFR (e.g.. Ha, IR) tend to trace the surfaces 
and bubbles of clouds rather than the young stars. By 
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the same token, much of the Ha and dust emission in 
galaxies, typicaUy 30 -60%, is emitted by diffuse ionized 
gas and dust located hundreds of parsecs or more from 
any young stars (e.g., Oey et al. 2007; Dale et al. 2007). 
Such emission produces a false positive signal of star for- 
mation, and care is needed to account for its effects when 
mapping the SFR within galaxies. 

Taken together, these factors complicate, but do not 
prevent, the construction of 2D maps of SFRs, so long as 
the methodology is adapted to the astrophysical applica- 
tion. For example, radial profiles of SFR distributions in 
galaxies may still be reliably derived using the integrated 
SFR calibrations, providing that large enough annuli arc 
used to assure that the IMF is fully populated in ag- 
gregate, and the SFR is regarded to be averaged over 
time scales of order 100 Myr. Likewise a 2D distribu- 
tion of star formation over the last 5 Myr or so can be 
derived from a short-lived tracer such as Ha emission, if 
the spatial covcragt; is limited to very young regions, and 
structure on scales smaller than individual H II regions is 
ignored. Maps of integrated Ha and IR emission can be 
used to study the population of star-forming regions and 
the star formation law, when restricted to young regions 
where the physical association of gas and stars is secure 
(§6.4). An alternative approach is to apply direct stellar 
photosphcric tracers, such as the ultraviolet continuum 
emission on large scales (bearing in mind the variable 
10-300 Myr time scales traced by the UV), or resolved 
stellar tracers such as YSOs or deep visible-wavelength 
color-magnitude diagrams, to map the distributions of 
young stars directly (§3.1). Likewise it should be pos- 
sible to use pixel-resolved SEDs of galaxies in the UV- 
visible to derive dust-corrected UV maps, and possibly 
apply local corrections for age in the maps. With such a 
multiplicity of approaches, we anticipate major progress 
in this area, which will be invaluable for imdcrstanding 
in detail the patterns of star formation in galaxies and 
connecting to detailed studies within the MW. 

4. THE LOCAL PERSPECTIVE: FROM THE 
INSIDE LOOKING OUT 

4.1. Outline of low-mass star formation 

The formation of low-mass stars can be studied in 
greatest detail because it occurs in relatively nearby 
clouds and sometimes in isolation from other forming 
stars. The established paradigm provides a point of com- 
parison for more massive, distant, and clustered star for- 
mation. 

Physically, an individual star forming event, which 
may produce a single star or a small number multiple 
system, proceeds from a prestellar core, which is a grav- 
itationally bound starless core (sec §2.2, di Francesco 
et al. 2007 for definitions), a dense region, usually within 
a larger molecular cloud. Prestellar cores arc centrally 
condensed, and can be modeled as Bonnor-Ebcrt spheres 
(Ward- Thompson et al. 1999; Evans et al. 2001; Kirk 
et al. 2005), which have nearly power-law (n(r) oc r~^) 
envelopes around a core with nearly constant density. 

Collapse leads to the formation of a first hydrostatic 
core in a small region where the dust continuum emis- 
sion becomes optically thick; this short-lived core is still 
molecular and grows in mass by continued infall from the 
outer layers (Larson 1969; Boss & Yorke 1995; Omukai 



2007; Stabler & Palla 2005). When the temperature 
reaches 2000 K, the molecules in the first hydrostatic 
core collisionally dissociate, leading to a further collapse 
to the true protostar, still surrounded by the bulk of the 
core, often referred to as the envelope. Rotation leads 
to flattening and a centrifugally supported disk (Terebey 
et al. 1984). Magnetic fields plus rotation lead to winds 
and jets (Shang et al. 2007; Pudritz et al. 2007), which 
drive molecular outflows from young stars of essentially 
all masses (Arce et al. 2007; Wu et al. 2004). 

Stage sources have protostars, disks, jets, and en- 
velopes with more mass in the envelope than in the star 
plus disk (Andre et al. 1993). Stage I sources are similar 
to Stage 0, but with less mass in the envelope than in star 
plus disk. Stage II sources lack an envelope, but have a 
disk, while Stage HI sources have little or no disk left. 
There are several alternative quantitative definitions of 
the qualitative terms used here (Robitaille et al. 2006; 
Crapsi et al. 2008). 

Stages I through HI are usually associated with SED 
Classes I through III defined by SED slopes rising, falling 
slowly, and falling more rapidly from 2 to 25 /um (Lada 
& Wilking 1984). A class intermediate between I and II, 
with a slope near zero (Flat SED) was added by Greene 
et al. (1994). In fact, orientation effects in the ear- 
lier stages can confuse the connection between class and 
stage; without detailed studies of each source, one cannot 
observationally determine the stage without ambiguity. 
For this reason, most of the subsequent discussion will 
use classes despite their questionable correspondence to 
physical configurations. Evans et al. (2009) give a his- 
torical account of the development of the class and stage 
nomenclature. 

For our present purposes, the main point is to estab- 
lish time scales for the changing observational signatures, 
so the classes arc useful. Using the boundaries between 
classes from Greene et al. (1994), Evans et al. (2009) 
found that the combined Class O/I phases last ~ 0.5 Myr, 
with a similar duration for Flat SEDs, assuming a con- 
tinuous flow through the classes for at least the Class II 
duration of 2 ± 1 Myr (e.g., Mamajek 2009), which is 
essentially the age when about half the stars in a dated 
cluster lack infrared excesses; hence all durations can be 
thought of as half-lives. 

To summarize, we would expect a single star-forming 
core to be dominated by emission at far-infrared and sub- 
millimeter wavelengths for about 0.5 Myr, near-infrared 
and mid-infrared radiation for about 1.5 Myr and near- 
infrared to visible light thereafter. However, the duration 
of emission at longer wavelengths can be substantially 
lengthened by material in surrounding clumps or clouds 
not directly associated with the forming star. 

The star formation rates and efficiencies can be calcu- 
lated for a set of 20 nearby clouds {{d) = 275 pc) with 
uniform data from Spitzer. The star formation rate was 
calculated from equation 9 and the assumptions in §3.1 
by Heiderman et al. (2010), who found a wide range of 
values for M*, with a mean value for the 20 clouds of 
M^, = 39 ± 18 Mq Myr~^. The star formation efficiency 
(e = Mi, /{Mi, + Mgas)) can only be calculated over the 
last 2 Myr because surveys are incomplete at larger ages; 
averaging over all 20 clouds, only 2.6% of the cloud mass 
has turned into YSOs in that period and (S(SFR)) = 1.2 
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M© yr~-'^kpc^^ (Hciderman et al. 2010). For individual 
clouds, e ranges from 2% to 8% (Evans ct al. 2009; Pe- 
terson et al. 2011). The final efficiency will depend on 
how long the cloud continues to form stars before being 
disrupted; cloud lifetimes are poorly constrained (McKee 
& Ostriker 2007; Bergin & Tafalla 2007), and they may 
differ for clouds where massive stars form. 

For comparison to extragalactic usage, the mean tdep = 
1/e' (§1.2) is about 82 Myr for seven local clouds, longer 
than most estimates of cloud lifetime, and much longer 
than either (tff) (1.4 Myr) or (tcross) (5.5 Myr) (Evans 
et al. 2009). The (tdep) for local clouds is, on the other 
hand, about 10% of the tdep for the MW molecular gas 
as a whole (§5.1). 

4.2. Formation of Clusters and High-mass Stars 

The study of young clusters provides some distinct op- 
portunities. By averaging over many stars, a character- 
istic age can be assigned, with perhaps more reliability 
than can be achieved for a single star. As noted (§4.1), 
the set point for ages of all SED classes is determined 
by the fraction of mid-infrared excesses in clusters. This 
approach docs implicitly assume that the cluster forms 
"coevally" , by which one really means that the spread in 
times of formation is small compared to the age of the 
cluster. Because objects ranging from prestellar cores to 
Class III objects often coexist (Rebull et al. 2007), this 
assumption is of dubious reliability for clusters with ages 
less than about 5 Myr, exactly the ones used to set the 
timescales for earlier classes. Note also that assuming 
coeval formation in this sense directly contradicts the 
assumption of continuous flow through the classes. We 
live with these contradictions. 

Most nearby clusters do not sample very far up the 
IMF. The nearest young cluster that has formed O stars 
is 415-430 pc away in Orion (Menten et al. 2007; Hirota 
et al. 2007), and most lie at much larger distances, mak- 
ing study of the entire IMF difficult. 

Theoretically, the birthplace of a cluster is a clump 
(Wilhams et al. 2000). The clumps identified by "CO 
maps in nearby clouds are dubious candidates for the 
reasons given in §2.5. The structures identified by Kain- 
ulainen et al. (2009) in the power-law tail of the probabil- 
ity distribution function are better candidates, but most 
are still unbound. Objects found in some surveys of dust 
continuum emission (§2.3) appear to be still better can- 
didates. For example, analysis (Schlingman et al. 2011; 
Dunham et al. 2011b) of the sources found in millimeter- 
wave continuum emission surveys of the MW (§2.3) indi- 
cates typical volume densities of a few x 10^ cm^'^, com- 
pared to about 200 cm~^ for ^'^CO structures. Mean 
surface densities are about 180 M© pc~^, higher than 
those of gas in the the power-law tail (Sgas ~ 40 — 80 
Mq pc^^). These structures have a wide range of sizes, 
with a median of 0.75 to 1 pc. 

Studies of regions with signposts of massive star for- 
mation, using tracers requiring higher densities (e.g., 
Plume et al. 1992; Beuther et al. 2002) or millimeter 
continuum emission from dust (e.g., Beuther et al. 2002; 
Mueller et al. 2002) have identified slightly smaller struc- 
tures (r ~ 0.5 pc) that are much denser, indeed denser 
on average than cores in nearby clouds. In addition to 
having a mass distribution consistent with that of clus- 



ters (§2.5), many have far-infrared luminosities consis- 
tent with the formation of clusters of stars with masses 
ranging up to those of O stars. The properties of the 
Plume et al. (1992) sample have been studied in a series 
of papers (Plume et al. 1997; Mueller et al. 2002; Shirley 
et al. 2003), with the most recent summary in Wu et al. 
(2010a). The properties of the clumps depend on the 
tracer used. Since the HCN J = 1 — > line is used in 
many extragalactic studies, we will give the properties 
as measured in that line. The mean and median FWHM 
are 1.13 and 0.71 pc; the mean and median infrared lumi- 
nosities are 4.7 x 10^ and 1.06 x 10^ L©; the mean and 
median virial mass are 5300 and 2700 M©; the mean and 
median siirface densities are 0.29 and 0.28 gm cni^^; and 
the mean and median of the average volume densities (n) 
are 3.2 x 10^ and 1.6 x 10* cm~^. As lines tracing higher 
densities are used, the sizes and masses decrease, while 
the surface densities and volume densities increase, as 
expected for centrally condensed regions. Very similar 
results were obtained from studies of millimeter contin- 
uum emission from dust (Faiindez et al. 2004) toward a 
sample of southern hemisphere sources surveyed in CS 
J = 2 ^ 1 (Bronfman et al. 1996). 

Embedded clusters provide important testing grounds 
for theories of star formation. Using criteria of 35 mem- 
bers with a stellar density of 1 M©pc^'^ for a cluster, 
Lada & Lada (2003) argue that most stars form in clus- 
ters, and 90% of those are in rich clusters with more than 
100 stars, but that almost all clusters (> 93%) dissipate 
as the gas is removed, a process they call "infant mor- 
tality." With more complete surveys enabled by Spitzer, 
the distributions of numbers of stars in clusters and stel- 
lar densities are being clarified. Allen et al. (2007) found 
that about 60% of young stars within 1 kpc of the Sun are 
in clusters with more than 100 members, but this num- 
ber is heavily dominated by the Orion Nebula cluster. 
Drawing on the samples from Spitzer surveys of nearly 
all clouds within 0.5 kpc, Bressert et al. (2010) found a 
continuous distribution of surface densities and no ev- 
idence for a bimodal distribution, with distinct "clus- 
tered" and "distributed" modes. They found that the 
fraction of stars that form in clusters ranged from 0.4 
to 0.9, depending on which definition of "clustered" was 
used. Even regions of low-mass star formation, often 
described as distributed star formation, are quite clus- 
tered and the youngest objects (Class I and Flat SED 
sources, see §4.1) are very strongly concentrated to re- 
gions of high extinction, especially after the samples have 
been culled of interlopers (Fig. 4, §6.4). Similarly, stud- 
ies of three clouds found that 75% of prestellar cores lay 
above thresholds in Ay of 8, 15, and 20 mag, while most 
of the cloud mass was at much lower extinction levels 
(Enoch ct al. 2007). 

A Spitzer study (Gutermuth et al. 2009) of 2548 YSOs 
in 39 nearby (d < 1.7 kpc), previously known (primar- 
ily from the compilation by Porras et al. 2003) young 
clusters, but excluding Orion and NGC2264, found the 
following median properties: 26 members, core radius of 
0.39 pc, stellar surface density of 60 pc~^, and embed- 
ded in a clump with Ak — 0.8 mag, which corresponds 
to Ay = 7.1 mag. Translating to mass surface density 
using (M*) =0.5 M©, the median stellar mass surface 
density would be 30 M© pc~^ and the gas surface den- 
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sity would be 107 M0 pc^^. The distributions arc often 
elongated, with a median aspect ratio of 1.82. The me- 
dian spacing between YSOs, averaged over all 39 clusters, 
is 0.072 ± 0.006 pc, comparable to the size of individual 
cores, and a plausible scale for Jeans fragmentation. The 
distributions are all skewed toward low values, with a tail 
to higher values. 

While various definitions of "clustered" have been 
used, one physically meaningful measure is a surface den- 
sity of ~ 200 YSOs pc-2 (Gutermuth et al. 2005), be- 
low which individual cores are likely to evolve in rela- 
tive isolation (i.e., the timescale for infall is less than the 
timcscalc for core collisions). With their sample, Bressert 
ct al. (2010) found that only 26% were likely to interact 
faster than they collapse. However, that statistic did not 
include the Orion cluster, which exceeds that criterion. 

The fact that the Orion cluster dominates the local 
star formation warns us that our local sample may be 
unrepresentative of the Galaxy as a whole. Leaving aside 
globular clusters, there are young clusters that are much 
more massive than Orion. Portegies Zwart et al. (2010) 
have cataloged massive (M+ > 10"* M0), young (age < 
100 Myr) clusters (12) and associations (13), but none of 
these are more distant than the Galactic Center, so they 
are clearly undercounted. The mean half-light radius of 
the 12 clusters is 1.7 ± 1.3 pc compared to 11.2 ± 6.4 
pc for the associations. For the clusters, (logM*(M0)) = 
4.2±0.3. Three of the massive dense clumps from the Wu 
ct al. (2010a) study have masses above 10^ Mq and arc 
plausible precursors of this class of clusters. Still more 
massive clusters can be found in other galaxies, and a 
possible precursor (Mdoud > 1 x 10^ Mq within a 2.8 pc 
radius) has recently been identified near the center of the 
MW (Longmore et al. 2012). 

The topic of clusters is connected to the topic of mas- 
sive stars because 70% of O stars reside in young clus- 
ters or associations (Gies 1987). Furthermore, most of 
the field population can be identified as runaways (de 
Wit et al. 2005), with no more than 4% with no evi- 
dence of having formed in a cluster. While there may be 
exceptions (see Zinnecker & Yorke 2007 for discussion), 
the vast majority of massive stars form in clusters. The 
most massive star with a dynamical mass (NGC 3603- 
Al) weighs in at 116 ± 31 Mq and is a 3.77 day binary 
with a companion at 89 ± 16 Mq (Schnurr et al. 2008). 
Still higher initial masses (105-170 Mq) for the stars in 
NGC3603 and even higher in R136 (165-320 Mq) have 
been suggested (Crowther ct al. 2010). Many of the most 
massive stars exist in tight (orbital periods of a few days) 
binaries (Zinnecker & Yorke 2007). For a recent update 
on massive binary properties, see Sana & Evans (2011). 

There is some evidence, summarized by Zinnecker & 
Yorke (2007), that massive stars form only in the most 

massive molecular clouds, with max(M^) oc Mdoud*''^^ 
suggested by Larson (1982). Roughly speaking, it takes 
Mcioud = 10^ Mq to make a 50 Mq star. Recognizing 
that clumps are the birthplaces of clusters and that effi- 
ciencies are not unity could make the formation of mas- 
sive stars even less likely. The question (§2.5) is whether 
the absence of massive stars in clumps or clusters of 
modest total mass is purely a sampling effect (Fumagalli 
et al. 2011; Calzetti et al. 2010a and references therein) 
or a causal relation (Weidner & Kroupa 2006). If causal, 



differences in the upper mass limit to clouds or clumps 
(§2.5) in a galaxy could limit the formation of the most 
massive stars. Larson (1982) concluded that his correla- 
tion could be due to sampling. As a concrete example, 
is the formation of a 50 Mq star as likely in an ensemble 
of 100 clouds, each with 10^ Mq, as it is in a single 10^ 
Mq cloud? Unbiased surveys of the MW for clumps and 
massive stars (§5.1) could allow a fresh look at this ques- 
tion, with due regard for the difficulty of distinguishing 
"very rarely" from "never." 

4.3. Theoretical Aspects 

The fundamental problem presented to modern theo- 
rists of star formation has been to explain the low effi- 
ciency of star formation on the scale of molecular clouds. 
Early studies of molecular clouds concluded that they 
were gravitationally bound and should be collapsing at 
free fall (e.g., Goldreich & Kwan 1974). Zuckerman & 
Palmer (1974) pointed out that such a picture would pro- 
duce stars at 30 times the accepted average recent rate 
of star formation in the Milky Way (§5.1) if stars formed 
with high efficiency. Furthermore, Zuckerman & Evans 
(1974) found no observational evidence for large scale 
collapse and suggested that turbulence, perhaps aided by 
magnetic fields, prevented overall collapse. This sugges- 
tion led eventually to a picture of magnetically subcriti- 
cal clouds that formed stars only via a redistribution of 
magnetic flux, commonly referred to as ambipolar diffu- 
sion (Shu et al. 1987; Mouschovias 1991), which resulted 
in cloud lifetimes about 10 times the free-fall time. If, 
in addition, only 10% of the cloud became supercritical, 
a factor of 100 decrease in star formation rate could be 
achieved. 

Studies of the Zeeman effect in OH have now provided 
enough measurements of the line-of-sight strength of the 
magnetic field to test that picture. While there are still 
controversies, the data indicate that most clouds (or, 
more precisely, the parts of clouds with Zeeman mea- 
surements) are supercritical or close to critical, but not 
strongly subcritical (Crutcher et al. 2010). Pictures of 
static clouds supported by magnetic fields are currently 
out of fashion (McKec & Ostrikcr 2007), but magnetic 
fields are almost certain to play a role in some way (for 
a current review, sec Crutcher 2012, this volume). Sim- 
ulations of turbulence indicate a fairly rapid decay, even 
when magnetic fields are included (Stone et al. 1998). 
As a result, there is growing support for a more dynam- 
ical picture in which clouds evolve on a crossing time 
(Elmegreen 2000). However, this picture must stiU deal 
with the Zuckerman-Palmer problem. 

There are two main approaches to solving this prob- 
lem at the level of clouds, and they are essentially exten- 
sions of the two original ideas, magnetic fields and tiirbu- 
lence, into larger scales. One approach, exemplified by 
Vazquez-Semadeni et al. (2011) argues that clouds are 
formed in colliding flows of the warm, neutral medium. 
They simulate the outcome of these flows with magnetic 
fields, but without feedback from star formation. Much 
of the mass is magnetically subcritical and star forma- 
tion happens only in the supercritical parts of the cloud. 
A continued flow of material balances the mass lost to 
star formation so that the star formation efficiency ap- 
proaches a steady value, in rough agreement with the 
observations. 
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A second picture, exemplified by Dobbs et al. (2011), 
is that most clouds and most parts of clouds are not 
gravitationally bound but arc transient objects. Their 
simulations include feedback from star formation, but 
not magnetic fields. During cloud collisions, material is 
redistributed, clouds may be shredded, and feedback re- 
moves gas. Except for a few very massive clouds, most 
clouds lose their identity on the timescales of a few Myr. 
In this picture, the low efficiency simply reflects the frac- 
tion of molecular gas that is in bound structures. 

It is not straightforward to determine observationally if 
clouds are bound, especially when they have complex and 
filamentary boundaries. Heyer et al. (2001) argued that 
most clouds in the outer galaxy with M > Mq were 
bound, but clumps and clouds with M < 10^ Mq were 
often not bound. Heyer et al. (2009) address the issues 
associated with determining accurate masses for molecu- 
lar clouds. In a study on the inner Galaxy, Roman-Duval 
et al. (2010) concluded that 70% of molecular clouds 
(both in mass and number) were bound. 

If molecular clouds arc bound and last longer than a 
few crossing times, feedback or turbulence resulting from 
feedback is invoked. For low-mass stars, outflows provide 
the primary feedback (e.g., Li & Nakamura 2006) while 
high-mass stars add radiation pressure and expanding 
HIT regions (e.g., Murray 2011). In addition, the effi- 
ciency is not much higher in most clumps, so the problem 
persists to scales smaller than that of clouds. 

Another longstanding problem of star formation the- 
ory has been to explain the IMF (§2.5). This topic has 
been covered by many reviews, so we emphasize only two 
aspects. The basic issue is that the typical conditions in 
star forming regions suggest a characteristic mass, either 
the Jeans mass or the Bonnor-Ebert mass, around 1 Mq 
(e.g., Bonnell et al. 2007). Recently, Krumholz (2011) 
has derived a very general expression for a character- 
istic mass of 0.15 Mq, with only a very weak depen- 
dence on pressure. However, we see stars down to the 
hydrogen-burning limit and a continuous distribution of 
brown dwarfs below that, extending down even to masses 
lower than those seen in extrasolar planets (AUers et al. 
2006). An extraordinarily high density would be required 
to make such a low mass region unstable. 

On the other end, making a star with mass over 100 
times the characteristic mass is challenging. Beuther 
et al. (2007a) provide a nice review of both observa- 
tions and theory of massive star formation. While dense 
clumps with mass much greater than 100 Mq are seen, 
they are likely to fragment into smaller cores. Fragmen- 
tation can solve the problem of forming low mass objects, 
but it makes it hard to form massive objects. Simula- 
tions of unstable large clumps in fact tend to fragment so 
strongly as the mean density increases that they overpro- 
duce brown dwarfs, but no massive stars (e.g., Klessen 
et al. 1998; Martel et al. 2006). This effect is caused 
by the assumed isothermality of the gas, because the 
Jeans mass is proportional to {Tk^ /n)^-^. Simulations 
including radiative feedback, acting on a global scale, 
have shown that fragmentation can be suppressed and 
massive stars formed (e.g.. Bate 2009; Urban et al. 2010; 
Krumholz et al. 2010). Krumholz k McKee (2008) have 
argued that a threshold clump surface density of 1 gm 
cm~^ is needed to suppress fragmentation, allowing the 



formation of massive stars. 

The radiative feedback, acting locally, can in principle 
also limit the mass of the massive stars through radiation 
pressure. This can be a serious issue for the formation 
of massive stars in spherical geometries, but more real- 
istic, aspherical simulations show that the radiation is 
channeled out along the rotation axis, allowing contin- 
ued accretion through a disk (Yorke & Sonnhalter 2002; 
Krumholz et al. 2009a; Kuiper et al. 2011). 

The picture discussed so far is basically a scaled-up 
version of the formation of low mass stars by accretion 
of material from a single core (sometimes called Core 
Accretion). An alternative picture, called Competitive 
Accretion, developed by Klessen et al. (1998), Bate et al. 
(2003), and Bonnell et al. (2003), builds massive stars 
from the initial low-mass fragments. The pros and cons 
of these two models are discussed in a joint paper by the 
leading protagonists (Krumholz & Bonnell 2009). They 
agree that the primary distinction between the two is 
the original location of the matter that winds up in the 
star: in Core Accretion models, the star gains the bulk 
of its mass from the local dense core, continuing the con- 
nection of the core mass function to the initial stellar 
mass function to massive stars; in Competitive Accre- 
tion models, the more massive stars collect most of their 
mass from the larger clump by out-competing other, ini- 
tially low mass, fragments. A hybrid picture, in which 
a star forms initially from its parent core, but then con- 
tinues to accrete from the surrounding clump, has been 
advanced (Myers 2009, 2011), with some observational 
support (e.g., Longmore et al. 2011). 

Observational tests of these ideas are difficult. The 
Core Accretion model requires that clumps contain a 
core mass function extending to massive cores. Some 
observations of massive dense clumps have found sub- 
structure on the scale of cores (Beuther et al. 2007b; 
Brogan et al. 2009), but the most massive cores identified 
in these works are 50-75 Mq. Discussions in those refer- 
ences illustrate the difficulties in doing this with current 
capabilities. ALMA will make this kind of study much 
more viable, but interpretation will always be tricky. The 
Competitive Accretion model relies on a continuous flow 
of material to the densest parts of the clump to feed the 
growing oligarchs. Evidence for overall inward flow in 
clumps is difficult; evidence for it is found in some sur- 
veys (Wu & Evans 2003; Fuller et al. 2005; Reiter et al. 
2011a), but not in others (Purcell et al. 2006). Since spe- 
cial conditions are required to produce an inflow signa- 
ture, the detection rates may underestimate the fraction 
with inflow. The dense clumps are generally found to be 
centrally condensed (Beuther et al. 2002; Mueller et al. 
2002; van der Tak et al. 2000), which can also suppress 
over- fragmentation. 

Other possible tests include the coherence of outflows 
in clusters, as these should be distorted by sufficiently 
rapid motions. The kinematics of stars in forming clus- 
ters, which should show more velocity dispersion in the 
competitive accretion models, provide another test. As- 
trometric studies are beginning to be able to constrain 
these motions (Rochau et al. 2010), along with cleaner 
separation of cluster members from field objects. Rel- 
ative motions of cores within clumps appear to be very 
low (< 0.1 /kms), challenging the Competitive Accretion 
model (Walsh et al. 2004). 
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Theoretical studies include evolutionary tracks of prc- 
main-sequcncc stars (e.g., Chabrier & Baraffe 2000 for 
low mass stars and substellar objects). These evolu- 
tionary calculations are critical for determining the ages 
of stars and clusters. For more massive stars, it is es- 
sential to include accretion in the evolutionary calcula- 
tions (Palla & Stabler 1992) and high accretion rates 
strongly affect the star's evolution (Zinnecker & Yorke 
2007; Hosokawa & Omukai 2009). Assumptions about 
accretion and initial conditions may also have substantial 
consequences for the usual methods of determining ages 
of young stars (Baraffe et al. 2009; Baraffe & Chabrier 
20i0; Hosokawa ct al. 2011). 

4.4. Summary Points from the Local Perspective 

The main lessons to retain from local studies of star 
formation as we move to the scale of galaxies are sum- 
marized here. 

1. Star formation is not distributed smoothly over 
molecular clouds but is instead highly concentrated 
into regions of high extinction or mass surface den- 
sity, plausibly associated with the theoretical idea 
of a cluster-forming clump. This is particularly ap- 
parent when prestellar cores or the youngest pro- 
tostars are considered (e.g.. Figure 4). In contrast, 
most of the mass in nearby molecular clouds is in 
regions of lower extinction. 

2. By counting YSOs in nearby clouds, one can ob- 
tain reasonably accurate measures of star forma- 
tion rate and efficiency without the uncertainties of 
extrapolation from the high-mass tail of the IMF. 
The main source of uncertainty is the ages of YSOs. 

3. There is no obvious bimodality between "dis- 
tributed" and "clustered" star formation, but dense 

clusters with massive stars tend to form in regions 
of higher mean density and tiirbulcncc, which arc 
centrally condensed. Plausible precursors of quite 
massive (up to 10** or perhaps 10^ Mq) clusters can 
be found in the MW. 

4. Despite much theoretical progress, the challenge of 
explaining the low efficiency of star formation, even 
in regions forming only low-mass stars, remains. 
Similarly, an understanding of the full IMF, from 
brown dwarfs to the most massive stars, remains 
elusive. 

5. THE GALACTIC PERSPECTIVE: FROM THE 

OUTSIDE LOOKING IN 

5.1. The Milky Way as a Star-forming Galaxy 

As our nearest example, the Milky Way (hereafter de- 
noted as MW) has obvious advantages in studies of star 
formation in galaxies. However, living inside the MW 
presents serious problems of distance determination and 
selection effects, compared to studies of other galaxies. 
We can "look in from the outside" only with the aid 
of models. After reviewing surveys briefly, we consider 
properties of gas and star formation within the MW, 
bearing in mind the issues raised in §2.4 and §3, first 
as a whole, and then the radial distribution, and finally 
some notes on non-axisymmetric structure. 



5.1.1. Surveys 

Recent and ongoing surveys of the Galactic Plane at 
multiple wavelengths are revitalizing the study of the 
MW as a galaxy. Surveys of H I from both northern (Tay- 
lor et al. 2003) and southern (McClure-Griffiths et al. 
2005) hemispheres, along with a finer resolution survey 
of parts of the Galactic Plane (Stil et al. 2006) have given 
a much clearer picture of the atomic gas in the MW (for a 
review of H I surveys, see Kalberla & Kerp 2009). Based 
on these siirvcys, properties of the cool atomic clouds 
have been analyzed by Dickey et al. (2003). 

Numerous surveys of the MW have been obtained in 
CO J = 1 (e.g.. Dame et al. 2001; Bronfman et al. 
1988; Clemens et al. 1988). A survey of the inner part of 
the MW in ^^CO (Jackson ct al. 2006) has helped with 
some of problems caused by optical depth in the main 
isotopologue. 

Surveys for tracers of star formation have been made 
in radio continuum (free-free) emission (Altenhoff et al. 
1970) and recombination lines (Anderson et al. 2011; 
Lockman 1989), in water masers (Cesaroni et al. 1988; 
Walsh ct al. 2011), and methanol masers (Pestalozzi ct al. 
2005; Green et al. 2009). In addition to being signposts 
of (mostly massive) star formation, the masers provide 
targets for astrometric studies using VLBI. 

Complementary surveys of clouds have been done using 
mid-infrared extinction (e.g., Perault et al. 1996; Egan 
et al. 1998; Peretto & Fuller 2009), which can identify In- 
frared Dark Clouds (IRDCs) against the Galactic back- 
ground emission. The MIPSGAL survey (Carey et al. 
2009) should also provide a catalog of very opaque ob- 
jects. At longer wavelengths, the dust is usually in emis- 
sion, and millimeter continuum emission from dust (e.g., 
Aguirre et al. 2011; Rosolowsky ct al. 2010; Schuller et al. 
2009) provides a different sample of objects. 

Each method has its own selection effects and efforts 
are ongoing to bring these into a common framework. A 
key ingredient is to determine the distances, using spec- 
tral lines and methods to break the distance ambiguity in 
the inner galaxy; initial work is underway (Russeil et al. 
2011; Schlingman et al. 2011; Dunham et al. 2011b; Fos- 
ter et al. 2011). Ultimately, these studies should lead to a 
better definition of the total amount and distribution of 
dense gas (§1.2). The Herschel HIGAL survey (Molinari 
et al. 2010) will add wavelengths into the far-infrared 
and lead to a far more complete picture. Quantities like 
LpiR/Mcioud and Lpm/Maense, as used in extragalactic 
studies, can be calculated for large samples. 

5.1.2. The Milky Way as a Whole 

The Milky Way is a barred spiral galaxy (Burton 1988; 

Dame et al. 2001; Benjamin et al. 2005). The number and 
position of spiral arms are still topics for debate, but the 
best currently available data favor a grand-design, two- 
armed, barred spiral with several secondary arms. A con- 
ception of what the MW would look like from the outside 
(Churchwell et al. 2009) is shown in Figure 5. We can 
expect continuing improvements in the model of the MW 
as VLBI astrometry improves distance determinations of 
star forming regions across the MW (Reid et al. 2009; 
Brunthaler et al. 2011). As we will be comparing the 
MW to NGC 6946, we also show images of that galaxy 
at Ha, 24 fj,m, H I, and CO J = 2 ^ 1 in Fig. 6. 
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Based on a model including dark matter (Kalberla & 
Kerp 2009; Kalberla et al. 2007; Kalberla & Dedes 2008), 
the MW mass within 60 kpc of the center is M{tot) = 
4.6 X IQii Mq, with M{baryon) = 9.5 x 10^° Mq. The 
total mass of atomic gas (H I plus He) is M (atomic) = 
8 X 10^ Mq, and the warm ionized medium contains 
M{WIM) = 2 X 10^ M0. The mass fraction of the HIM 
(§2.1) is negligible. With a (perhaps high) estimate for 
the molecular mass of M{mol) = 2.5 x 10^ M©, Kalberla 
& Kerp (2009) derive a gas to baryon ratio of 0.13. As 
explained in §5.1.3, for consistent comparison to the ra- 
dial distribution in NGC 6946, wc will use a constant 
X{CO) = 2.0 X 10^°, correct for helium, and assume an 
outer radius of the star forming disk of 13.5 kpc to define 
masses, surface densities, etc. Within Rgai = 13.5 kpc, 
M{mol) = 1.6 X 10^ Mq and M (atomic) = 5.0 x 10^ 
Mq with these conventions. 

The volume filling factor of the CNM is about 1% (Cox 
2005) and that of molecular clouds, as traced by the ^^CO 
survey (Roman-Duval et al. 2010), has been estimated 
at about 0.5% (M. Heyer, personal communication) in 
the inner galaxy, but much lower overall. Denser (n ^ 
few X 10^ cm~^) structures found in millimeter contin- 
uum surveys, roughly corresponding to clumps, appear 
to have a surface filling factor of 10~ and a volume filling 
factor of about 10^® (M. K. Dunham, personal commu- 
nication), but better estimates should be available soon. 

The star formation rate of the MW has been estimated 
from counting H II regions, which can be seen across 
the MW and extrapolating to lower mass stars (Mezger 
1987; McKee & WiUiams 1997; Murray & Rahman 2010). 
These methods average over the effective lifetime of mas- 
sive stars, about 3-10 Myr (Table 1). Estimates of M* 
from a model of the total far-infrared emission of the 
MW (e.g., Misiriotis et al. 2006) are somewhat less sen- 
sitive to the high end of the IMF (§3) and average over 
a longer time. An alternative approach, based on count- 
ing likely YSOs in the GLIMPSE survey of the Galactic 
plane (Robitaille & Whitney 2010) is much less biased to- 
ward the most massive stars, but is limited by sensitivity, 
issues of identification of YSOs, and models of extinc- 
tion. Chomiuk & Povich (2011) have recently reviewed 

all methods of computing M* for the MW and conclude 
that they are consistent with M* = 1.9 ± 0.4 M© yr~^. 
However they conclude that resolved star counts give 
SFRs that are factors of 2-3 times higher (sec §6.4 for 
further discussion). 

With a radius of active star formation of 13.5 kpc 
(§5.1.3), = 1.9 Mq yr-i yields (S(SFR)) = 

3.3x10-3 Mq yr-ikpc-2. Taking M(mol,Rgai < 
13.5) = 1.6 X 10^ Mq, tdep = 0.8 Gyr, 10 times longer 
than (tdep) in local clouds (§4.1). For all the gas in the 
(60 kpc) MW, tdep = 6.6 Gyr, or 5.5 Gyr if the WIM is 
excluded. If we attribute all non-gaseous baryons in the 
Kalberla & Kerp (2009) model (8.25 x lO^" M©) to stars 
and stellar remnants, ignore recycling, and take an age 
of 1 x 10^° yr, we would derive an average star formation 
rate of 8.25 M© yr"^, about 4 times the current rate. 

5.1.3. Radial Distributions 

The MW also offers a wide range of conditions, from 
the far outer Galaxy to the vicinity of the Galactic Cen- 



ter, for detailed study. The hurdle to overcome, espe- 
cially in the inner MW, is to assign accurate distances. 
The radial distributions of atomic and molecular gas, 
along with the star formation rate surface density (here 
in imits of M© Gyr^^ pc;^^). arc plotted in Figure 7, 
along with a similar plot for NGC 6946 (Schruba et al. 
2011). For consistency, we have used the same A'(CO) 
and correction for helium (for both molecular and atomic 
gas) as did Schruba et al. (2011). These distributions for 
the MW have considerable uncertainties, and should be 
taken with appropriate cautions, especially in the inner- 
most regions, as discussed below. Misiriotis et al. (2006) 
presents a collection of estimates of S(SFR) and a model. 
They all show a steady decline outward from a peak at 
F^gal ~ 5 kpc, but the situation in the innermost galaxy 
is unclear. The data on H II regions that forms the basis 
of most estimates are very old, and a fresh determination 
of j:(SFR)(Rgai) is needed. 

In comparison, NGC 6946 has a similar, rather flat, dis- 
tribution of atomic gas, with S (atomic) ~ 10 M© pc"^. 
The molecular gas distribution is similar to that of the 
MW, except for a clearer and stronger peak within 
Rgai = 2 kpc. The distribution of E(SFR) follows the 
molecular gas in both galaxies, more clearly so in NGC 
6946. 

We discuss the MW distributions from the outside, 
moving inward. The average surface density of molec- 
ular gas drops precipitously beyond 13.5 kpc, after a lo- 
cal peak in the Perseus arm (Heyer et al. 1998), as does 
the stellar density (Ruphy et al. 1996), thereby defining 
the far outer Galaxy, and the radius used in our whole 
Galaxy quantities (§5.1.2). The atomic surface density 
also begins to drop around 13-17 kpc, following an ex- 
ponential with scale length 3.74 kpc (Kalberla & Kerp 
2009). Molecular clouds in the far outer galaxy are rare, 
but can be found with large-scale surveys. Star formation 
does continue in the rare molecular clouds (Wouterloot & 
Brand 1989), and Figure 7 would suggest a higher ratio 
of star formation to molecular gas, but there are issues of 
incompleteness in the outer galaxy. A study of individual 
regions (Snell et al. 2002) found a value of Lpjji/M sim- 
ilar to that for the inner Galaxy. They conclude that the 
star formation process within a cloud is not distinguish- 
able from that in the inner Galaxy; the low global rate 
of star formation is set by the inefficient conversion of 
atomic to molecular gas in the far outer Galaxy. This re- 
sult is consistent with the results of Schruba et al. (2011) 
for other galaxies, as exemplified by NGC 6946. 

Inside about 13 kpc, the atomic surface density is 
roughly constant at 10-15 M© pc~^, while the average 
molecular surface density increases sharply, passing 1 
M© pc~^ somewhere near the solar neighborhood (Dame 
et al. 2001), to a local maximum around Rgai = 4-5 kpc 
(Nakanishi & Sofue 2006). After a small decrease inside 
4 kpc, the surface density rises sharply within Rgai = 1 
kpc. Nakanishi & Sofue (2006) and thus Figure 7 use a 
constant X(CO), which could overestimate the molecular 
mass near the center (§2.4). 

The innermost part of the MW, Rgai < 250 pc, known 
as the Central Molecular Zone (CMZ), potentially pro- 
vides an opportunity for a close look at conditions in 
galactic nuclei without current AGN activity, provided 
that the issues of non-circular motions, foreground and 
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background confusion, and possible changes in X{CO) 
(Oka et al. 1998) can be overcome. Despite their conclu- 
sion that X(CO) is lower, Oka et al. (1998) argue that the 
CMZ has a molecular mass of 2-6 x lO'^ Mq . Analysis of 
dust emission yields a mass of 5 x 10^ Mq (Pierce-Price 
et al. 2000), which leads to Sga^ = 250 Mq pc'^ for the 
CMZ. A TIR luminosity of 4 x 10^ Lq (Sodroski et al. 
1997) for the CMZ would imply Af* - 0.06 M© yr'^ us- 
ing the conversion in Table 1, roughly consistent with 
other recent estimates (Immer et al. 2012), based on 
analysis of point infrared sources, but new data from 
the surveys mentioned above should allow refinement of 
these numbers. Evidence for variations in star forma- 
tion rate on short timescales is discussed by Yusef-Zadeh 
et al. (2009). Using = 0.06 Mq yi'^ and Rgai = 0.25 
kpc for the CMZ, S(SFR) = 300 Mq Gyr"! pc'^. These 
values of T,moi and S(SFR) for the CMZ are plotted sep- 
arately in Figure 7 and identified as "CMZ." If correct, 
the CMZ of the MW begins to look a bit more like the 
inner few kpc of NGC 6946. The CMZ provides an ex- 
cellent place to test scaling relations, including those for 
dense gas, in detail, if the complicating issues can be 
understood. Preliminary results suggest that S(SFR) is 
similar to that expected from T,moi (see Fig. 7), but 
lower than expected from dense gas relations (S. Long- 
more, personal communcation) . 

The properties of clouds and clumps also may vary 
with Rgai- Heyer et al. (2009) reanalyzed cloud proper- 
ties in the inner Galaxy based on CO and ^^CO sur- 
veys. For the area of the cloud defined roughly by 
the 1 K CO detection threshold, he found a median 
^mol{cloud) = 42 Mq pc~^ for clouds, substantially less 
than originally found by Solomon et al. (1987), and only 
a few clouds have Y,moi{doud) > 100 Mq pc^^. Analysis 
based on detection thresholds for ^'^CO yield a median 
S™oi(^'^CO clump) = 144 Mq pc^^ for "i3C0 clumps" 
(Roman-Duval et al. 2010, who refer to them, however as 
clouds). The mean in 0.5 kpc bins of surface density per 
^^CO clump is roughly constant at 180 M© pc~^ from 
Rgai = 3 to 6.6 kpc, beyond which it drops sharply to- 
ward the much lower values in the solar neighborhood. 
Assuming the region sampled is representative, Roman- 
Duval et al. (2010) have plotted the azimuthally averaged 
surface density of ^^CO clumps versus Galactocentric ra- 
dius; it peaks at 2.5 M© pc^^ at 4.5 kpc, declining to 
< 0.5 M© pc^^ beyond 6.5 kpc. The clouds do appear 
to be associated with spiral arms, so azimuthal averaging 
should be taken with caution. The cloud mass function of 
246 clouds in the far outer galaxy has a power-law slope 
of dcioud = —1.88, similar to but slightly steeper than 
that found inside the solar circle (§2.5) and a maximum 
mass of about 10^ M© (Snell et al. 2002). The lower 
value for the maximum mass of a cloud seems to result 
in a concomitant limit on the number of stars formed 
in a cluster, but no change is inferred for the intrinsic 
IMF (Casassus et al. 2000). A piece of a spiral arm at 
Rgai = 15 kpc has been recently identified in H I and 
CO, containing a molecular cloud with M = 5 x 10** M© 
(Dame & Thaddeus 2011). 

5.1.4. Non-axisymmetric Structure 



The molecular gas surface density shows a strong local 
maximum around Rgai = 4.5 kpc in the northern surveys 
(Galactic quadrants I and II), while southern surveys 
(quadrants III and IV) show a relatively flat distribution 
of T,moi from about 2 to 7 kpc (Bronfman et al. 2000). 
This asymmetry is likely associated with spiral struc- 
ture (Nakanishi & Sofue 2006) and a long (half-length 
of 4.4 kpc) bar (Benjamin et al. 2005). The star for- 
mation rate, measured from LpjR, does not, however, 
reflect this asymmetry; in fact, it is larger in quadrant 
IV than in I, suggesting a star formation efficiency up 
to twice as high, despite the lower peak ^moi (Bronf- 
man et al. 2000). Interestingly, the distribution of denser 
gas, traced by millimeter-wave dust continuum emission 
(Beuther et al. 2012), is also symmetric on kpc scales, 
but asymmetries appear within the CMZ region (Bally 
et al. 2010). Major concentrations of molecular clouds 
with very active star formation appear to be associated 
with regions of low shear (Luna et al. 2006), and with 
the junction of the bar and spiral arms (Nguyen Luong 
et al. 2011). 

5.2. Demographics of Star-Forming Galaxies Today 

The recent influx of multi-wavelength data has ex- 
panded the richness of information available on global 
star formation properties of galaxies, and transfornicxl 
the interpretive framework from one based on mor- 
phological types to a quantitative foundation based on 
galaxy luminosities, masses, and other physical proper- 
ties. 

The role of galaxy mass as a fundamental determinant 
of the star formation history of a galaxy has been long 
recognized (e.g., Gavazzi & Scodeggio 1996), but data 
from SDSS and subsequent surveys have reshaped our 
picture of the population of star-forming galaxies (e.g., 
Kauffmann et al. 2003; Baldry et al. 2004, 2006; Brinch- 
mann et al. 2004). The integrated colors of galaxies, 
which are sensitive to their star formation histories, show 
a strongly bimodal dependence on stellar mass, with 
a relatively tight "red sequence" populated by galaxies 
with little or no current star formation, and a somewhat 
broader "blue sequence" or "blue cloud" of actively star- 
forming galaxies. The dominant population shifts from 
blue to red near a transition stellar mass of ~ 3 x 10^° M© 
(Kauffmann et al. 2003). The relative dearth of galaxies 
in the "green valley" between the red and blue sequences 
suggests a deeper underlying physical bimodality in the 
galaxy population and a rapid evolution of galax;ies from 
blue to red sequences. 

A similar bimodality characterizes the mass depen- 
dence of the SFR per unit galaxy mass (SSFR). The mass 
dependence of the SSFR has been explored by numerous 
investigators (e.g., Brinchmann et al. 2004; Salim et al. 
2007; Schiminovich et al. 2007; Lee et al. 2007). Fig- 
ure 8 shows an example from Schiminovich et al. (2007), 
based in this case on dust-corrected FUV measurements 
of SDSS galaxies. 

A clear separation between the blue and red sequences 
is evident, and the dispersion of SSFRs within the blue 
sequence is surprisingly small, suggesting that some kind 
of self-regulation mechanism may be at work among the 
actively star-forming galaxies. The bimodality is not ab- 
solute; there is a clear tail of less active but significantly 
star forming galaxies between the two sequences. These 
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represent a combination of relatively inactive (usually 
early-type) disk galaxies and unusually active spheroid- 
dominated systems (§5.4). The sharp increase in the frac- 
tion of inactive galaxies above stellar masses of order a 
few times 10^'^ Mq is also seen. 

The blue sequence in Figure 8 is not horizontal; the 
SSFR clearly increases with decreasing galaxy mass. The 
slope (—0.36 for the data in Figure 8) varies somewhat 
between different studies, possibly reflecting the effects 
of different sampling biases. The negative slope implies 
that lower-mass galaxies are forming a relatively higher 
fraction of their stellar mass today, and thus must have 
formed relatively fewer of their stars (compared to more 
massive galaxies) in the past. The most straightforward 
explanation is that the dominant star-forming galaxy 
population in the Universe has gradually migrated from 
more massive to less massive galaxies over cosmic time. 
Direct evidence for this "downsizing" is seen in obser- 
vations of the SSFR vs mass relation in high-redshift 
galaxies (e.g., Noeske et al. 2007). 

Another instructive way to examine the statistical 
properties of star-forming galaxies is to compare abso- 
lute SFRs and SFRs normalized by mass or area. Sev- 
eral interesting trends can be seen in Figure 9, which 
plots integrated measurements of the SFR per unit area 
as a function of the absolute SFRs. The first is the ex- 
traordinary range in SFRs, more than seven orders of 
magnitude, whether measured in absolute terms or nor- 
malized per unit area or (not shown) galaxy mass. Much 
of this range is contributed by non-equilibrium systems 
(starbursts). Normal galaxies occupy a relatively tight 
range of SFRs per unit area, reminiscent of the tightness 
of the blue sequence when expressed in terms of SSFRs. 
The total SFRs of the quiescent star-forming galaxies 
are also tightly bounded below a value of ~20Mq yr~^. 
Starburst galaxies, and infrared-luminous and ultralumi- 
nous systems in particular, comprise most of the galaxies 
in the upper 2-3 decades of absolute SFRs and 5](SFR). 

The distribution of SFRs along the X-axis of Figure 9 
(after correcting for volume completeness biases) is sim- 
ply the SFR distribution function (e.g., Gallego et al. 
1995; Martin et al. 2005b). Figure 10 shows a recent 
determination of this distribution from Bothwell et al. 
(2011), based on flux-limited UV and TIR samples of 
galaxies and SFRs derived using the methods of Hao 
et al. (2011), and corrections for AGN contamination 
following Wu et al. (2010b, 2011). The SFR function is 
well fitted by a Schechter (1976) exponentially truncated 
power-law, with a faint-end slope a — —1.5 and char- 
acteristic SFR, SFR* = 9Mq. Although SFRs as high 
as ~1000M© yr~^ are found in present-day ULIRGs, the 
contribution of LIRGs and ULIRGs to the aggregate star 
formation today is small (< 10% also see Goto et al. 
2011). Steady-state star formation in galaxies dominates 
today. The value of SFR* increases rapidly with z, and 
galaxies with Ltoi > 10^^ Lq (i.e., LIRGs) become domi- 
nant by redshifts z > 1 (e.g., Le Floc'h et al. 2005). This 
change partly reflects an increase in merger-driven star 
formation at higher redshift, but at early epochs even 
the steady-state star formation in massive galaxies at- 
tained levels of order tens to hundreds of solar masses per 
year, thus placing those galaxies in the LIRG and ULIRG 
regime. This suggests that most of the decrease in the 
cosmic SFR in recent epochs has been driven by downsiz- 



ing in the level of steady-state star formation (e.g., Bell 
et al. 2005; Jogee et al. 2009). 

5.3. The High-Density Regime: Starbursts 

As highlighted in §1.2, the starburst phenomenon en- 
compasses a wide range of physical scales and host galaxy 
properties, and no precise physical definition of a star- 
burst has been placed into wide use. Part of the expla- 
nation can be seen in Fig. 9; although the most extreme 
starburst galaxies have properties that are well separated 
from those of normal star- forming galaxies, there is no 
clear physical break in properties; starbursts instead de- 
fine the upper tails of the overall distributions in SFRs 
within the general galaxy population. 

Much of the attention in this area continues to focus 
on the infrared-luminous and ultraluminous systems, be- 
cause they probe star formation in the most extreme 
high-density circumnuclear environments seen in the lo- 
cal Universe. A major breakthrough in recent years 
has been the use of mid-infrared spectroscopy to distin- 
guish dust heated by massive stars from that heated by a 
buried AGN (e.g., Genzel et al. 1998; Laurent et al. 2000; 
Dale et al. 2006; Armus et al. 2007). Regions heated by 
AGNs arc distinguished by the appearance of highly ion- 
ized atomic species, as well as suppressed mid-infrared 
PAH emission relative to stellar- heated dust. Applica- 
tion of these diagnostics has made it possible to construct 
clean samples of LIRGs and ULIRGs that are dominated 
by star formation. 

With the Spitzer and Herschel observatories, it has 
been possible to extend imaging and spectroscopy of 
the most luminous infrared-emitting galaxies to inter- 
mediate redshifts (e.g., Elbaz et al. 2005; Yan et al. 
2007). As mentioned earlier the fraction of star for- 
mation in LIRGs and ULIRGs increases sharply with 
redshift, but imaging and spectroscopy of these objects 
reveals a marked shift in the physical characteristics of 
the population. Whereas at low redshift, ULIRGs (and 
many LIRGs) are dominated by very compact circum- 
nuclear starbursts triggered by mergers, at higher red- 
shifts the LIRG/ULIRG population becomes increasingly 
dominated by large disk galaxies with extended star for- 
mation. As a consequence of cosmic downsizing (i.e., 
"upsizing" with increasing redshift), by redshifts z ^ 1 
the populations of LIRGs and ULIRGs become increas- 
ingly dominated by the progenitors of present-day nor- 
mal galaxies, rather than by the transient merger-driven 
starbursts which dominate the present-day popiilations 
of ULIRGs, though examples of the latter are still found, 
especially in the population of submillimeter-luminous 
galaxies (SMGs) (e.g.. Chapman et al. 2005). 

Visible-wavelength observations of these high-redshift 
galaxies with IFU instruments, both with and without 
adaptive optics correction, have revealed many insights 
into the physical nature of this population of starburst 
galaxies (e.g., Genzel et al. 2008, 2011; Forster Schreiber 
et al. 2009 and references therein). The Ha kinemat- 
ics show a wide range of properties, from normal, dif- 
ferentially rotating disks to disturbed disks, and a sub- 
set with signatures of ongoing or recent mergers. The 
disks tend to be characterized by unusually high velocity 
dispersions, which have been interpreted as reflecting a 
combination of dynamical instability and possibly energy 
injection into the ISM from young stars. Further discus- 
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sion of high-rcdshift galaxies is beyond the scope of this 
article, but wc shall return to the cold ISM properties 
and SFRs of these galaxies in §6. 

5.4. The Low SFR and Low Density Regimes 

One of the most important discoveries from the 
GALEX mission was the detection of low levels of star 
formation in environments which were often thought 
to have been devoid of star formation. These include 
early-type galaxies, dwarf galaxies, low surface bright- 
ness galaxies, and the extreme outer disks of many nor- 
mal galaxies. 

Although star formation had been detected occasion- 
ally in nearby elliptical and SO galaxies (e.g., Pogge & 
Eskridge 1993), most of these galaxies have historically 
been regarded as being "red and dead" in terms of recent 
star formation. Deep GALEX imaging of E/SO galaxies 
(e.g., Kaviraj et al. 2007), however, has revealed that 
approximately 30% of these early-type galaxies exhibit 
near-ultraviolet emission in excess of what could reason- 
ably arise from an evolved stellar population (e.g., Kavi- 
raj et al. 2007). Confirmation of the star formation has 
come from visible-wavelength IFU spectroscopy, as ex- 
ploited for example by the SAURON survey (Shapiro 
et al. 2010). The SFRs in these galaxies tend to be very 
low, with at most 1 3% of the stellar mass formed over 
the past Gyr. High-resolution UV imaging with HST of- 
ten reveals extended star-forming rings or spiral arms in 
these galaxies (Salim & Rich 2010). Follow-up CO ob- 
servations of these galaxies reveals significant detections 
of molecular gas in a large fraction of the galaxies with 
detected star formation, with molecular gas masses of 
10^ - 10^° Mq (Combes et al. 2007; Crocker et al. 2011), 
and with SFRs roughly consistent with the Schmidt law 
seen in normal galaxies. 

The GALEX images also reveal that star formation 
in more gas-rich disk galaxies often extends much far- 
ther in radius than had previously been appreciated. In 
exceptional cases (e.g., NGC5236=M83, NGC4625) the 
star formation extends to 3-4 times the normal R25 op- 
tical radius, to the edge of the HI disk (Thilker et al. 
2005; Gil de Paz et al. 2005). A systematic study by 
Thilker et al. (2007) reveals that extended "XUV disks" 
are found in ~20% of spiral galaxies, with less distinct 
outer UV structures seen in another ~10% of disk galax- 
ies. Follow-up deep Hq imaging and/or spectroscopy re- 
veals extended disks of HII regions that trace the UV 
emission in most cases (e.g., Zaritsky & Christlein 2007; 
Christlein et al. 2010; Goddard et al. 2010), confirming 
the earlier detection of HII regions at large radii (e.g., 
Ferguson et al. 1998; Ryan- Weber et al. 2004). Individ- 
ual examples of extended UV features without Ha coun- 
terparts arc sometimes found, however (Goddard et al. 
2010). The GALEX images also have led to a break- 
through in measurements of star formation in low sur- 
face brightness spiral galaxies (LSBs). These provide a 
homogeneous body of deep measurements of the star for- 
mation and the nature of the star formation law at low 
surface densities (§6; Wyder et al. 2009). Comparison 
to the outer MW (§5.1) suggests that deep searches for 
molecular gas will find few cloiids, which will nonetheless, 
be the sites of star formation. 

The other low-density star formation regime can be 
found in dwarf irregular galaxies. Drawing general in- 



ferences about the star formation properties of these 
galaxies requires large samples with well-defined (ide- 
ally volume-limited) selection criteria, and studies of this 
kind (based on UV, visible, and/or Ha observations) 
have been carried out by several groups (e.g., Karachent- 
sev et al. 2004; Blanton et al. 2005; Meurer et al. 2006; 
Kennicutt et al. 2008; Dalcanton et al. 2009; Hunter et al. 

2010) . One perhaps surprising result is the near ubiq- 
uity of star formation in the dwarf galaxies. Kennicutt 
et al. (2008) observed or compiled Ha luminosities for 
galaxies within the local 11 Mpc, which included ~300 
dwarf galaxies (Mb > —17). Excluding a handful of 
dwarf spheroidal galaxies which have no cold gas, only 
10 of these (^^3%) were not detected in Ha, meaning that 
star formation has taken place over the last 3-5 Myr in 
the other 97% of the systems. Moreover many of the 
handful of Ha non-detections show knots of UV emis- 
sion, demonstrating that even fewer of the galaxies have 
failed to form stars over the last ^100 Myr (Lee et al. 

2011) . Recent star formation is seen in all of the galaxies 
with Mb < -13 and M{HI) > 5 x lO'^M©. For less 
massive galaxies it is possible for star formation to cease 
for timescales that are longer than the ionization lifetime 
of an HII region (up to 5 Myr) , but examples of galaxies 
with extended periods of no star formation are extremely 
rare. 

These data also provide a fresh look at the temporal 
properties of the star formation in dwarf galaxies gener- 
ally. The distribution of SSFRs shows a marked increase 
in dispersion for galaxies with Mb > —14.5 or circular ve- 
locity of 50 km s-i (Lee et al. 2007; Bothwell et al. 2009). 
The corresponding neutral gas fraction does not show a 
similar increase in dispersion (Bothwell et al. 2009), so 
this suggests an increased short-term fiuctuation of in- 
tegrated SFRs in low-mass galaxies. Lee et al. (2009b) 
analyzed the statistics of the SFRs in more depth and 
confirmed a larger fraction of stars formed in bursts in 
low- mass galaxies (~25%), but these still represent only 
a small fraction of the total stars formed- even in the 
dwarfs the majority of stars appear to form in extended 
events of duration longer than ~10 Myr. Recent mod- 
elling of the statistics of Ha and UV emission by Weisz 
et al. (2012) and Fumagalli et al. (2011) as well as analy- 
sis of resolved stellar populations (e.g., Weisz et al. 2008) 
confirm the importance of fluctuations in the SFR. 

6. STAR FORMATION RELATIONS 

The immense dispersion in SFR properties seen in Fig- 
ure 9 collapses to a remarkably tight scaling law when 
the SFR surface densities (I](SFR)) are plotted against 
mean gas surface densities C^gas) (Kennicutt 1998b). 
This emergent order reflects the fact that gas is the in- 
put driver for star formation. The concept of a power- 
law relation between SFR density and gas density dates 
to Schmidt (1959, 1963), and relations of this kind are 
commonly referred to as "Schmidt laws." On physical 
grounds we might expect the most fundamental relation 
between the volume densities of star formation and gas, 
but since most observations of external galaxies can only 
measure surface densities integrated along the line of 
sight, the most commonly used relation, often called a 
Kennicutt-Schmidt (KS) law, is in terms of surface den- 
sities. 
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^SFR = AT, gas (13) 

The precise form of this relation depends on assumptions 
about how Egas is derived from the observations (§2.4), 
but a strong correlation is clearly present. 

6.1. The Disk-Averaged Star Formation Law 

K98 presented a review of observations of the Schmidt 
law up to the time, and nearly all of that work charac- 
terised the relation between the disk-averaged SFR and 
gas surface densities in galaxies. The upper panel of Fig- 
ure 11 presents an updated version of the global Schmidt 
law in galaxies. Each point is an individual galaxy (color 
coded as explained in the caption) , with the surface den- 
sity defined as the total gas mass (molecular plus atomic) 
or SFR normalized to the radius of the main star-forming 
disk, as measured from Ha, Paa, or IR maps. For sim- 
phcity, a constant X{CO) factor [2.3 x 10^° cm'^ (K km 
s"^)^^ with no correction for helium, see §2.4] has been 
applied to all of the galaxies; the consequences of a pos- 
sible breakdown in this assumption are discussed later. 
The sample of galaxies has been enlarged from that stud- 
ied in Kennicutt (1998b), and all of the Ha-bascd SFR 
measurements have been improved by incorporating in- 
dividual (IR-based) corrections for dust attenuation and 
[N II] contamination. 

The form of this integrated Schmidt law appears to be 
surprisingly insensitive to SFR environment and param- 
eters such as the atomic versus molecular fraction, but 
some metal-poor galaxies (defined as Z < 0.3^0 ) devi- 
ate systematically from the main relation, as shown by 
the blue open circles in Figure 11 (upper panel). These 
deviations could arise from a physical change in the star 
formation law itself, but are more likely to reflect a break- 
down in the application of a constant X{CO) factor (§2.4; 
Lcroy et al. 2011, and references therein). Adopting 
higher values of X(CO) for metal-poor galaxies brings 
the galaxies much more into accord with the main rela- 
tion in Figure 11. 

Recent observations of LSBs by Wyder et al. (2009) ex- 
tend the measurements of the integrated star formation 
law to even lower mean surface densities, as shown by the 
purple crosses in the upper panel of Figure 11. A clear 
turnover is present, which is consistent with breaks seen 
in spatially-resolved observations of the star formation 
law (§6.2). 

The slope of the integrated Schmidt law is non-linear, 
with N ~ 1.4-1.5 (Kennicutt 1998b), when a constant 
X{CO) factor is applied. This uncertainty range does 
not include all possible systematic errors, arising for ex- 
ample from changes in X(CO) or the IMF with increas- 
ing siirfacc density or SFR. Major systematic changes in 
cither of these could easily change the derived value of 
N by as much as 0.2-0.3. For example, if X{CO) were 
five times lower in the dense starburst galaxies (§2.4), 
the slope of the overall Schmidt law would increase from 
1.4-1.5 to 1.7-1.9 (Narayanan et al. 2012). 

Usually the Schmidt law is parametrized in terms of 
the total (atomic plus molecular) gas siirface density, 
but one can also explore the dependences of the disk- 
averaged SFR densities on the mean atomic and molecu- 
lar surface densities individually. Among normal galaxies 
with relatively low mean surface densities, the SFR den- 



sity is not particularly well correlated with either compo- 
nent, though variations in X{CO) could partly explain 
the poor correlation between SFR and derived H2 den- 
sities (e.g., Kennicutt 1998b). In starburst galaxies with 
high gas surface densities, however, the gas is overwhelm- 
ingly molecular, and a strong non-linear Schmidt law is 
observed (upper panel of Figure 11). 

A similar non-linear dependence is observed for total 
SFR (as opposed to SFR surface density) versus total 
molecular gas mass (e.g., Solomon & Sage 1988; Gao & 
Solomon 2004), and presumably is another manifesta- 
tion of the same underlying physical correlation. The 
dependence of the SFR on dense molecular gas mass is 
markedly different, however. The lower panel of Figure 
11, taken from Gao & Solomon (2004), shows the relation 
between the integrated SFRs and the dense molecular 
gas masses, as derived from HON J = 1 — >■ measure- 
ments (§2.4) for a sample of normal and starburst galax- 
ies. In contrast to the correlation with total molecular 
mass from CO J = 1 — )> 0, this relation is linear, implying 
a strong coupling between the masses of dense molecular 
clumps and stars formed, which is largely independent of 
the galactic star-forming environment. Wu et al. (2005b) 
have subsequently shown that this linear relation extends 
down to the scales of individual star-forming molecular 
clouds and dense clumps in the Galaxy. Combined with 
the MW studies (§4), these dense gas relations for galax- 
ies suggest that dense clumps are plausible fundamental 
star-forming units. If so, the mass fraction of the ISM 
(and fraction of the total molecular gas) residing in dense 
clumps must itself increase systematically with the SFR. 

Because Lhcn in the J = 1 — > line is only an ap- 
proximate tracer for dense gas (§2.4.3) and because the 
definition of a dense clump is by no means precise (§2.5), 
one should be careful not to overinterpret these trends. 
Similar comparisons of the SFR with the emission from 
other molecular tracers suggest that the slope of the SFR 
versus gas mass relation changes continuously as one pro- 
ceeds from lower-density tracers such as CO J = 1 — >■ 
to higher-excitation (and higher-density) CO transitions, 
and further to high-density tracers such as HCN and 
HCO+ and to higher transitions of those molecules (e.g., 
Juneau et al. 2009). Multi-transition studies, along with 
realistic modeling will help to refine the interpretation of 
line luminosities of dense gas tracers (e.g., Gracia-Carpio 
et al. 2008; Juneau et al. 2009). 

The data shown in Figure 11 all come from observa- 
tions of nearby galaxies (z < 0.03), but recently a num- 
ber of studies have addressed the form of the molecular 
gas Schmidt law for starburst galaxies extending to red- 
shifts z>2 (e.g., Bouche et al. 2007; Daddi et al. 2010; 
Genzel et al. 2010). The interpretation of these results 
is strongly dependent on the assumptions made about 
A (CO) in these systems. When a Galactic A (CO) con- 
version is applied, the high-redshift galaxies tend to fall 
roughly on the upper parts of the Schmidt law seen lo- 
cally (e.g.. Figure 11). However if a lower A(CO) factor 
is applied to the most compact starburst and submillime- 
ter galaxies (SMGs), as suggested by many independent 
analyses of A(CO) (§2.4), the Schmidt relations shift 
leftwards by the same factor, forming a parallel relation 
(Daddi et al. 2010; Genzel et al. 2010). 

Taken at face value, these results suggest the pres- 
ence of two distinct modes of star formation with diff'er- 
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ent global efficiencies, which separate the extended star- 
forming disks of normal galaxies from those in the densest 
circumnuclear starbursts; this distinction is likely to be 
present both in the present-day Universe and at early cos- 
mic epochs. This inferred bimodality, however, is a direct 
consequence of the assumption of two discrete values for 
X(CO) in the two modes. A change in the interpretation 
of CO emission is certainly plaiisiblc when the derived 
molecular surface density is similar to that of an indi- 
vidual cloud (Smoi > 100 Mq pc~^), but the behavior of 
X{CO) may be complex (§2.4). Variation of X{CO) over 
a continuous range would indicate a steeper Schmidt law, 
rather than a bimodal law (Narayanan et al. 2012). 

Daddi et al. (2010) also note that the higher concentra- 
tion of gas in the ULIRGs/SMGs is consistent with an en- 
hanced fraction of dense gas and the dense gas relations 
of Gao & Solomon (2004). However, the star formation 
rate also appears to be larger for a given mass of dense 
gas in extreme starbursts (Garci'a-Burillo et al. 2012 and 
references therein) , especially when a lower value of ancN 
is used. As with CO, one can interpret these as bimodal 
relations or as steeper than linear dependences on the 
gas tracer lines. These results, together with evidence 
for lower SFR per mass of dense gas in the CMZ of the 
MW (§5.1) warn against an overly simplistic picture of 
dense clumps as the linear building blocks for massive 
star formation, with no other variables in the picture. 

The correlations between SFR and gas surface densi- 
ties (and masses) are not the only scaling laws that are 
observed. Kennicutt (1998b) pointed out that the SFR 
surface densities also correlate tightly with the ratio of 
the gas surface density to the local dynamical time, de- 
fined in that case to be the average orbit time. This pre- 
scription is especially useful for numerical simulations 
and semi-analytical models of galaxy evolution. Inter- 
estingly, Daddi et al. (2010) and Genzel et al. (2010) 
found that the bifurcation of Schmidt laws between nor- 
mal galaxies and ULIRGs/SMGs described above does 
not arise in the dynamical form of this relation. 

Blitz & Rosolowsky (2006) have discovered another 
strong scaling relation between the ratio of molecular to 
atomic hydrogen in disks and the local hydrostatic pres- 
sure. The relation extends over nearly three orders of 
magnitude in pressure and H2 /H I ratio and is nearly lin- 
ear (slope = 0.92). Technically speaking this scaling re- 
lation only applies to the phase balance of cold gas rather 
than the SFR, but it can be recast into a predicted star 
formation law if assumptions are made about the scal- 
ing between the SFR and the molecular gas components 
(e.g.. Blitz & Rosolowsky 2006; Leroy et al. 2008). Re- 
cent work by Ostriker et al. (2010) and Ostriker & Shetty 
(2011) provides a theoretical explanation for these rela- 
tions. 

Finally, a number of workers have explored the scal- 
ing between SFR surface density and a combination of 
gas and stellar surface densities. For example Dopita 
(1985) and Dopita & Ryder (1994) proposed a scaling 
between the SFR density and the product of gaseous 
and stellar surface densities; the latter scales with the 
disk hydrostatic pressure and hence bears some relation 
to the picture of Blitz & Rosolowsky (2006). More re- 
cently Shi et al. (2011) show that the scatter in the star 
formation law is minimized with a relation of the form 

^SFR OC llgas 



Determining which of these different formulations of 
the star formation law is physical and which are mere 
consequences of a more fundamental relation is difficult 
to determine from observations alone. Some of the degen- 
eracies between these various relations can be understood 
if most gas disks lie near the limit of gravitational sta- 
bility (Q'-^l; Kennicutt 1989), and the critical column 
densities for the formation of cold gas phase, molecule 
formation, and gravitational instability lie close to each 
other (e.g., Elmegreen & Parravano 1994; Schaye 2004). 
In such conditions it can be especially difficult to identify 
which physical process is most important from observa- 
tions. We return to this topic later. 

We conclude this section by mentioning that there arc 
other SFR scaling laws that can be understood as aris- 
ing from an underlying Schmidt law. The best known of 
these is a strong correlation between characteristic dust 
attenuation in a star-forming galaxy and the SFR itself, 
with the consequence that galaxies with the highest abso- 
lute SFRs arc nearly all dusty infrared-luminous and ul- 
traluminous galaxies (e.g., Wang & Heckman 1996; Mar- 
tin et al. 2005b; Bothwell et al. 2011). This opacity ver- 
sus SFR relation is partly a manifestation of the Schmidt 
law, because we now know that the most intense star 
formation in galaxies takes place in regions with abnor- 
mally high gas surface densities, and thus also in regions 
with abnormally high dust surface densities. The other 
factor underlying the SFR versus opacity correlation is 
the prevalence of highly concentrated circumnuclear star 
formation in the most intense starbursts observed in the 
present-day universe (§5.3); this may not necessarily be 
the case for starburst galaxies at early cosmic epochs. 

6.2. Radial Distributions of Star Formation and Gas 

Over the past few years major progress has been 
made in characterizing the spatially-resolved star for- 
mation law within individual galaxies. This work has 
been enabled by large multi- wavelength surveys of nearby 
galaxies such as the Spitzer Infrared Nearby Galaxies 
Survey (SINGS) (Kennicutt et al. 2003), the GALEX 
Nearby Galaxies Survey (Gil de Paz et al. 2007), the 
Spitzer/GALEX Local Volume Legacy survey (LVL; 
Dale et al. 2009; Lee et al. 2011), and the Herschel 
KINGFISH Survey (Kennicutt et al. 2011). The re- 
sulting datascts provide the means to measure spatially- 
resolved and dust-corrected SFRs across a wide range 
of galaxy properties (§3). These surveys in turn have 
led to large spinoff surveys in H I (e.g., THINGS, Wal- 
ter et al. 2008; FIGGS, Begum et al. 2008; Local Vol- 
ume HI Survey; Koribalski 2010; LITTLE THINGS, 
Hunter et al. 2007) and in CO (e.g., BIMA-SONG, Heifer 
et al. 2003; IRAM HERACLES, Leroy et al. 2009; JCMT 
NGLS, Wilson et al. 2009; CARMA/NRO Survey, Koda 
& Nearby Galaxies CO Survey Group 2009; CARMA 
STING, Rahman et al. 2011, 2012). 

The capabilities of these new datasets arc illustrated in 
Figure 6, which shows Spitzer 24 fj,m, Ha, VLA H I, and 
CO (in this case from HERACLES) maps of NGC 6946 
(see Figure 7 for the corresponding radial distributions 
of gas and star formation) . 

The first step in exploiting the spatial resolution of the 
new observations is to analyze the azimuthally- averaged 
radial profiles of the SFR and gas components (as illus- 
trated for example in Figure 7) and the resulting SFR 



Star Formation 



27 



versus gas surface density correlations. This approach 
has the advantage of spatially averaging over large phys- 
ical areas, which helps to avoid the systematic effects 
that are introduced on smaller spatial scales (§3.9). Ra- 
dial profiles also have limitations arising from the fact 
that a single radial point represents an average over sub- 
regions with often wildly varying local gas and SFR den- 
sities, and changes in other radially varying physical pa- 
rameters may be embedded in the derived SFR versus 
gas density relations. 

Large-scale analyses of the star formation law derived 
in this way have been carried out by numerous authors 
(e.g., Kcnnicutt 1989; Martin & Kcnnicutt 2001; Wong 
& Blitz 2002; Boissier ct al. 2003; Hcyer et al. 2004; Ko- 
mugi et al. 2005; Schuster et al. 2007; Leroy et al. 2008; 
Schruba et al. 2011; Gratier et al. 2010). A strong cor- 
relation between SFR and gas surface density is seen in 
nearly all cases, with a non- linear slope when plotted in 
terms of total (atomic -|-molecular) density. The best fit- 
ting indices N vary widely, however, ranging from 1.4-3.1 
for the dependence on total gas density and 1.0-1.4 for 
the dependence on molecular gas surface density alone 
(for constant X{CO)). The measurements of the MW 
range from = 1.2 ± 0.2, when only molecular gas is 
used (Luna et al. 2006), to 2.18±0.20 when total (atomic 
and molecular) gas is used (Misiriotis et al. 2006). 

Some of the differences between these results can be 
attributed to different schemes for treating dust atten- 
uation, different CO line tracers, different metallicities 
(which may influence the choice of X(CO)), and different 
fitting methods. The influence of a metallicity-dcpcndcnt 
X{CO) factor is investigated explicitly by Boissier et al. 
(2003). Since any variation in X{CO) is likely to flat- 
ten the radial molecular density proflles, a variable con- 
version factor tends to steepen the slope of the derived 
Schmidt law; for the prescription they use, the best flt- 
ting slope can increase to as high as iV = 2.1 — 3.6, 
underscoring once again the critical role that assump- 
tions about X(CO) play in the empirical determination 
of the star formation law. Despite these differences in 
methodology, real physical variation in the SFR versus 
gas density relation on these scales cannot be ruled out. 

This body of work also conflrms the presence of a 
turnover or threshold in the star formation relation at 
low surface densities in many galaxies. Early work on 
this problem (e.g., Kcnnicutt 1989; Martin & Kcnnicutt 
2001; Boissier et al. 2003) was based on radial proflles in 
Ha, and some questions have been raised about whether 
these thresholds resulted from breakdowns in the SFR 
versus Ha calibration in low surface brightness regimes, 
as opposed to real thresholds in the SFR. Subsequent 
comparisons of UV and Ha profiles appear to confirm 
the presence of radial turnovers in the SFR in most disks, 
however, even in cases where lower levels of star forma- 
tion persist to much larger radii (e.g., Thilker et al. 2007; 
Christlcin et al. 2010; Goddard et al. 2010). This ex- 
tended star formation can be seen in Figure 7 (i? > 7 — 8 
kpc). Recent studies of nearby galaxies with unusually 
extended XUV disks also show that the UV-based SFR 
falls off much more rapidly than the cold gas surface den- 
sity, with global star formation efficiencies (e') an order 
of magnitude or more lower than those in the inner disks 
(Bigiel et al. 2010). Taken together, these recent results 
and studies of the outer MW (§5.1) confirm the presence 



of a pronounced turnover in the Schmidt law for total 
gas at surface densities of order a few Mq pc~^ (Figure 
12). 

6.3. The Star Formation Law on Sub-Kiloparsec Scales 

The same multi-wavelength data can be used in prin- 
ciple to extend this approach to point-by-point studies 
of the star formation law. The angular resolution of the 
currently available datasets (optical/UV, mid-far IR, H I, 
CO) allows for extending this analysis down to angular 
scales of ~10", which corresponds to linear scales of 30- 
50 pc in the Local group and 200- 1000 pc for the galaxies 
in the local supercluster targeted by SINGS, THINGS, 
KINGFISH, and similar surveys. Recent analyses have 
become available (Kcnnicutt et al. 2007; Bigiel et al. 
2008, 2010, 2011; Leroy et al. 2008; Blanc et al. 2009; 
Bales et al. 2010; Verley et al. 2010; Liu et al. 2011; Rah- 
man et al. 2011, 2012; Schruba et al. 2011), and several 
other major studies arc ongoing. 

The most comprehensive attack on this problem to 
date is based on a combination of SINGS, THINGS, and 
HERACLES observations (papers above by Bigiel, Leroy, 
and Schruba), and Figure 12, taken from Bigiel et al. 
(2008), nicely encapsulates the main results from this 
series of studies. The colored regions are the loci of indi- 
vidual sub-kiloparsec measurements of SFR surface den- 
sities (measured from a combination of FUV and 24 /xm 
infrared fiuxes) and total gas densities (H I plus H2 from 
C0(2-l) with a constant X(CO) factor), for 18 galaxies 
in the SINGS/THINGS sample. Other data (including 
Kcnnicutt 1998b) arc ovcrplotted as described in the fig- 
ure legend and caption. The measurements show at least 
two distinct regimes, a low-density sub-threshold regime 
where the dependence of the SFR density on gas density 
is very steep (or uncorrelated with gas density), and a 
higher density regime where the SFR is strongly corre- 
lated with the gas density. The general character of this 
relation, with a high-density power law and a low-density 
threshold confirms most of the results presented earlier. 

One can also examine separately the dependence of the 
SFR density on the atomic and molecular surface densi- 
ties and the two relations arc entirely different. The local 
SFR density is virtually uncorrelated with H I density, in 
part because the range of H I column densities is trun- 
cated above Shi ~ lOMgpc'^ or iV(H I) ~ 10^^ cm-^. 
This upper limit to the local H I densities corresponds to 
the column density where the H I efliciently converts to 
molecular form. This phase transition also roughly coin- 
cides with the turnover in the SFRs in Figure 12, which 
hints that the SFR threshold itself may be driven in part, 
if not entirely, by an atomic-molecular phase transition 
(see §7). 

In contrast to the H I, the SFR surface density is 
tightly correlated with the H2 surface density, as inferred 
from CO. A recent stacking analysis of the HERACLES 
CO maps made it possible to statistically extend this 
comparison to low H2 column densities, and it shows 
a tight correlation between SFR and CO surface bright- 
ness that extends into the H I-dominated (sub-threshold) 
regime (Schruba ct al. 2011), reminiscent of studies in the 
outer MW (Fig. 7, §5.1). 

The strong local correlation between I](SFR) and T,„ioi 
is seen in all of the recent studies which probe linear 
scales of '~200pc and larger. This is hardly surpris- 
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ing, given the strong coupling of star formation to dense 
molecular gas in local clouds (§4, §6.4). However, on 
smaller linear scales the measurements show increased 
scatter for reasons discussed in §3.9. On these scales we 
expect the scaling laws to break down, as the stars and 
gas may arise from separate regions; this effect has been 
directly observed in high-resolution observations of M33 
(Onodera et al. 2010). 

The main results described above - the presence of 
a power-law SFR relation with a low-density threshold, 
the lack of correlation of the local SFR with H I column 
density, the strong correlation with H2 column density, 
and the rapidly increasing scatter in the Schmidt law 
on linear scales below 100 200 pc arc seen consistcaitly 
across most, if not all, recent spatially-resolved studies. 
Most of the recent studies also derive a mildly non-linear 
slope to the SFR density versus total gas density Schmidt 
law in the high-density regime, with indices N falling 
in the range 1.2-1.6. Less clear from the recent work 
is the linearity and slope of the Schmidt law on small 
scales, especially the depcndcincc of SFR surface density 
on molecular gas surface density. The results from the 
HERACLES/THINGS studies have consistently shown 
a roughly linear T,sfr versus T,moi relation (N — 1.0- 
1.1), and similar results have been reported by Blanc 
et al. (2009), Ealcs et al. (2010), and Rahman et al. 
(2012). Other authors however have reported steeper 
dependences {N = 1.2-1.7), closer to those seen in inte- 
grated measurements (e.g., Kennicutt et al. 2007; Verley 
et al. 2010; Rahman et al. 2011; Liu et al. 2011; Momose 
2012). The discrepancies between these results arc much 
larger than the estimated fitting errors, and probably 
arise in part from systematic differences in the observa- 
tions and in the way the data are analyzed (see discussion 
in e.g., Blanc et al. 2009). One effect which appears to be 
quite important is the way in which background diffuse 
emission is treated in measuring the local SFRs (§3.9). 
Liu et al. (2011) show that they can produce either a lin- 
ear or non-linear local molecular SFR law depending on 
whether or not the diffuse emission is removed. Another 
CARMA-based study by Rahman et al. (2011), however, 
suggests that the effects of diffuse emission may not be 
sufficient to account for all of the differences between dif- 
ferent analyses. Another important factor may be the ex- 
citation of the CO tracer used (e.g., Juneau et al. 2009). 

It is important to bear in mind that nearly all of our 
empirical knowledge of the form of the local star forma- 
tion law is based on observations of massive gas-rich spi- 
ral galaxies with near-solar gas-phase metal abundances. 
Most studies of the star formation relations in metal-poor 
dwarf galaxies have been limited to H I data or marginal 
detections in CO at best (e.g., Bigiel et al. 2008). A 
recent study of the SMC (- 1/5 Zq) by Bolatto ct al. 
(2011) offers clues to how these results may change in 
low-metallicity environments. They find that the atomic- 
dominated threshold regime extends up to surface densi- 
ties that are an order of magnitude higher than in spirals. 

6.4. Local Measurements of Star Formation Relations 

With the ongoing large-scale surveys of the Galaxy and 
the Magellanic Clouds in recent years it is becoming pos- 
sible to investigate star formation rate indicators, the 
scaling laws, and other star formation relations for local 
samples. These can provide valuable external checks on 



the methods applied on larger scales to external galax- 
ies and probe the star formation relations on physical 
scales that are not yet accessible for other galaxies. Local 
studies of the conversion of CO observations into column 
density or mass were discussed in §2.4. 

Comparison of various star formation rate indicators 
for the MW (§5.1), as would be used by observers in 
another galaxy, against more direct measures, such as 
young star counts, suggests that star formation rates 
based on the usual prescriptions (mid-infrared and ra- 
dio continuum) may be underestimating absolute M* 
by factors of 2-3 (Chomiuk & Povich 2011). The au- 
thors suggest that changes to the intermediate mass IMF, 
timescale issues, models for O stars, and stochastic sam- 
pling of the upper IMF can contribute to the discrep- 
ancy. As discussed earlier, measures of star formation 
rate that use ionizing photons (Ha or radio continuum) 
require regions with a fully populated IMF to be reli- 
able, and they systematically underestimate the SFR in 
small star-forming clouds (§3.9). For example, star for- 
mation in clouds near the Sun would be totally invisible 
to these measures and the star formation rate would be 
badly underestimated in the Orion cluster (§3.9). Tests 
of other star formation rate tracers, such as 24 /im omis- 
sion, within the Milky Way would be useful. Changes 
in the IMF as large as those proposed by Chomiuk & 
Povich (2011) should also be readily observable in more 
evolved Galactic star clusters, if not in the field star IMF 
itself. 

Using a YSO-counting method for five GMCs in the 
N159 and N44 regions in the LMC, Chen ct al. (2010) 
were able to reach stellar masses of about 8 Mq, below 
which they needed to extrapolate. The resulting ratios 
of M,(YSOs) to M,(Ha + 24/im) ranged from 0.37 to 
11.6, with a mean of 3.5. 

Recent studies have begun to probe the form of the 
star formation relations on scales of clouds and clumps. 
Using the YSO star counting method to get M* and ex- 
tinction maps to get mean mass surface density of indi- 
vidual, nearby clouds, Evans et al. (2009) found that the 
local clouds all lay well above the K98 relation. Taken 
in aggregate, using the mean T,gas, they lay a factor of 
20 above the Kennicutt et al. (1998b) relation and even 
farther above the relation of Bigiel et al. (2008). 

Subsequent studies show evidence for both a surface 
density threshold for star-forming clumps and a E(SFR)- 
T,moi relation in this high-density sub-cloud regime. Lada 
ct al. (2010) found a threshold surface density for efficient 
star formation in nearby clouds; the star formation rate 
per cloud mass scatters widely (Fig. 13), but is linearly 
proportional to the cloud mass above a surface density 
contour of 116 ± 28 Mq pc^^ (Fig. 14), and the coef- 
ficient agrees well with the dense gas relations of Gao 
& Solomon (2004). Heiderman et al. (2010) studied the 
behavior of S(SFR) on smaller scales using contours of 
extinction. They limited the YSOs to Class I and Flat 
SED objects to ensure that they were still closely related 
to their surrounding gas and foimd a steep increase in 
E(SFR) with increasing Sgos up to about 130 Mq pc , 
above which a turnover was suggested. By adding the 
dense clumps data from Wu et al. (2010a), they identi- 
fied a turnover at 129 ± 14 M© pc'^, where the S(SFR) 
began to match the dense gas relation, but was far above 
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the K98 relation for total gas (Figure 15). The agreement 
between these two independent approaches is encourag- 
ing, suggesting that a contour of about 125 M© pc~^ 
is a reasonable defining level for a star-forming clump 
(§2.2). Other studies have found similar thresholds for 
efficient star formation or the presence of dense cores 
(Onishi et al. 1998: Enoch et al. 2007; Johnstone et al. 
2004; Andre et al. 2010; Li et al. 1997; Lada 1992). The- 
oretical explanations for such thresholds can be found 
by considering magnetic support (Mouschovias & Spitzer 
1976) or regulation by photo-ionization (McKee 1989). 
Alternatively, this threshold may simply correspond to 
the part of the cloud that is gravitationally bound (cf. 
§4.3). 

The star formation rate density is even higher within 
clumps. Gutermuth et al. (2011) find a continuation of 
the steep increase in E(SFR) to higher ^gas in a study 
including embedded clusters. They find that E(SFR) oc 
Sgas^ up to several 100 Mq pc~^, with no evidence of a 
threshold. Near the centers of some centrally condensed 
clumps, T,gas reaches 1 gm cm~^, or about 4800 M© pc~^ 
(Wu et al. 2010a), a threshold for the formation of mas- 
sive stars suggested by theoretical analysis (Krumholz & 
McKee 2008). 

Other studies have identified possible scaling relations 
in which SFR surface densities are not simply propor- 
tional to surface densities of dense gas. One analysis 
found that clouds forming stars with masses over 10 M© 
satisfied the following relation: m{r) > 870MQ(r/pc)'^ '^^ 
(Kauffmann et al. 2010; Kauffmann & PiUai 2010), where 
m(r) is the enclosed mass as a function of radius. This 
is essentially a criterion they call "compactness" ; it can 
be thought of as requiring a central condensation, with 
p(r) oc rP, with p > 1.67. 

Dunham et al. (2011b) have compared a subset of 
clumps in the BGPS catalog (Aguirre et al. 2011; 
Rosolowsky et al. 2010) with kinematic distances to the 
Heiderman-Lada (HL) criterion for efficient star forma- 
tion and to the Kauffman-Pillai (KP) criterion for mas- 
sive star formation. About half the clumps satisfy both 
criteria. Interestingly, Dunham et al. (2011a) found that 
about half the BGPS sample contained at least one mid- 
infrared source from the GLIMPSE survey. For the 
clumps most securely identified with star formation, 70% 
to 80% satisfy the HL or KP criteria. 

All these studies consistently show that S(SFR), espe- 
cially for massive stars, is strongly localized to dense gas, 
and is much higher (for Y,gas > 100 M© pc^^) than in the 
extragalactic Schmidt relations (Kennicutt 1998b; Bigiel 
et al. 2008). This offset between the extragalactic and 
dense clump "Schmidt laws" can be straightforwardly 
understood as refiecting the mass fraction (or surface fill- 
ing factor) of dense clumps in the star-forming ISM. If 
most or all star formation takes place in dense clumps, 
but the clumps contain only a small fraction of the; total 
molecular gas, we would expect the characteristic star 
formation tdep measured for clumps to be proportion- 
ally shorter than tdep of the total cloud mass and the 
efficiency (e) to be higher. Using a threshold of G^*0 
emission to define clumps, Higuchi et al. (2009) found 
that the star formation efficiency (e) in clumps varies 
widely but averages 10%, about 4 times that in clouds 
as a whole (§4.1). If this interpretation is correct, we 



would also expect the relations between SFR and dense 
gas mass for galaxies (Gao-Solomon relation) to be sim- 
ilar for galaxies and the clumps, and they indeed appear 
to be roughly consistent (Wu et al. 2005b). 

This conclusion is sulDject to some caveats, how- 
ever. Various authors (Krumholz & Thompson 2007; 
Narayanan et al. 2008; Juneau et al. 2009) have pointed 
out that lines like those of HCN J = 1 ^ are not ther- 
malized at lower densities, so can result in a linear re- 
lation even if the underlying star formation relation is 
a local version of a non-linear KS law that extends to 
much lower densities. The SFR vs dense gas relation is 
based on a correlation between total FIR and HCN lu- 
minosities, and Lpin is likely to underestimate the SFR 
in young clusters (Urban et al. 2010; Gutermuth et al. 
2011; Krumholz & Thompson 2007), perhaps by up to 
factors of 3-30. However, Lada et al. (2012) find a sim- 
ilar continuity between the Gao-Solomon starbursts and 
SFRs measured by counting YSOs in gas above a thresh- 
old surface density. 

7. TOWARD A SYNTHESIS 
7.1. Summary: Clues from Observations 

Before we embark on an interpretation of the observed 
star formation law, it is useful to collect the main conclu- 
sions which can be drawn from the observations of star 
formation both outside and inside the Galaxy. 

Beginning on the galactic scale (> 1 kpc), we can iden- 
tify at least two and probably three distinct star forma- 
tion regimes (Table 3). Whole galaxies may lie in one 
of these regimes, but a single galaxy may include two or 
three regimes. 

7.1.1. The Low-density Regime 

The lowest-density regime (sometimes referred to as 
the sub-threshold regime), is most readily observed in 
the outer disks of spiral galaxies, but it also can be 
found in the interarm regions of some spiral galaxies and 
throughout the disks of some gas-poor galaxies. The so- 
lar neighborhood lies near the upper end of this low- 
density regime and the outer disk of the MW is clearly 
in this regime (§5.1). 

The cold gas in these regions is predominantly atomic, 
though local concentrations of molecular gas are often 
found. Star formation is highly dispersed, with young 
clusters and H II regions only observed in regions of un- 
usually high cold gas densities. The global "efficiency" of 
star formation, e' (defined in §1.2), or S(SFR)/Egos, is 
very low, and it is uncorrelated with T,gas- The steep, 
nearly vertical "Schmidt" relation seen in this regime 
(e.g., Fig. 12) mainly refiects lack of correlation over a 
region where the SFR has a much larger dynamic range 
than the local gas density; any apparent correlation is not 
physical. Recent stacking analysis of CO maps and stud- 
ies of the outer MW suggest, however, that there may be 
a strong correlation with molecular surface density (§6, 
Schruba et al. 2011). 

7.1.2. The Intermediate- density Regime 

The next, intermediate-density regime is roughly char- 
acterized by average gas surface densities of T^gas > 10 
M© pc~^, corresponding roughly to A''(H) ~ lO^^cm"^, 
or Ay ~ 1 mag, for solar metallicity. The upper limit 
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to this regime is about Sgas ^ 100 — 300 M0 pc~^, as 
discussed in the next section. The intermediate regime 
applies within the main optical radii (R25) of most gas- 
rich, late-type spiral and irregular galaxies. The "Galac- 
tic Ring" region of the Milky Way lies at the low end of 
this regime, while the CMZ may lie near the high end 

The transition from the low-density regime roughly 
corresponds to the transition between H I-dominated and 
H2-dominated ISMs. Above Egas = 10 Mq pc~^, both 
the phase balance of the ISM and the form of the star 
formation law begin to change. The intermediate range 
features an increasing filling factor of molecular clouds, 
and star formation becomes more pervasive. The SFR 
surface density is strongly and tightly correlated with 
the cold gas surface density, whether expressed in terms 
of the total (atomic plus molecular) or only the molecular 
surface density. In most massive spiral galaxies, the cold 
gas in this regime is molecular-dominated. In low-mass 
and irregular galaxies, atomic gas can dominate, though 
this is somewhat dependcint on the value of X{CO) that 
is assumed. The characteristic depletion time {tdep §1-2) 
for the interstellar gas is 1-2 Gyr (e.g., Bigiel et al. 2011). 

When expressed in terms of e' (§1.2), the star forma- 
tion rate per unit total gas mass, nearly all studies sug- 
gest that e' increases with gas surface density, with an 
exponent of 0.2 — 0.5, (the Schmidt law exponent, N—1). 
When measured against molecular mass (or surface den- 
sity), some studies suggest e' is constant, but others sug- 
gest e' increasing systematically with surface density. 

7.1.3. The High-density Regime 

The two regimes discussed above were able to repro- 
duce all of the early observations of the large-scale star 
formation law by Kennicutt (1989, 1998b) and Martin 
& Kennicutt (2001). However a number of recent ob- 
servations suggest the presence of a third regime (and 
second transition) around T,gas > 100 — 300 M© pc~^, 
into what one might call the high-density, or starburst, 
regime. Around this value of Sgasi the interpretation of 
7(C0) is likely to change (§2.4), and studies of local MW 
clouds indicate that a similar 'Egas value may correspond 
to the theoretical notion of a cluster-forming clump, in 
which the SFR is much higher than in the rest of the 
cloud (§6.4). In the most extreme starburst galaxy en- 
vironments, if standard values of X{CO) are used, the 
average surface densities of gas (virtually all molecular) 
reach 1000 to 10^ Mq pc^^, and the volume filling factor 
of clumps could reach unity (Wu et al. 2009). While the 
interpretation of molecular emission in these conditions 
warrants skepticism (§2.4, Garcia-Burillo et al. 2012), the 
highest inferred surface density also corresponds to the 
densest parts of cluster forming clumps and the theo- 
retical threshold of 1 gm cm~^ for efficient formation of 
massive stars (§4.3). 

Some recent observations strongly hint at the existence 
of such a transition. As discussed in §6, luminous and ul- 
traluminous starburst galaxies (and high-redshift SMGs) 
have characteristic ratios oi LifiJ L{CO) that are as much 
as 1-2 orders of magnitude higher than in normal galax- 
ies, implying that at some point the SFR per molecu- 
lar mass must increase dramatically. This could be ex- 
plained by a break in the slope of the Schmidt law at 



high densities, a continuous non-linear Schmidt law slope 
extending from the intermediate to high-density regimes, 
or a second mode of star formation in extreme starbursts 
with much higher "efficiency" (e'). As discussed in §6.3, 
the constancy of the molecular e' at intermediate sur- 
face densities is uncertain. Observations of CO in high- 
redshift galaxies have been interpreted in terms of just 
such a bimodal Schmidt law (Daddi et al. 2010; Genzel 
et al. 2010). This interpretation rests on the assumption 
of a bimodality in X(CO) (§6.1), and a change in e' that 
follows these changes in X{GO); this is not entirely im- 
plausible, because the same physical changes in the ISM 
environment in the densest starbursts could affect both 
the CO conversion factor and e' if an increasing fraction 
of the gas is in dense clumps. 

Both a Galactic value of X{CO) in the starbursts and 
a much lower value have their discomforting aspects. Ap- 
plying a Galactic conversion factor produces total molec- 
ular masses which often exceed dynamical mass limits for 
the regions, whereas adopting values of X(CO) which are 
factors of several lower produces gas consumption times 
as short as 10 Myr (Daddi et al. 2010; Genzel et al. 2010), 
with implications for the triggering and duty cycles of 
these massive starbursts. 

The relations between SFR and ISM properties s\im- 
marizecl above strictly refer to the (;orrelations with the 
total cold gas surface density or in some cases the total 
H I and total H2 densities; these are the relations of most 
interest for applications to galaxy modelling and cosmol- 
ogy. However the observational picture is quite differ- 
ent when we correlate the SFR with the supply of dense 
gas, as traced by HCN J = 1 — > and other dense clump 
tracers. Here there seems to be a single linear relation 
which extends across all of the SFR regimes described 
above, and which even extends to star-forming clouds in 
the Milky Way (Gao & Solomon 2004; Wu et al. 2005b; 
Gao et al. 2007; Lada et al. 2012), However, there is 
still some evidence of bimodality even if tracers of dense 
gas are used (Garcia-Burillo et al. 2012), and physically 
based models, rather than simple conversion factors, are 
needed. 

7.1.4. Clues from Studies of Molecular Clouds 

The biggest hurdle one confronts when attempting to 
interpret these observations of galaxies is the severe in- 
fluences of spatial averaging, both across the sky and 
along the line of sight. The star formation law relates 
surface densities of young stars and interstellar gas- al- 
ready smoothed along the line of sight averaged over 
linear dimensions ranging from order 100 pc to 50 kpc. 
These are 2-4 orders of magnitude larger than the sizes of 
the dense clump regions in which most stars form, and 4- 
8 orders of magnitude larger in terms of the surface areas 
being measured. This difference in scales means that the 
"surface densities" measured in extragalactic studies are 
really characterizing the filling factors of gas clumps and 
star-forming regions, rather than any measure of physical 
densities. Understanding how this "active component" of 
the star-forming ISM works is essential to even an em- 
pirical understanding of large-scale star formation, much 
less understanding its underlying physics. Detailed stud- 
ies within the MW provide a way to "zoom in" further 
than is possible for other galaxies. 

A common observed feature across all of the density 
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regimes observed in galaxies is that star formation takes 
place in molecular clouds. In the lowest-density regimes, 
clouds are rare and widely separated, but within individ- 
ual molecular clouds in the MW, which can be probed in 
detail, the star formation seems to be indistinguishable 
from that in regions of somewhat higher average surface 
density (§5.1). 

These observations also reveal that nearly all star for- 
mation within molecular clouds is highly localized, tak- 
ing place in clumps, roughly defined by Smoi > 125 
M© pc~^, or n > 10^ cm"^ (§6.4). The clumps host 
young stars, YSOs, and pre-stellar cores, the sites of in- 
dividual star formation. This scale is as close to a deter- 
ministic environment for star formation as can be found. 
Once a pre-stellar core reaches substantial central con- 
densation, about one-third of its mass will subsequently 
turn into young stars within a few Myr (Alves et al. 2007; 
Enoch et al. 2008). The relatively low global efficiencies 
(both e and e') within GMCs are largely a reflection of 
the low mass fraction in clumps and cores (§4.1). The 
star formation rate surface density in the clumps is 20- 
40 times that predicted for the mass surface density from 
the extragalactic Schmidt relation (§6.4), and this like- 
wise can be largely understood as reilecting the low mass 
fraction of molecular gas in clumps and cores. 

One might be tempted to identify the dense clumps 
within molecular clouds as a possible fourth density 
regime, in addition to the three regimes already discussed 
from observations of galaxies. However this would be 
very misleading, because so far as is currently known 
the formation of most stars in dense clumps is a com- 
mon feature of star formation across all of the ISM en- 
vironments in galaxies, extending from low-density sub- 
threshold disks to the most intense starbursts. The dense 
clump may well be the fundamental unit of massive, clus- 
tered star formation (Wu et al. 2005b). If this picture is 
correct, then the three regimes identified in galaxies and 
the order-of-magnitude increases in gas-to-star formation 
conversion rate across them must reflect changes in the 
fraction of the ISM that is converted to dense clumps, 
approaching 100% in the densest and most intense star- 
bursts. 

7.2. Some Speculations 

The observations summarized in this review have stim- 
ulated a rich literature of theoretical ideas, models, and 
simulations aimed at explaining the observed star for- 
mation relations and constructing a coherent picture of 
galactic-scale star formation. Unfortunately we have nei- 
ther the space nor the expertise to review that large body 
of theoretical work here. A review of many of the theo- 
retical ideas can be found in McKee & Ostriker (2007), 
and an informative siimmary of the main models in the 
literature up to 2008 can be found in Leroy et al. (2008). 
Here we offer some speculations aimed toward explain- 
ing the observations and connecting the extragalactic and 
MW studies, and identifying directions for future work. 
Many of these speculations reflect ideas being discussed 
in the literature, and we make no claims for originality 
in the underlying concepts. 

As mentioned in the introduction, the formation of 
stars represents the endpoint of a chain of physical pro- 
cesses that begins with cooling and infall of gas from 
the intergalactic medium onto disks, followed by the for- 



mation of a cool atomic phase, contraction to gravita- 
tionally bound clouds, the formation of molecules and 
molecular clouds, the formation of dense clumps within 
those molecular clouds, and ultimately the formation of 
pre-stellar cores, stars, and star clusters. Although the 
beginning and end points of this process are relatively 
clear, the sequence of intermediate steps, in particular 
the respective roles of forming cool atomic gas, molecular 
gas, and bound clouds is unclear, and it is possible that 
different processes dominate in different galactic envi- 
ronments. By the same token, a variety of astrophysical 
time scales may be relevant: e.g., the free-fall time of the 
gas within clouds, the crossing time for a cloud, the free- 
fall time of the gas layer, or the dynamical time scales for 
the disc and spiral arm passages. Any of these timescales 
may be invoked for setting the time scale for star forma- 
tion and the form of the star formation law. However 
one can construct two scenarios to explain the observa- 
tions outlined in §7.1, which illustrate the boundaries 
of fully locally-driven versus globally-driven approaches. 
In some sense, these approaches are complementary, but 
they need to be brought together. 

The first approach, which might be called a bottom- 
up picture, assumes that star formation is controlled lo- 
cally within molecular clouds, (e.g., Krumholz & McKee 
2005), building on what we observe in well-studied lo- 
cal regions of star formation. In this picture, one can 
identify three distinct regimes (§7.1), and associate the 
transitions between them to the crossing of two physi- 
cally significant thresholds, the threshold for conversion 
of atomic to molecular gas, and the threshold for efficient 
star formation in a molecular cloud, identified with the 
theoretical notion of a clump (§2.2) . In the purest form of 
this picture the SFR is driven completely by the amount 
and structure of the molecular gas. Current evidence 
from studies of nearby clouds indicates that the star for- 
mation rate scales linearly with the amount of dense gas 
in clumps (Lada et al. 2010), and this relation extends to 
starburst galaxies if the HCN emission is used as a proxy 
for the dense gas (Wu et al. 2005b; Lada et al. 2012). In 
this picture, the non-linear slope (N 1.5) of the global 
Schmidt relation would arise from either a decrease in 
the characteristic timescale (Krumholz et al. 2012) or 
by an increase in the fraction of gas above the clump 
threshold (fdense), from its typical value in local clouds 

[fdense ~ 0.1 (Lada et al. 2012)] with fdense CX Egas" ''. 

The second approach, which we could call a top-down 
picture, assumes that star formation is largely controlled 
by global dynamical phenomena, such as disk instabil- 
ities (e.g.. Silk 1997), and the dynamical timescales in 
the parent galaxy. In this picture the transition be- 
tween the low-density and higher-density SFR regimes 
is mainly driven by gravitational instabilities in the disk 
rather than by cooling or molecular formation thresh- 
olds, and the; non-linear increase in the SFR relative 
to gas density above this threshold reflects shorter self- 
gravitational timescales at higher density or the shorter 
dynamical timescales. In this picture there is no particu- 
lar physical distinction between the intermediate-density 
and high-density regimes; in principle the same dynam- 
ical processes can regulate the SFR continuously across 
this wide surface density regime. Likewise it is the total 
surface density of gas, whether it be atomic or molecular, 
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that drives the SFR. The asymptotic form of this pictme 
is a self-regulated star formation model: the disk adjusts 
to an equilibrium in which feedback from massive star 
formation acts to balance the hydrostatic pressure of the 
disk or to produce an equilibrium porosity of the ISM 
(e.g., Cox 1981; Dopita 1985; Silk 1997; Ostriker et al. 
2010). 

When comparing these pictures to the observations, 
each has its particular set of attractions and challenges. 
The bottom-up picture has the attraction of simplic- 
ity, associating nearly all of the relevant SFR physics 
with the formation of molecular gas and molecular cloud 
clumps. It naturally fits with a wide range of observa- 
tions including the tight correlation of of the SFR and 
molecular gas surface densities (e.g., Schruba et al. 2011 
and references therein) , the concentration of star forma- 
tion within clouds in regions of dense gas (§4), and the 
observation of a linear relation between the dense gas 
traced by HCN emission and the total SFR in galaxies 
(§6). If the resolved star formation relation at interme- 
diate surface densities is linear (§6, Bigiel et al. 2008), 
f dense would be Constant in that regime, while increasing 
monotonically with T,gas in the higher-density starburst 
regime with the transition occuring where "Sgas derived 
from CO is similar to the threshold for dense gas. This 
picture does not explain why a particiilar galaxy or part 
of a galaxy lies in one of these regimes, a question per- 
haps best answered by the top-down picture. 

Aspects of the top-down picture date back to early 
dynamical models of the ISM and the first observations 
of star formation thresholds in disks (e.g.. Quirk & Tins- 
ley 1973; Larson 1987; Zasov & Stmakov 1988; Kennicutt 
1989; Ehncgrccn 1991; Silk 1997). There is some observa- 
tional evidence for associating the observed low-density 
thresholds in disks with gravitational instabilities in the 
disk (e.g., Q instabilities), rather than with atomic or 
molecular phase transitions (e.g., Kennicutt 1989; Martin 
& Kennicutt 2001), but recent observations and theoret- 
ical analyses have raised questions about this interpreta- 
tion (e.g., Schaye 2004; Leroy et al. 2008). At higher sur- 
face densities the global relation between E(SFR) and the 
ratio of gas density to local dynamical time {^gas/Tdyn) 
shows a correlation that is nearly as tight as the con- 
ventional Schmidt law (Kennicutt 1998b), and it also re- 
moves the double sequence of disks and starbursts that 
results if X{GO) is systematically lower in the starbursts 
(§6, Daddi et al. 2010; Genzel et al. 2010). In this pic- 
ture, the higher star formation rate for a given gas surface 
density in mergers is caused by the compaction of the gas 
and the resulting shorter rotation period. However, the 
efficiency per orbital period is not clearly explained in 
this picture. Theories of feedback-regulated star forma- 
tion show promise in explaining the low efficiency per 
orbit in normal galaxies (Kim et al. 2011), and they may 
also explain the less effective role of negative feedback in 
merger-driven starbursts (e.g., Ostriker & Shetty 2011). 
This model does not extend to cloud-level star formation. 
If feedback from massive stars is fundamental, the simi- 
larity between the efficiency in local clouds forming only 
low mass stars and regions with strong feedback from 
massive stars (cf. Evans et al. 2009 and Murray 2011) 
remains a mystery. 

Both of these scenarios can claim some successes, but 
neither provides a complete explanation that extends 



over all scales and environments. As but one example, 
the large changes in the present-day SFRs and past star 
formation histories of galaxies as functions of galaxy mass 
and type may well be dictated mainly by external in- 
fluences such as the accretion history of cold gas from 
the cosmic web and intergalactic medium. Neither the 
bottom-up or top-down pictures as articulated above in- 
corporate these important physical processes. On smaller 
scales a complete model for star formation may combine 
features of both scenarios. The dynamical picture is quite 
attractive on the scales of galaxies, and it would be in- 
teresting to extend it to smaller scales. In doing so, the 
question of what to use for Tdyn arises. The galaxy ro- 
tation period cannot control star formation in individual 
clouds, so a more local dynamical time is needed. One 
option is the cloud or clump crossing time (essentially 
the size over the velocity dispersion) (Elmegreen 2000), 
which may be particularly relevant in regions of triggered 
star formation. Since clumps arc the; star forming units, 
their crossing times may be the more relevant quantities. 

A popular option is the free-fall time (e.g., Krumholz 
& McKee 2005). Since tff oc p~^'^, a volumetric star 
formation law, psFR p/tff oc p^'^, where p is the gas 
density, is a tempting rule. With this rule, the roughly 
1.5 power of the KS relation appears to be explained if 
we ignore the difference between volume and surface den- 
sity. Indeed, Krumholz et al. (2012) argue that such a 
volumetric law reproduces observations from the scale of 
nearby clouds to starburst galaxies. However, the defi- 
nition of t J f changes for compact starbursts (see eq. 9 
in their paper), essentially at the boundary between the 
intermediate-density and high-density regimes discussed 
in §7.1. Including the atomic- molecular threshold, as 
treated in Krumholz et al. (2009b), clarifies that this 
"scale-free" picture implicitly recognizes the same three 
regimes discussed in §7.1. 

Theories that rely on tff face two problems. First, 
no evidence has been found in well-studied molecular 
clouds for collapse at tff (§4.3, e.g., Zuckerman & Evans 
1974). To match observations, an "efficiency" of about 
0.01 must be inserted. Why the star formation density 
should remain proportional to tff, while being slowed 
by a factor of 100, is a challenging theoretical question. 
Studies of the role of turbulence show some promise in 
explaining this paradox (Krumholz & McKee 2005; Hen- 
nebelle & Chabrier 2011), and magnetic fields may yet 
play a role. 

The second problem is how to calculate a relevant tff 
in a cloud, much less a galaxy, with variations in p by 
many orders of magnitude, and much of the gas unlikely 
to be gravitationally bound. For example, Krumholz 
et al. (2012) calculate tff from the mean density of the 
whole cloud {tff oc l/{p)). A computation of {tff) from 
{1/ p) would emphasize the densest gas, where star forma- 
tion is observed to occur. In a model that accounts for 
this, Hennebelle & Chabrier (2011) reproduce at some 
level the observations of surface density thresholds for 
efficient star formation observed in nearby clouds (Lada 
et al. 2010; Heiderman et al. 2010). An alternative view 
to the models involving tff and emphasizing instead the 
critical role of the dense gas threshold can be found in 
Lada et al. (2012). Because both models claim to apply 
to cloud-level star formation, observations of MW clouds 
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should be able to test them. 

7.3. Future Prospects 

We hope that this review has conveyed the tremendous 
progress over the last decade in understanding the sys- 
tematic behavior of star formation, both in the MW and 
in other galaxies. As often happens when major obser- 
vational advances arc made, observational pictures that 
once seemed simple and certain have proven to be more 
complex and uncertain. We have attempted to inject a 
dose of skepticism about some well-accepted truisms and 
to highlight important questions where even the obser- 
vations do not present a completely consistent picture. 

As we return to the key questions outstanding in this 
subject (§1.3), a few clear themes emerge, many of which 
lie at the cusp between MW and extragalactic observa- 
tions, and between theory and simulation on the sub- 
cloud scale on the one hand, and the galactic and cos- 
mological scales on the other. The recent observations 
within the MW of a near-universal association of stars 
with dense molecular clumps (§6.4) offers the potential 
key of a fundamental sub-unit of high-efficiency star for- 
mation in all galactic environments, from the low-density 
and quiescent environments of outer disks and dwarf 
galaxies to the most intense starbursts. However this hy- 
pothesis needs to be validated obscrvationally in a wider 
range of environments. Restated from another perspec- 
tive, we need much better information on the structure 
(physical structure, substructure) and dynamics of star- 
forming clouds across the full range of star-forming en- 
vironments found in galaxies today. 

Fortunately we are on the brink of major progress on 
multiple fronts. One approach is to assemble more com- 
plete and "zoomed-out" views of the MW, while pre- 
serving the unparalleled spatial resolution and sensitiv- 
ity of MW observations to fully characterise the statis- 
tical trends in cloud structure, kinematics, mass spec- 
tra, and associated star formation for complete, unbi- 
ased, and physically diverse samples. In the near-term, 
Herschel surveys will deliver images of nearby clouds, the 
plane of the MW, and many galaxies in bands from 60 to 
500 /im. When the MW plane data are combined with 
higher resolution surveys of the MW at 0.87 to 1.1 mm 
from ground-based telescopes, spectroscopic follow-up, 
and improved distances from ongoing VLBA studies, we 
will have a much improved and nuanced picture of the gas 
in the MW, which can then provide a more useful tem- 
plate for understanding similar galaxies. Recent surveys 
have doubled the number of known H II regions in the 
MW (Anderson et al. 2011), allowing study of a wider 
range of star formation outcomes. Tests of the limita- 
tions on star formation rate tracers used in extragalactic 
work will be possible, as will comparisons to KS relations 
on a variety of spatial scales. The Magellanic Clouds also 
offer great potential for extending this approach to two 
sets of environments with different metallicities, ISM en- 
vironments, and star formation properties (exploiting for 
example the 30 Doradus region). 

This expanded information from MW surveys will then 
need to be combined with more sensitive and "zoomcd- 
in" observations of other galaxies. For such studies the 
weakest link currently is the knowledge of the molecular 
gas. We can trace H I and SFRs to much lower lev- 
els than we can detect CO, even with stacking of the 



CO maps. When fully commissioned, ALMA will have a 
best resolution of about 13 mas at 300 GHz (0.5 pc at 
the distance of M51), and about 8 times the total col- 
lecting area of any existing millimeter facility, allowing 
us to trace molecular emission to deeper levels and/or to 
obtain resolution at least 10 times better than the best 
current studies. It may be possible to study dense struc- 
tures within molecular clouds in other galaxies for com- 
parison to dense clumps in the MW clouds. To be most 
effective, these programs will need to expand beyond the 
typical surveys in one or two CO rotational transitions to 
include the high-density molecular tracers and ideally a 
ladder of tracers with increasing excitation and/or criti- 
cal density. Continuum observations will gain even more 
because of large bandwidths and atmospheric stability. 
Maps of mm-wave dust emission may become the pre- 
ferred method to trace the gas in other galaxies, as they 
are becoming already in the MW, bypassing the issues of 
A" (CO). Observations of radio recombination lines may 
also provide alternative tracers of star formation rate in 
highly obscured regions. It is difficult to exaggerate the 
potential transformational power of ALMA for this sub- 
ject. 

For molecular line observations of nearby galaxies, 
large single dishes (the Nobeyama 45-m, the IRAM 30-m, 
and the future CCAT), along with the smaller millime- 
ter arrays, IRAM PdB and CARMA, will provide com- 
plementary characterization of the molecular and dense 
gas on larger scales than can be efficiently surveyed with 
Alma, a substantial expansion of the IRAM interferom- 
eter (NOEMA) will provide complementary advances in 
sensitivity in the northern hemisphere. This subject will 
also advance with a substantial expansion of H I map- 
ping (by e.g., the EVLA) to expand the range of envi- 
ronments probed and ideally to attain spatial resolutions 
comparable to the rapidly improving molecular line ob- 
servations. Although it is likely that star formation on 
the local scale concentrates in molecular-dominated re- 
gions, the formation of this molecular gas from atomic 
gas remains a critical (and possibly controlling) step in 
the entire chain that leads from gas to stars. 

Future opportunities also await for improving our mea- 
surements of SFRs in galaxies, though the truly transfor- 
mational phase of that subject may already be passing 
with the end of the ISO, Spitzer, GALEX, and Herschel 
eras. HST continues to break new ground, especially 
in direct mapping of young stars and their ages via re- 
solved color-magnitude diagrams. Herschel surveys of 
nearby galaxies will provide far-infrared data, comple- 
mentary to the other wavelength regions that trace star 
formation (§3), and multiband imaging with the EVLA 
will allow better separation of non-thermal from the free- 
free continuum emission, which directly traces the ion- 
ization rate and massive SFR. SOFIA offers the opportu- 
nity for mapping the brightest regions of star formation 
with higher spatial and spectral resolution than Spitzer. 
Further into the future, major potential lies with JWST, 
which will vastly improve our capabilities for tracing star 
formation via infrared radiation, both on smaller scales 
in nearby galaxies and in very distant galaxies. SPICA 
may provide complementary capability at longer wave- 
lengths. As highlighted in §1.3, key unknowns in this 
subject include robust constraints on the ages and life- 
times of star- forming clouds; observations of resolved star 
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clusters, both in the MW and nearby galaxies, offer po- 
tential for considerable inroads in this problem. 

A number of the questions we have listed also can be 
attacked with groundbased OIR observations. As dis- 
cussed in §3.9, despite the availability now of spatially- 
resolved multi-wavelength observations of nearby galax- 
ies, we still do not have an absolutely reliable way to 
measure dust-corrected SFRs, especially from short-lived 
tracers that are critical for exploring the local Schmidt 
law. For observations of normal galaxies where local at- 
tenuations are modest, integral- field mapping of galaxies 
offers the means to produce high-quality Ha maps cor- 
rected for dust using the Balmer decrement (e.g., Blanc 
et al. 2009, 2010; Sanchez et al. 2012). A "gold stan- 
dard" for such measurements also is available in the hy- 
drogen recombination lines of the Paschen and Brackett 
series in the near-infrared. The emissivities of these lines 
are directly related to ionizing luminosities in the same 
way as are the more widely applied Balmer lines (albeit 
with somewhat stronger density and temperature depen- 
dences), and these lines suffer much lower dust attenu- 
ation (readily calibrated by comparing to Ha or other 
shorter wavelength lines). These lines are much fainter, 
however, and are subject to strong interference from tel- 
luric OH emission (and thermal emission at longer wave- 
lengths). Continuing advances in near-infrared detectors 
have now brought many of these lines into the accessible 
range, and soon we should begin to see large-scale sur- 
veys that will provide high-resolution emission maps and 
robust "SFR maps" , at least on spatial scales larger than 
individual H 11 regions. These maps in turn can be used 
to test and hopefully recalibrate other dust-free tracers 
such as combinations of infrared dust emission with UV 
and optical emission-line maps. Looking further ahead, 
spectral imaging in the near-infrared (and with Ha in 
the visible) with massively parallel integral-field spec- 
trometers will provide more precise and full views of star 
formation. 

The current, pioneering studies of star formation at 
2; ~ 1 — 3 will become much easier with new instru- 
mentation. Advances in submillimeter array technology 
(e.g., SCUBA-2 on the JCMT) and larger dishes at very 
high altitude, such as CCAT, will allow deep and wide 
searches for dust continuum emission from distant galax- 
ies. Many of the biases in current samples can be alle- 
viated and a more complete picture of star formation 
through cosmic time can be constructed. Huge surveys 
of Lya emitting galaxies will be undertaken to constrain 
dark energy, providing as a by-product, nearly a million 
star- forming galaxies at 1.9 < z < 3.5 and a large num- 
ber of [OH] emitting galaxies for z < 0.5 (e.g.. Hill et al. 
2008). 

Since the theoretical side of the subject lies outside 
the scope of this review, we comment only briefly on 
this area. Numerical simulations are poised to make 
major contributions to the subject over the next sev- 
eral years. Simulations with higher resolution and more 
sophisticated treatments of the heating, cooling, phase 
balance, and feedback are providing deeper insights into 



the physical nature of the star formation scaling laws 
and thresholds and the life cycles of molecular clouds 
(e.g., Robertson & Kravtsov 2008; Dobbs 2008; Dobbs & 
Pringie 2010; Tasker & Tan 2009; Tasker 2011; Brooks 
et al. 2011). Inclusion of realistic stellar feedback may 
obviate the need for artificial constraints on efficiency 
(e.g., Hopkins et al. 2011). 

Analytical models which incorporate the wide range of 
relevant physical processes on the scales of both molec- 
ular clouds and galactic disks are also leading to deeper 
insights into the triggering and regulation of star forma- 
tion on the galactic scale (e.g., Ostriker et al. 2010, and 
papers cited in §7.2). Exploring different prescriptions 
for cloud-scale star formation in galaxy evolution mod- 
els will illuminate the effects of the cloud-scale prescrip- 
tion on galaxy-scale evolution. At the "micro-scale'' of 
molecular clouds, inclusion of radiative and mechanical 
feedback from star formation is producing more realistic 
models, and perhaps the relative importance of feeding 
from the local core and from the greater clump may be 
quantified (e.g., references in §4.3). 

We conclude with a hope and a prediction. From the 
start, this review was designed to exchange ideas be- 
tween those studying star formation in the Milky Way 
and those studying star formation in other galaxies. The 
authors have benefited immensely from this exchange, 
and we hope for more of this cross-talk between our two 
communities: people studying star formation within the 
Milky Way can put their results in a larger context; 
those studying other galaxies can appreciate that our 
home galaxy offers unique advantages for understand- 
ing galaxies. A common theme across both arenas has 
been the rapid pace of advance, both in obscirvations and 
theory. When the next review of this subject is written, 
we predict that it will focus on observational tests of de- 
tailed physical theories and simulations rather than on 
empirical star formation laws. We conclude with heart- 
felt thanks to all of those who are working toward this 
common goal. 
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TABLE 1 
Star Formation Rate Calibrations 



Band 


Age Range (Myr)'' 


Units 




M«/M*(K98) 


References 


FUV 


- 10 - 100 


ergss^"'" (uLv) 


43.35 


0.63 


1, 2 


NUV 


- 10 - 200 


crgss^^ (vLv) 


43.17 


0.64 


1, 2 


Ho 


0-3-10 


ergs s^-"- 


41.27 


0.68 


1, 2 


TIR 


- 5 - 100*^ 


crgss^-*- (3— 1100 ^m) 


43.41 


0.86 


1, 2 


24 fj.m 


- 5 - 100*^ 


crgss^-"- (uL^) 


42.69 




3 


70 ^J,m 


- 5 - 100*^ 


crgss^-"- (uLv) 


43.23 




4 


1.4 GHz 


- 100 : 


ergss"'^ Hz~^ 


28.20 




1 


2-10 keV 


- 100 : 


ergss"'^ 


39.77 


0.86 


5 



Note. — References: (1) Murphy et al. (2011); (2) Hao et al. (2011); (3) Rieke et al. (2009); (4) Calzetti et al. (2010a); 
(5) Ranalli et al. (2003) 

'^Second number gives mean age of stellar population contributing to emission, third number gives age below which 90% of emission is 
contributed. 

''Numbers are sensitive to star formation history, those given are for continuous star formation over 0-100 Myr. For more quiescent regions 
(e.g., disks of normal galaxies) the maximum age will be considerably longer. 



TABLE 2 

Multi-Wavelength Dust-Corrections 



Composite Tracer 


Reference 


L{FUV)^orr 


= L[FUV)obs + 0A6L(TIR) 


1 


L{FUV)corr 


= L(FUV)obs + 3.89L(25Mm) 


1 


L{FUV)corr 


= L(FUV)obs + 7.2 X 10i''L(1.4GHz)'' 


1 


L{NUV)corr 


= L{NUV)obs + 0.27L(T/i?) 


1 


L{NUV)rorr 


= L{NVV)„b, + 2.26L(25//ni) 


1 


L{Ai:\-),„r, 


= L(A"(;i ")„(,,, + 1.2 X lO-^-iL(l.lGliz)" 


1 


L{Ha)corr = 


L{Ha)obs + 0.0024L(TIR) 


2 


L(Ha)corr = 


L(Ha)obs + 0.020L(25Mm) 


2 


L{Ha)corr = 


L{Ha)obs +0.011L(8^/in) 


2 


L(Ho!.)corr = 


L{Ha)obs +0.39 x 10i*L(1.4GHz)'' 


2 



Note. — References: (1) Hao et al. (2011); (2) Kennicutt et al. (2009) 
''Radio luminosity in units of ergss"' Hz"' 



TABLE 3 
Star Formation Regimes 



Name 


^gas 

(Mq pc-2) 


Gas Properties 


Star Formation 


Examples 


Low 


< 10 


Mostly atomic 


low, sparse 


outer disks 


Density 








early type galaxies 
LSB galaxies 
dwarf galaxies 
solar nbd in MW 


Intermediate 


10 ^dense 


atomic-^ molecular 


moderate 


normal disks 


Density 








inner MW 
molecular clouds 


High 


^ ^dense 


molecular— >■ dense 


high, concentrated 


some nuclear regions 


Density 








starbursts 
mergers 

clumps and cores 



Note. — The dividing line between Intermediate and High-density regimes (^dense) ranges from 100 to 300 Mq pc 
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Fig. 1. — Compilation of the conversion factor (X(CO)) from the CO J = 1 — >■ integrated intensity [Ico (K km s ^)] or luminosity [L'qq 
(K km pc^)] to H2 column density (left vertical scale) and total (H2 and He) gas mass (right vertical scale), derived in various Galactic 
and extragalactic targets. Blue circles denote measurements in the disk and center of the Milky Way, based on various virial, extinction, and 
isotopomeric analyses. Crossed green squares denote measurements in starbursts and (U)LIRGs, mainly based on dynamical constraints. 
Filled triangles denote conversion factors as a function of decreasing metallicity (vertical arrow) from (bottom) to 8.2 (top), derived mainly 
from global (large scale) dust mass measurements in nearby galaxies and dwarfs by several groups. In contrast, red filled squares mark 
X-factor measurements toward individual clouds, over the same range in metallicity. Taken from Tacconi et al. (2008); references to the 
original work are given there. Reproduced by permission of the AAS. 
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NGC 6946 (6.8 Mpc, SABcd) 
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Fig. 2. a montage of infrared images of NGC 6946 from Spitzer (SINGS) and Herschel (KINGFISH). Top panels: SpitzerlRAC images 
at 3.6 fini and 8.0 fini, and MIPS image at 24 /im. The emission at these wavelengths is dominated by stars, small PAH dust grains, and 
small dust grains heated by intense radiation fields, respectively. Middle panels: Herschel PACS images at 70 ^m, 100 /im, and 160 /im, 
processed with the Scanamorphos map making package. Note the excellent spatial resolution despite the longer wavelengths, and the 
progressive increase in contributions from diffuse dust emission ("cirrus") with increasing wavelength. Bottom panels: Herschel SPIRE 
images at 250 /um, 350 //m, and 500 /um. These bands tra<;e increasingly cooler components of the main thermal dust emission, with possible 
additional contributions from "submillimeter excess" emission at the longest SPIRE wavelengths. FWHM beam sizes for the respective 
Herschel bands are shown in the lower left corner of each panel. This figure originally appeared in the Publications of the Astronomical 
Society of the Pacific (Kennicutt et al. 2011). Copyright 2011, Astronomical Society of the Pacific; reproduced with permission of the 
Editors. 
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Fig. 3. — Top: relation between observed 24 /im IR luminosity and dust-corrected Ha luminosity for nearby galaxies. The dust corrections 
were derived from the absorption-corrected Ha/H/3 ratios in optical spectra. The dotted line shows a linear relation for comparison, while 
the other lines show published fits to other samples of galaxies. Bottom: linear combination of (uncorrected) Ho and 24 fim luminosities 
compared to the same Balmer-corrected Ha luminosities. Note the tightness and linearity of the relation over nearly the entire luminosity 
range. Taken from Kennicutt et al. (2009); reproduced by permission of the AAS. 
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Fig. 4. — Example of the strong concentration of star formation in regions of high extinction, or mass surface density in the Perseus 
molecular cloud. The gray-scale with black contours is the extinction map ranging from 2 to 29 mag in intervals of 4.5 mag. The yellow 
filled circles are Flat SEE) sources and the red filled circles are Class I sources. Sources with an open star were not detected in HCO+ 
J = 3 — f 2 emission and are either older sources that may have moved from their birthplace or background galaxies. Essentially all truly 
young objects lie within contours of Ay > 8 mag. Taken from Heiderman et al. (2010); reproduced by permission of the AAS. 
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Fig. 5. — Sketch of approximately how the Galaxy is likely to appear viewed face-on, with Galactic coordinates overlaid and the locations 
of spiral arms and the Sun indicated. This sketch was originally made by Robert Hurt of the Spitzer Science Center in consultation with 
Robert Benjamin at the University of Wisconsin- Whitewater. The image is based on data obtained from the literature at radio, infrared, 
and visible wavelengths. As viewed from a great distance our Galaxy would appear to be a grand-design two-armed barred spiral with 
several secondary arms: the main arms being the Scutum-Centaurus and Perseus arms and the secondary arms being Sagittarius, the outer 
arm, and the 3 kpc expanding arm. Adapted from a figure in Churchwell et al. (2009) by R. Benjamin. The original figure appeared in the 
Publications of the Astronomical Society of the Pacific, Copyright 2009, Astronomical Society of the Pacific; reproduced with permission 
of the Editors. 
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Fig. 6. — High-resolution maps of star formation tracers and cold gas components in NGC 6946. Top left: Spitzcr 24 fim dust emission 
from the SINGS/KINGFISH project (Kennicutt et al. 2011); Top right; Ha emission from Knapen et al. (2004); Bottom left; H I emission 
from the VLA THINGS survey (Walter et al. 2008); Bottom right: CO J = 2 1 emission from the HERACLES survey (Leroy et al. 
2009). Note that the H I map extends over a much wider area than the CO, Ha, and 24 /nm observations. 
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Fig. 7. — Top: The radial distribution of surface densities of atomic gas, molecular gas, and star formation rate for the MW. The atomic 
data were supplied by P. Kaberla, (Kalberla & Dedes 2008), but we corrected for helium. The discontinuity around 8.5 kpc reflects the 
inability to model H I near the solar circle. The molecular data are taken from Table 1 in Nakanishi & Sofue (2006), but scaled up to be 
consistent with the bottom panel, using X(CO) = 2.0 X 10^" and multiplying by 1.36 to include helium. The molecular surface density 
in the innermost bin is highly uncertain and may be overestimated by using a constant value of X{CO) (§2.4). The S(SFR) is read from 
Figure 8 in Misiriotis et al. (2006), but it traces back, through a complex history, to Guesten & Mezger (1982), based on radio continuum 
emission from H II regions. That result was in turn scaled to the latest estimate of total star formation rate in the MW of 1.9 M© yr~^ 
(Chomiuk & Povich 2011). Separate points for ^mol B.nd S(SFR) are plotted for the CMZ {Rgai < 250 pc). An uncertainty of 0.3 was 
assigned arbitrarily to the values of log S(SFR). Bottom: Same as the top panel, but for NGC 6946. It is based on a figure in Schruba 
et al. (2011) but modified by A. Schruba to show radius in kpc for ready comparison to the plot for the MW. Reproduced by permission 
of the AAS. All gas distributions for both galaxies include a correction for helium. 
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log M* 

Fig. 8. — Relation between specific star formation rate (Af*/A/, ) and galaxy mass for galaxies in the SDSS spectroscopic sample, and 
with SFRs measured from GALEX ultraviolet luminosities. Gray contours indicate the ( 1 /Vmax- weighted) number distribution of galaxies 
in this bivariate plane. The blue solid line shows the fit to the star-forming sequence fit, and the red line shows approximate position of 
non-star- forming red sequence on this diagram. The dotted line shows the locus of constant SFR = 1 Mq yr~^. Figure adapted from 
Schiminovich et al. (2007). Reproduced by permission of the AAS. 
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Fig. 9. — Distribution of integrated star formation properties of galaxies in the local Universe. Each point represents an individual galaxy 
or starburst region, with the average SFR per unit area (SFR intensity) plotted as a function of the absolute SFR. The SFR intensities 
are averaged over the area of the main star-forming region rather than the photometric area of the disk. Diagonal lines show loci of 
constant star-formation radii, from 0.1 kpc (top) to 10 kpc (bottom). Several sub-populations of galaxies are shown: normal disk and 
irregular galaxies measured in Ha and corrected for dust attenuation (solid black points) from the surveys of Gavazzi et al. (2003), James 
et al. (2004), Hameed & Devereux (2005), and Kennicutt et al. (2008); luminous star-formation dominated infrared galaxies (LIRGs) and 
ultraluminous infrared galaxies (ULIRGs), plotted as red points, from Dopita et al. (2002) and Scoville et al. (2000); blue compact starburst 
galaxies measured in Ha (blue points) from Perez-Gonzalez et al. (2003) and Gil de Paz et al. (2005); and circumnuclear star- forming 
rings in local barred galaxies measured in Paa (green points) as compiled by Kormendy & Kennicutt (2004). The MW (magenta square) 
lies near the high end of normal spirals in total star formation rate, but it is near the average in S(SFR) (§5.1). The relative numbers of 
galaxies plotted does not reflect their relative space densities. In particular the brightest starburst galaxies are extremely rare. 
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log v [Mo/yr ] 

Fig. 10. — Volume-corrected SFR distribution function for disk and irregular galaxies in the local Universe, as derived from a combination 
of large flux-limited ultraviolet and infrared catalogs and a multi-wavelength survey of the local 11 Mpc volume. Vertical lines indicate 
uncertainties due to finite sampling; these are especially important for low SFRs, where the statistics are dominated by dwarf galaxies, 
which can only be observed over small volumes. Statistics on gas-poor (elliptical, dwarf spheroidal) galaxies are not included in this study. 
The blue line shows a maximum-likelihood Schechter function fit, as described in the text. Figure is original, but similar to one in Bothwell 
et al. (2011). 
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Fig. 11. — Upper: Relationship between the disk-averaged surface densities of star formation and gas (atomic and molecular) for different 
classes of star- forming galaxies. Each point represents an individual galaxy, with the SFRs and gas masses normalized to the radius of the 
main star-forming disk. The line shows the original N = 1.4 fit from Kennicutt (1998b) superimposed on the data. Most of the galaxies 
form a tight relation (exceptions discussed below), and with the improved dataset, even the normal galaxies (shown with black points) 
follow a well-defined Schmidt law on their own. The dispersion of the normal galaxies overall from the average relation (±0.30 dex rms) 
is considerably higher than can be attributed to observational uncertainties, which suggests that much of the dispersion is physical. The 
Milky Way (magenta square) fits well on the main trend seen for other nearby normal galaxies. The purple crosses show data for low surface 
brightness galaxies. The sample of galaxies has been enlarged from that studied in Kennicutt (1998b), with many improved measurements 
as described in the text. Lower: Corresponding relation between the total (absolute) SFR and the mass of dense molecular gas as traced 
in HCN. The line is a linear fit, which contrasts with the non-linear fit in the upper panel. Figure adapted from Gao & Solomon (2004). 
Reproduced by permission of the A AS. 
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Fig. 12. — Relation between SFR surface densities and total (atomic and molecular) gas surface densities for various sets of measurements, 
from Bigiel et al. (2008). Colored contours show the distribution of values from measurements of sub-regions of SINGS galaxies. Overplotted 
as black dots are data from measurements in individual apertures in M51 (Kennicutt et al. 2007). Data from radial profiles from M51 
(Schuster et al. 2007), NGC 4736, and NGC 5055 (Wong k. Blitz 2002), and NGC 6946 (Crosthwaite & Turner 2007) are shown as black 
filled circles. The disk-averaged measurements from 61 normal spiral galaxies (filled gray stars) and 36 starburst galaxies (triangles) from 
Kennicutt (1998b) are also shown. The black filled diamonds show global measurements from 20 low surface brightness galaxies (Wyder 
et al. 2009). In all cases, calibrations of IMF, X(CO), etc. were placed on a common scale. The three lines extending from lower left to 
upper right show lines of constant global star formation efficiency. The two vertical lines denote regimes which correspond roughly to those 
discussed in §6 of this review. Figure taken from Bigiel et al. (2008). Reproduced by permission of the AAS. 
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Fig. 13. — Plot of the ratio of the total YSO content of a cloud to the total cloud mass vs. total cloud mass. This is equivalent to a 
measure of the star formation efficiency as a function of cloud mass for the local sample. It is also equivalent to the measure of the SFR 
per unit cloud mass as a function of the cloud mass. The plot shows large variations in the efficiency and thus the SFR per unit mass for 
the local cloud sample. Taken from Lada et al. (2010). Reproduced by permission of the AAS. 
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Fig. 14. — Taken from Lada et al. (2010). Relation between N(YSOs), the number of YSOs in a cloud, and Mo.g, the integrated cloud 
mass above the threshold extinction of = 0-8 mag. For these clouds, the SFR is directly proportional to N(YSOs), and thus this graph 
also represents the relation between the SFR and the mass of highly extincted and dense cloud material. A line representing the best-fit 
linear relation is also plotted for comparison. There appears to be a strong linear correlation between N(YSOs) (or SFR) and Mo.g, the 
cloud mass at high extinction and density. Taken from Lada et al. (2010). Reproduced by permission of the AAS. 
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Fig. 15. — Comparison of Galactic total c2d and GB clouds, YSOs, and massive clumps to extragalactic relations. SFR and gas surfaces 
densities for the total c2d and GB clouds (cyan squares), c2d Class I and Flat SED YSOs (green and magenta stars), and Lm > lO'*'^ Lq 
massive clumps (yellow diamonds) are shown. The range of gas surface densities for the spirals and circumnuclear starburst galaxies in the 
Kennicutt (1998b) sample is denoted by the gray horizontal lines. The gray shaded region denotes the range for = 129 ± 14 Mq pc"'^. 
Taken from Heiderman et al. (2010). Reproduced by permission of the AAS. 



